WA了7个点，求大佬改错。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#define re register
#define LL long long
#define N 1000005
using namespace std;
const LL Mod=1e9+7; const int lim=1e6;
LL T,n,m,f[N],g[N],mul[N],inv_mul[N],p_cnt,mu[N],p[100005]; bool vis[N];
LL power(LL x,LL P){
	LL ans=1,c=x;
	while(P){
		if(P&1) ans=(ans*c%Mod+Mod)%Mod;
		P>>=1; c=(c*c%Mod+Mod)%Mod;
	}
	return ans;
}
void get_f(){
	f[1]=f[2]=1;
	for(re int i=3;i<=lim;++i) f[i]=(f[i-1]+f[i-2])%Mod;
	for(re int i=1;i<=lim;++i) g[i]=power(f[i],Mod-2);
}
void get_mu(){
	mu[1]=1;
	for(re int i=2;i<=lim;++i){
		if(!vis[i]) p[++p_cnt]=i,mu[i]=-1;
		for(re int j=1;j<=p_cnt;++j){
			if(i*p[j]>lim) break; vis[i*p[j]]=1;
			if(i%p[j]==0) break; mu[i*p[j]]=-mu[i];
		}
	}
}
void get_mul(){
	for(re int i=0;i<=lim;++i) mul[i]=1;
	for(re int i=1;i<=lim;++i)
	for(re int j=1;i*j<=lim;++j){
		if(mu[j]==1) mul[i*j]=(mul[i*j]*f[i]%Mod+Mod)%Mod;
		else if(mu[j]==-1) mul[i*j]=(mul[i*j]*g[i]%Mod+Mod)%Mod;
	}
	for(re int i=1;i<=lim;++i) mul[i]=(mul[i]*mul[i-1]%Mod+Mod)%Mod;
	for(re int i=0;i<=lim;++i) inv_mul[i]=power(mul[i],Mod-2);
}
int main(){
	scanf("%lld",&T); get_f(); get_mu(); get_mul();
	for(re int i=1;i<=T;++i){
		scanf("%lld%lld",&n,&m); LL ans=1;
		if(n>m){ LL t=n; n=m; m=t; }
		for(re LL l=1,r;l<=n;l=r+1){
			r=min(n/(n/l),m/(m/l));
			ans=(ans*power((mul[r]*inv_mul[l-1]%Mod+Mod)%Mod,((n/l)*(m/l)%Mod+Mod)%Mod)%Mod+Mod)%Mod;
		}
		printf("%lld\n",ans);
	}
	return 0;
}