#include <bits/stdc++.h>
using namespace std;
#define int long long
#define maxn 100005
const int p=1e9+7;
int dp[80][1<<8][20];
int N[80],M[80];
vector<int>h[1<<8];
int fac[maxn],inv[maxn];

int qpow(int a,int b){
    int ans=1;
    int tmp=a;
    while(b){
        if(b&1){
            ans=ans*tmp%p;
        }
        tmp=tmp*tmp%p;
        b>>=1;
    }
    return ans;
}

int C(int a,int b){
    if (a < b) return 0;
    return fac[a]*inv[b]%p*inv[a-b]%p;
}

void solve(){
    
    int n,m,k;cin>>n>>m>>k;
    for(int i=60;i>=0;i--){
        N[i]=(n>>i&1);
        M[i]=(m>>i&1);
    }
    dp[60][0][k]=1;
    int ans=0;
    for(int i=60;i>=0;i--){
        for(int x=0;x<=k;x++){
            int y=k-x;
            int v=x*y;
            for(int cnt=0;cnt<=k;cnt++){
                for(int u=0;u<(1<<8);u++){
                    if(dp[i][u][cnt]==0||u<v)continue;
                    int nxt=(u-v)*2+(i>0?N[i-1]:0);
                    if(nxt>=(1<<8))continue;
                    for(int cnt2=0;cnt2<=cnt;cnt2++){
                        if(M[i]==0&&cnt2!=cnt)continue;
                        if(cnt-cnt2<=y){
                            int tp=C(k-cnt,y-cnt+cnt2)*C(cnt,cnt2)%p;
                            if(i>0){
                                dp[i-1][nxt][cnt2]+=dp[i][u][cnt]*tp%p;
                                if(dp[i-1][nxt][cnt2]>=p)dp[i-1][nxt][cnt2]-=p;
                            }
                            else if(nxt==0){
                                ans+=dp[i][u][cnt]*tp%p;
                            }
                        }
                    }
                }
            }
        }
    }

    cout<<ans%p;
}

signed main(){
    ios::sync_with_stdio(0);
    cin.tie(0);
    fac[0]=1;
    for(int i=1;i<maxn;i++){
        fac[i]=fac[i-1]*i%p;
    }
    inv[maxn-1]=qpow(fac[maxn-1],p-2);
    for(int i=maxn-2;i>=0;i--){
        inv[i]=inv[i+1]*(i+1)%p;
    }

    solve();
}