#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 ;

}