有多组测试数据。

鸣谢infinityedge上传

 1 #include<iostream> 2 #include<algorithm> 3 #include<cstring> 4 #include<cstdio> 5 #include<cmath> 6 #define LL long long 7 using namespace std; 8 const int mxn=1000005; 9 const int mod=1e9+7;10 int read(){11     int x=0,f=1;char ch=getchar();12     while(ch<‘0‘ || ch>‘9‘){if(ch==‘-‘)f=-1;ch=getchar();}13     while(ch>=‘0‘ && ch<=‘9‘){x=x*10+ch-‘0‘;ch=getchar();}14     return x*f;15 }16 int pri[mxn],mu[mxn],f[mxn],cnt=0;17 LL g[mxn];18 bool vis[mxn];19 LL ksm(LL a,LL k){20     LL res=1;21     while(k){22         if(k&1)res=res*a%mod;23         a=a*a%mod;24         k>>=1;25     }26     return res;27 }28 void init(){29     f[1]=1;g[0]=1;mu[1]=1;30     for(int i=2;i<mxn;i++)f[i]=((LL)f[i-1]+f[i-2])%mod;31     for(int i=2;i<mxn;i++){32         if(!vis[i]){pri[++cnt]=i;mu[i]=-1;}33         for(int j=1;j<=cnt && (LL)pri[j]*i<mxn;j++){34             int tmp=pri[j]*i;35             vis[tmp]=1;36             if(i%pri[j]==0){mu[tmp]=0;break;}37             mu[tmp]=-mu[i];38         }39     }40     for(int i=1;i<mxn;i++)g[i]=1;41     for(int i=1;i<mxn;i++){42         LL invi=ksm(f[i],mod-2);43         for(int j=i,c=1;j<mxn;j+=i,c++){44             if(!mu[c])continue;45             if(mu[c]==1)g[j]=(LL)g[j]*f[i]%mod;46             else g[j]=(LL)g[j]*invi%mod;47         }48     }49     for(int i=2;i<mxn;i++)g[i]=(LL)g[i]*g[i-1]%mod;50     return;51 }52 void calc(int n,int m){53     if(n>m)swap(n,m);54     LL ans=1;55     for(int i=1,pos;i<=n;i=pos+1){56         LL x=n/i,y=m/i;57         pos=min(n/x,m/y);58         LL inv=ksm(g[i-1],mod-2);59         ans=ans*ksm(g[pos]*inv%mod,(LL)x*y%(mod-1))%mod;60     }61     printf("%lld\n",ans);62     return;63 }64 int n,m;65 int main(){66 //  freopen("product.in","r",stdin);67 //  freopen("product.out","w",stdout);68     int i,j;69     int T=read();70     init();71     while(T--){72         n=read();m=read();73         calc(n,m);74     }75     return 0;76 }