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

https://ac.nowcoder.com/acm/contest/5674/E

考虑x和y的每个共有的质因子在答案中的幂次q，ans*=p[k]^q

r[k]为p[k]在x中的幂次，s[k]为p[k]在y中的幂次

枚举i=a->b，然后i*r[k]

由于q可能很大，注意幂次是要%phi(mod)也就是%(mod-1)

注意特判当x=mod,y=mod的情况

#include
using namespace std;
const int mod=998244353,phi=mod-1;
int a,b,c,d;
int x,y;
int p[105],r[105],s[105],tot;
long long S(long long n){
	return n*(n+1)/2;
}
long long qpow(long long x,int n){
	long long ans=1;
	for(;n;n>>=1,x=x*x%mod)if(n&1)ans=ans*x%mod;
	return ans;
}
long long f(int k){
	if(!s[k])return 1;
	long long q=0;
	for(int i=a;i<=b;i++){
		int v=i*r[k]/s[k];
		if(vy)swap(x,y);
	for(int i=2;i*i<=x;i++){
		if(x%i==0){
			tot++;
			p[tot]=i;
			while(x%i==0)r[tot]++,x/=i;
			while(y%i==0)s[tot]++,y/=i;
		}
	}
	if(x!=1){
		tot++;
		p[tot]=x;
		r[tot]=1;
		while(y%x==0)s[tot]++,y/=x;
	}
	long long ans=1;
	for(int i=1;i<=tot;i++)ans=ans*f(i)%mod;
	printf("%lld",ans);
	return 0;
}