【SDOI2017】数字表格

题目大意

给定n,m,求 Πni=1Πmj=1fib[gcd(i,j)]
10% 1n,m100
30% 1n,m1000
30% T3
100% T1000,1n,m106

题目解析

一看见题,我就信心满满的打了发莫比乌斯反演正解,然后0分滚粗,还没有暴力容斥60分高。。。后来发现是对某些幂指数取了模的原因。

好的,开始正解:
Πni=1Πmj=1fib[gcd(i,j)]
我们们不妨令 f(x)=Πd|xg(d)
则有 g(x)=f(x)Πd|x,dxg(d)
于是我们可以递推求出 f(x) g(x) 。。。
然后再看原来的方程
Πni=1Πmj=1fib[gcd(i,j)]
=Πni=1Πmj=1Πd|i,d|jg(d)
=Πmin(n,m)d=1g(d)nd×md
又因为 nd×md 仅有 n+m 个取值,搞搞就行了。

代码

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
#include<queue>
#include<stack>
using namespace std;

#define MAXN 1000000
#define MAXM
#define INF 0x3f3f3f3f
typedef long long int LL;

template<class T>
void Read(T &x){
    x=0;char c=getchar();bool flag=0;
    while(c<'0'||'9'<c){if(c=='-')flag=1;c=getchar();}
    while('0'<=c&&c<='9'){x=x*10+c-'0';c=getchar();}
    if(flag)x=-x;
}

const int MOD = 1000000007;
int n,m;
int fib[MAXN+10],g[MAXN+10];
int sum[MAXN+10],inv[MAXN+10];

int ksm(int a,int p,int mod){
    int ret=1;
    while(p){
        if(p&1)ret=1ll*ret*a%mod;
        a=1ll*a*a%mod;
        p>>=1;
    }
    return ret;
}

void init(){
    g[1]=g[2]=fib[1]=fib[2]=1;
    for(int i=3;i<=MAXN;++i)g[i]=fib[i]=(fib[i-2]+fib[i-1])%MOD;
    for(int i=1;i<=MAXN;++i){
        int inv=ksm(g[i],MOD-2,MOD);
        for(int j=i*2;j<=MAXN;j+=i)g[j]=1ll*g[j]*inv%MOD;
    }

    inv[0]=sum[0]=1;
    for(int i=1;i<=MAXN;++i)
        sum[i]=1ll*sum[i-1]*g[i]%MOD,inv[i]=ksm(sum[i],MOD-2,MOD);
}

int main(){
    //freopen("product.in","r",stdin);
    //freopen("product.out","w",stdout);

    init();

    int T;
    Read(T);
    while(T--){
        Read(n),Read(m);
        if(n>m)swap(n,m);

        int ans=1,last=0;
        for(int i=1;i<=n;i=last+1){
            last=min(n/(n/i),m/(m/i));
            ans=1ll*ans*ksm(1ll*sum[last]*inv[i-1]%MOD,1ll*(n/i)*(m/i)%(MOD-1),MOD)%MOD;
        }

        printf("%d\n",ans);
    }
}
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值