累乘函数线性逆元打表，阶乘反演——bzoj4816

学了一种新套路，倒序打表函数的逆元可以直接线性完成

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define LL long long
const long long mod = 1e9+;
#define maxn 1000001

LL Pow(LL a,LL b){
    if (!b)
        return 1LL;
    LL x=Pow(a,b/);
    x=x*x%mod;
    if (b&1LL)
        x=x*a%mod;
    return x;
}
/*ll Pow(ll a,ll b){
    ll res=1;
    while(b){
        if(b%2)
            res=res*a%mod;
        b>>=1;a=a*a%mod;
    }
    return res;
}*/
ll n,m;
ll F[maxn],pre[maxn],invF[maxn];
void init1(){
    F[]=;F[]=;
    for(int i=;i<maxn;i++)
        F[i]=(F[i-]+F[i-])%mod;
    pre[]=;
    for(int i=;i<maxn;i++)//临时数组算累乘
        pre[i]=pre[i-]*F[i]%mod;

    ll tmp=Pow(pre[maxn-],mod-);
    for(int i=maxn-;i>=;i--)
        invF[i]=tmp*pre[i-]%mod,tmp=tmp*F[i]%mod;
}

bool vis[maxn];
ll prime[maxn],mm,mu[maxn];
void init2(){
    mu[]=;
    for(int i=;i<maxn;i++){
        if(!vis[i]){
            mu[i]=-;
            prime[++mm]=i;
        }
        for(int j=;j<=mm;j++){
            if(i*prime[j]>=maxn)break;
            vis[i*prime[j]]=;
            if(i%prime[j]==){
                mu[i*prime[j]]=;
                break;
            }
            else mu[i*prime[j]]=-mu[i];
        }
    }
}

ll mul[maxn],invmul[maxn],g[maxn];
void init3(){
    for(int i=;i<maxn;i++)g[i]=;
    for(int i=;i<maxn;i++)
        for(int j=;j*i<maxn;j++){
            if(mu[j]==)
                g[i*j]=g[i*j]*F[i]%mod;
            else if(mu[j]==-)
                g[i*j]=g[i*j]*invF[i]%mod;
        }
    mul[]=;
    for(int i=;i<maxn;i++)
        mul[i]=mul[i-]*g[i]%mod;
    invmul[maxn-]=Pow(mul[maxn-],mod-);
    for(int i=maxn-;i>=;i--)
        invmul[i]=invmul[i+]*g[i+]%mod;
}

int main(){
    int t;cin>>t;
    init1();
    init2();
    init3();
    while(t--){
        cin>>n>>m;
        if(n>m)swap(n,m);
        ll ans=;
        for(int l=,r;l<=n;l=r+){
            r=min(n/(n/l),m/(m/l));
            ll tmp=mul[r]*invmul[l-]%mod;
            ans=ans*Pow(tmp,(n/l)*(m/l)%(mod-))%mod;
        }
        cout<<ans<<'\n';
    }
}