原题：ZOJ 3774 http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3774
----------------------------------------------------------------------------------------------------------------------
这题比较复杂，看这篇比较详细：http://blog.****.net/acdreamers/article/details/23039571
结论就是计算：
充分利用了快速幂及求逆元。
代码：
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>
#include <algorithm>
#define Mod 1000000009
#define ll long long
using namespace std;
#define N 100007 ll fac[N],A[N],B[N]; void init()
{
int i;
fac[] = ;
for(i=;i<N;i++)
fac[i] = fac[i-]*i%Mod;
A[] = B[] = ;
for(i=;i<N;i++)
{
A[i] = A[i-]* % Mod;
B[i] = B[i-]* % Mod;
}
} ll fastm(ll n,ll k,ll MOD)
{
ll res = 1LL;
n %= MOD;
while(k)
{
if(k&1LL)
res = (res*n)%MOD;
k >>= ;
n = n*n%MOD;
}
return res;
} ll Inv(ll n,ll MOD)
{
return fastm(n,MOD-,MOD);
} int main()
{
int cs;
ll n,k;
init();
ll ans,r;
scanf("%d",&cs);
while(cs--)
{
scanf("%lld%lld",&n,&k);
ans = ;
for(r=;r<=k;r++)
{
ll t = A[k-r]*B[r] % Mod;
ll x = fac[k]; // k!
ll y = fac[k-r]*fac[r] % Mod; // (k-r)!*(r)!
ll c = x*Inv(y,Mod) % Mod; // c = C(k,r) = x/y = x*Inv(y)
ll tmp = t*(fastm(t,n,Mod)-)%Mod*Inv(t-,Mod)%Mod; //t(t^n-1)/(t-1) = t(t^n-1)*Inv(t-1)
if(t == )
tmp = n%Mod;
tmp = tmp*c%Mod;
if(r&1LL) // (-1)^r
ans -= tmp;
else
ans += tmp;
ans %= Mod;
}
//ans = ans*(1/sqrt(5))^k
ll m = Inv(,Mod)%Mod;
ans = ans*fastm(m,k,Mod)%Mod;
ans = (ans%Mod+Mod)%Mod;
printf("%lld\n",ans);
}
return ;
}