// #include<bits/stdc++.h>
// using namespace std;
// int n,b[1005];
// vector<pair<int,int> >ans;
// void solve(int l,int r){
// 	if(l==r) return ;
// 	int mid=(l+r)>>1;
// 	for(int i=1;i<=n;++i) if(b[i]>mid&&b[i]<=r) ans.emplace_back(make_pair(2,i));
// 	for(int i=l+1;i<=mid;++i) ans.emplace_back(make_pair(1,i));
// 	solve(l,mid);solve(mid+1,r);
// }
// int main(){
// 	scanf("%d",&n);
// 	for(int i=1;i<=n;++i) scanf("%d",&b[i]);
// 	solve(0,n);
// 	printf("%d\n",ans.size());
// 	for(auto a:ans) printf("%d %d\n",a.first,a.second);
// 	return 0;
// }
//A

// #include<bits/stdc++.h>
// #define M 100005
// using namespace std;
// int n,_4,_5,fa[M],Fa[M],dep[M],du[M],U,V;
// long long ans;
// vector<int> e[M];
// void dfs(int u,int fa){
// 	for(int v:e[u]) if(v!=fa){
// 		dep[v]=dep[u]+1;Fa[v]=u;
// 		dfs(v,u);
// 	}
// }
// void del(int x){
// 	_4-=(du[x]==4);_5-=(du[x]==5);
// 	du[x]--;
// 	_4+=(du[x]==4);_5+=(du[x]==5);
// }
// void add(int x){
// 	_4-=(du[x]==4);_5-=(du[x]==5);
// 	du[x]++;
// 	_4+=(du[x]==4);_5+=(du[x]==5);
// }
// int find(int x){return x==fa[x]?x:fa[x]=find(fa[x]);}
// int main(){
// 	scanf("%d",&n);
// 	for(int i=1;i<=n;++i) fa[i]=i;
// 	for(int i=1,u,v;i<=n;++i){
// 		scanf("%d%d",&u,&v);du[u]++;du[v]++;
// 		if(find(u)!=find(v)) e[u].emplace_back(v),e[v].emplace_back(u),fa[find(u)]=find(v);
// 		else U=u,V=v;
// 	}
// 	// printf("%d %d\n",U,V);
// 	for(int i=1;i<=n;++i) _4+=(du[i]==4),_5+=(du[i]==5);
// 	dfs(1,0);
// 	del(U);del(V);if(!_5) ans+=n-_4;add(U);add(V);
// 	while(U!=V){
// 		if(dep[V]>dep[U]) swap(U,V);
// 		del(U);del(Fa[U]);
// 		if(!_5) ans+=n-_4;
// 		add(U);add(Fa[U]);U=Fa[U];
// 	}
// 	printf("%lld",ans);
// 	return 0;
// }
//F

#include<bits/stdc++.h>
#define mod 1000000009
using namespace std;
long long m,k;
int mi(int a,int b=mod-2){
	int ans=1;
	while(b){
		if(b&1) ans=1ll*ans*a%mod;
		b>>=1;a=1ll*a*a%mod;
	}
	return ans;
}
int main(){
	scanf("%lld%lld",&m,&k);
	long long nw=0,lf=0;
	for(int i=1;i<=k;++i){
		nw=(1ll*(lf+1)*m%mod+1ll*nw*(1+(i-1)*m)%mod)*mi((1ll*i*m+1)%mod)%mod;
		lf+=1ll*m*mi((1ll+1ll*i*(m-1))%mod)%mod;lf%=mod;
	}
	nw=nw*(1+1ll*k*m)%mod;
	printf("%lld",nw);	
	return 0;
}