2020牛客多校E- Groundhog Chasing Death

原题

题意

求$$\prod _{i=a}^b\prod _{j=c}^d\gcd (x^i,y^j)$$

分析

拆一下质因子。先枚举$i$，然后$i$相同时，用$p$表示当前的质因数，设$x$中有$u$个$p$，$y$中有$v$个$p$，比较一下$i\ast u$和$j\ast v$的大小即可。把小的乘到答案里去。定义变量$P$来记录循环时$p$需要乘的次数。
设$tmp$为分割点，则$tmp=k\ast u/v$。若$tmp<c$，说明$c$到$d$之间$i\ast u$一直小于$j\ast v$，把$d-c+1$个当前的$i\ast u$加入到$P$中，若$tmp>d$，说明$c$到$d$之间$i\ast u$一直大于$j\ast v$，等差数列把所有的j加起来乘$v$加入到$P$中，如果$tmp$在$c$和$d$中间则说明$i\ast u$和$j\ast v$的大小关系并不恒定，老老实实来就行。每次循环结束时答案乘$p^P$即可。
剩下就是瞎搞了。

#include<iostream>
#define ll long long
using namespace std;
const int mod=998244353;
ll qpow(ll a,ll x)
{
  int res=1;
  while(x)
  {
    if(x&1)res=res*a%mod;
    a=a*a%mod;
    x>>=1;
  }
  return res;
}
int p[205],p2[205],q[205],q2[205],ln,ln2;
int main()
{
  int a,b,c,d,x,y;
  cin>>a>>b>>c>>d>>x>>y;
  for(int i=2;i<=x;i++)
  {
    if(x%i==0)p[++ln]=i;
    while(x%i==0)
    {
      q[ln]++;
      x/=i;
    }
  }
  for(int i=2;i<=y;i++)
  {
    if(y%i==0)p2[++ln2]=i;
    while(y%i==0)
    {
      q2[ln2]++;
      y/=i;
    }
  }
  ll ans=1;
  for(int i=1;i<=ln;i++)
  {
    for(int j=1;j<=ln2;j++)
    {
      if(p[i]==p2[j])
      {
        int u=q[i],v=q2[j];
        ll P=0;
        for(int k=a;k<=b;k++)//枚举i
        {
          ll tmp=1ll*k*u/v;
          if(tmp<c)P=(P+1ll*(d-c+1)*k%(mod-1)*u%(mod-1))%(mod-1);
          else if(tmp>d)P=(P+1ll*(d-c+1)*(c+d)/2%(mod-1)*v%(mod-1))%(mod-1);//等差数列求和
          else
          {
            P=(P+1ll*(tmp-c+1)*(c+tmp)/2%(mod-1)*v%(mod-1))%(mod-1);
            P=(P+1ll*(d-tmp)*k%(mod-1)*u%(mod-1))%(mod-1);
          }
        }
        ans=ans*qpow(p[i],P)%mod;
      }
    }
  }
  cout<<ans;
  return 0;
}