样例过了，交上去全$WA$

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
#define int long long
#define mod 1000000007

int t,n,m,f[1000002],g[1000002],p[1000002],mu[1000002],res[1000002],sum[1000002];
bool vis[1000002];

void getmu() {
	vis[1]=mu[1]=1;
	for(int i=2;i<=1000000;i++) {
		if(!vis[i]) p[++p[0]]=i,mu[i]=-1;
		for(int j=1;j<=p[0] && i*p[j]<=1000000;j++) {
			vis[i*p[j]]=1;
			if(i%p[j]==0) break;
			mu[i*p[j]]=-mu[i];
		}
	}
}

int power(int x,int y) {
	if(y==0) return 1;
	if(y==1) return x;
	if(y%2==0) return power((x*x)%mod,y/2);
	return (x*power(x,y-1))%mod;
}

void getf() {
	g[1]=f[1]=res[1]=1;
	for(int i=2;i<=1000000;i++) f[i]=(f[i-1]+f[i-2])%mod,g[i]=power(f[i],mod-2),res[i]=1;
	for(int i=1;i<=1000000;i++) {
		if(!mu[i]) continue;
		for(int j=1;i*j<=1000000;j++) res[i*j]=(res[i*j]*(mu[i]==1?f[j]:g[j]))%mod;
	}
	sum[0]=1;
	for(int i=1;i<=1000000;i++) (sum[i]=1ll*sum[i-1]*res[i])%mod;
}

signed main() {
	getmu();
	getf();
	scanf("%lld",&t);
	while(t--) {
		scanf("%lld%lld",&n,&m);
		if(n>m) swap(n,m);
		int ans=1,l=1,r=1;
		while(l<=n) {
			r=min(n/(n/l),m/(m/l));
			int inv=(sum[r]*power(sum[l-1],mod-2))%mod;
			ans=(ans*power(inv,(n/l)*(m/l)%(mod-1)))%mod;
			l=r+1;
		}
		printf("%lld\n",(ans+mod)%mod);
	}
	return 0;
}