2020牛客暑期多校训练营(第九场)—— Groundhog Chasing Death

2020牛客暑期多校训练营(第九场)—— Groundhog Chasing Death

输入

样例1

1 2 1 2 8 4

样例2

1 2 3 4 120 180

输出

样例1

2048

样例2

235140177

题目大意

\displaystyle\prod_{i=a}^b\prod_{j=c}^d\gcd(x^i,y^j) mod 998244353

题解

对给出的数分解质因数,然后对于每个质因数分别讨论其幂次。那么问题就从求 乘积问题转化为o(log x + log y)个子问题。对于每个子问题,在幂次上,形如给出两个数x' , y',求\sum_{i=a}^b\sum_{j=c}^dmin(x'^i,y'^j)。 那么可以枚举这个最小值是多少,是x'^i取到还是在y'^j取到,并可以𝑂(1)求出每种情况的方案数。注意处理好x'^i=y'^j的情况,注意如果要对幂次取模,按照欧拉/费马小定理,模数应该取 \phi(998244353)。令n=max\left ( b,d \right ),这种算法的最终复杂度是O(n(log x + log y))的,通过本题已经绰绰有余。

AC Code

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int MAXN=1e5+10;
const int MOD=998244353;
const int mod=MOD-1;
ll ksm(ll a,ll b)
{
	ll ret=1;
	while(b)
	{
		if(b&1)ret=ret*a%MOD;
		a=a*a%MOD,b>>=1;
	}
	return ret;
}
ll a,b,c,d,x,y,ans=1;
void solve(int p)
{
	ll aa=0,bb=0;
	while(x%p==0)aa++,x/=p;
	while(y%p==0)bb++,y/=p;
	if(aa==0||bb==0)return;
	ll tp=0;
	for(ll i=a;i<=b;i++)
	{
		ll j=aa*i/bb,tmp,now;
		if(j<c)tmp=(d-c+1)*i%mod*aa%mod;
		else if(j>=d)tmp=(c+d)*(d-c+1)/2%mod*bb%mod;
		else now=(c+j)*(j-c+1)/2%mod*bb%mod,tmp=(now+(d-j)*i%mod*aa%mod)%mod;
		tp=(tp+tmp)%mod;
	}
	ans=1ll*ans*ksm(p,tp)%MOD;
}
int main()
{
	scanf("%lld%lld%lld%lld%lld%lld",&a,&b,&c,&d,&x,&y);
	for(int k=2;k<=x;k++)if(x%k==0)solve(k);
	if(x>1)solve(x);
	printf("%lld\n",ans%MOD);
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值