hdu 6390 GuGuFishtion(欧拉函数+容斥)

稍加推导,就可以转化为一个非常经典非常简单的式子,求1-n,1-m中gcd(a,b)=k,k的出现次数.
容斥即可,注意反向.
但是我对这个容斥不太清晰,所以先不写.
非常对不起,这个题目我做不出来真的是非常非常不应该
附上巨神的代码,怕自己到时候不会了..

#include <bits/stdc++.h>
using namespace std;
#define rep(i,a,n) for (int i=a;i<n;i++)
#define per(i,a,n) for (int i=n-1;i>=a;i--)
#define pb push_back
#define mp make_pair
#define all(x) (x).begin(),(x).end()
#define fi first
#define se second
#define SZ(x) ((int)(x).size())
typedef vector<int> VI;
typedef long long ll;
typedef pair<int,int> PII;
ll gcd(ll a,ll b) { return b?gcd(b,a%b):a;}
// head

const ll N=1000000;
int pr[N/5],p[N+100],tot,phi[N+100];
ll f[N+10],inv[N+10];
int n,m,mod,_;

void init() {
    phi[1]=1;p[1]=1;
    rep(i,2,N+1) {
        if (!p[i]) p[i]=i,pr[++tot]=i,phi[i]=p[i]-1;
        for (int j=1;j<=tot&&pr[j]*i<=N;j++) {
            p[i*pr[j]]=pr[j];
            if (p[i]==pr[j]) {
                phi[i*pr[j]]=phi[i]*pr[j];
                break;
            } else phi[i*pr[j]]=phi[i]*(pr[j]-1);
        }
    }
}

int main() {
    init();
    for (scanf("%d",&_);_;_--) {
        scanf("%d%d%d",&n,&m,&mod);
        if (n>m) swap(n,m);
        inv[1]=1;
        for (int i=2;i<=n;i++) inv[i]=(mod-mod/i)*inv[mod%i]%mod;
        ll ans=0;
        for (int i=n;i>=1;i--) {
            f[i]=(ll)(n/i)*(m/i);
            for (int j=i+i;j<=n;j+=i) f[i]=f[i]-f[j];
            ans=(ans+f[i]%mod*i%mod*inv[phi[i]])%mod;
        }
        printf("%lld\n",ans);
    }
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值