HDU6390多校GuGuFishtion

Pandauncle

于 2018-10-23 20:14:58 发布

阅读量214

点赞数 1

分类专栏：数论

本文链接：https://blog.csdn.net/Pandauncle/article/details/83314002

版权

数论专栏收录该内容

30 篇文章 0 订阅

订阅专栏

题意：就是要求那个公式和。
首先：对一个数中的一个素因子进行分析贡献。
设a=p^a1， b=p^a2;
则phi(ab)=(p-1)*p^(a1+a2-1);
phi(a)=(p-1)*p^(a1-1);
phi(b)=(p-1)*p^(a2-1);
当a1!=0,a2!=0时；
phi(ab)/(phi(a)*phi(b))=p/(p-1)
当a1=0时；
phi(ab)/(phi(a)*phi(b))=1;
当a2=0时;
phi(ab)/(phi(a)phi(b))=1;
所以
phi(ab)/(phi(a) * phi(b))=gcd(a,b)/phi(gcd(a,b));
所以该问题就转换为求1~ n, 1 ~ m区间内有多少对gcd(a,b)=k(1<=a<=n，1<=b<=m)；
设f(d)表示d|gcd(a,b)的对数;
g(d)表示d=gcd(a,b)的对数;
则f(d)=(n/d) * (m/d);
所以进行容斥：
g(d)=f(d)-sum(g(id))(i>=2);

#include<bits/stdc++.h>
using namespace std;
typedef long long int ll;
const ll N=1000000;
int phi[N+10],prime[N+10],tot,ans;
bool mark[N+10];
void getphi()
{
    int i,j;
    phi[1]=1;
    for(i=2;i<=N;i++)
    {
        if(!mark[i])
        {
            prime[++tot]=i;
            phi[i]=i-1;
        }
        for(j=1;j<=tot;j++)
        {
            if(i*prime[j]>N)break;
            mark[i*prime[j]]=1;
            if(i%prime[j]==0)
            {
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
            else
            phi[i*prime[j]]=phi[i]*(prime[j]-1);
        }
    }
}
ll inv[N+10];
ll g[N+10];
int main()
{
    getphi();
    int t;
    scanf("%d",&t);
    while(t--)
    {
        ll n,m,p;
        scanf("%lld%lld%lld",&n,&m,&p);
        memset(inv,0,sizeof(inv));
        memset(g,0,sizeof(g));
        if(n>m)swap(n,m);
        inv[1]=1;
        ll ans=0;
        for(int i=2;i<=n;i++)inv[i]=inv[p%i]*(p-p/i)%p;
        for(ll i=n;i>=1LL;i--)
        {
            g[i]=(ll)(n/i)*(m/i);
            for(ll j=i+i;j<=n;j+=i)
            g[i]=g[i]-g[j];
            ans=(ans+g[i]%p*inv[phi[i]]%p*i%p)%p;
        }
        printf("%lld\n",ans);
    }
    return 0;
}