#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+9;
int n,m,ans;
int qpow(int a,int b){
	int c=1;
	for(;b;b>>=1){
		if(b&1)c=1ll*a*c%mod;
		a=1ll*a*a%mod;
	}
	return c;
}
int av[10100000],jc[10010000],ij[10100000];
int main(){
	scanf("%d%d",&n,&m);
	jc[0]=1;
	for(int i=1;i<=m;i++)
	av[i]=(n-1)*i+1,jc[i]=1ll*jc[i-1]*av[i]%mod;
	ij[m]=qpow(jc[m],mod-2);
	for(int i=m;i;i--)ij[i-1]=1ll*av[i]*ij[i]%mod;
	for(int i=1;i<=m;i++){
		//后面还会接入n(m-i)个点，当前共(n-1)i+1个叶子
		//(n(m-i)/((n-1)i+1)+1)*n
		int xp=1ll*n*(m-i)%mod*jc[i-1]%mod*ij[i]%mod+1;
		(ans+=xp)%=mod;
	}
	return printf("%lld",(1ll*ans*n%mod+mod)%mod),0;
}