HDU 6050 Funny Function（快速幂）

Description
函数F(x,y)满足：

给出n,m，求F(m,1)%1e9+7
Input
第一行一整数T表示用例组数，每组用例输入两个整数n,m (1<=T<=1e4,1<=n,m < 2^63)
Output
对于每组用例，输出F(m,1)%1e9+7
Sample Input
2
2 2
3 3
Sample Output
2
33
Solution

Code

#include<cstdio>
using namespace std;
typedef long long ll;
const int mod=1000000007,inv3=333333336;
int T;
ll n,m;
int mod_pow(int a,int b)
{
    int ans=1;
    while(b)
    {
        if(b&1)ans=(ll)ans*a%mod;
        a=(ll)a*a%mod;
        b>>=1;
    }
    return ans;
}
int main()
{
    scanf("%d",&T);
    while(T--)
    {
        scanf("%I64d%I64d",&n,&m);
        if(m==1)
        {
            printf("1\n");
            continue;
        }
        int ans=0;
        int a=mod_pow(2,n%(mod-1))-1;
        ans=(ans+2ll*mod_pow(a,(m-1)%(mod-1)))%mod;
        if(n&1)ans++;
        if(ans>=mod)ans-=mod;
        ans=(ll)ans*inv3%mod;
        printf("%d\n",ans);
    }
    return 0;
}