2020牛客多校第九场 E-Groundhog Chasing Death

Groundhog Chasing Death

题意

给定a、b、c、d、x、y，算
比如a=c=1,b=d=2,x=4,y=8
那么就算4个gcd的乘积
即gcd(4,8)* gcd(4,64)* gcd(16,8)* gcd(16,64)=4 *4 *8 *16=2048

思路

因为gcd的本质就是把二者的质因子提出来之后，取各质因子的最小次数。
如gcd(2²* 3²,2³* 3)=2^min(2,3) * 3^min(2,1)=2² * 3=12
那么也就是说只要i固定，那么就可以计算j的一个临界点，小于这个临界点就计算i的质因子次数贡献，否则就计算j的质因子次数贡献，把这些质因子次数加起来答案就是这个质因子的次数。
因为质因子次数很大，所以要用到费马小定理
只要a不是p的倍数，则

也就是说
a^p(mod p)=a¹(mod p)
a^p+1(mod p)=a²(mod p)
所以质因子次数只需要模p-1即可（只在p为质数条件下）

代码

#include<bits/stdc++.h>
using namespace std;
const int P=998244353;

long long Pow(long long a,long long p)
{
    static long long ret;
    for(ret=1;p;p>>=1,a=a*a%P)if(p&1)ret=ret*a%P;
    return ret;
}

int a,b,c,d,x,y;

map<int,int>A[2];
void Make_Map(int x,int id)
{
    A[id].clear();
    for(int i=2;i*i<=x;i++)
    {
        if(x%i)continue;
        int c=0;
        while(x%i==0)c++,x/=i;
        A[id][i]=c;
    }
    if(x>1)A[id][x]=1;
}

void Solve()
{
    Make_Map(x,0);
    Make_Map(y,1);
    long long ret=1;
    for(map<int,int>::iterator it=A[0].begin();it!=A[0].end();++it)
    {
        long long f=(*it).first,px=(*it).second,py=A[1][f];
        if(py==0)continue;
        long long ptot=0;
        for(int i=a;i<=b;i++)
        {
            long long ano=max(c-1ll,i*px/py);
            (ptot+=max(0ll,d-ano)*i*px)%=P-1;
        }
        for(int i=c;i<=d;i++)
        {
            long long ano=max(a-1ll,(i*py-1)/px);
            (ptot+=max(0ll,b-ano)*i*py)%=P-1;
        }
        (ret*=Pow(f,ptot))%=P;
    }
    printf("%lld\n",ret);
}

int main()
{
    scanf("%d%d%d%d%d%d",&a,&b,&c,&d,&x,&y);
    Solve();
    return 0;
}