[质因数分解+费马小定理] 牛客多校2020 第九场 E-Groundhog Chasing Death

题目

题目链接：https://ac.nowcoder.com/acm/contest/5674/E

思路

简单说下思路
算了不会写式子还是不说了吧
这里有大佬博客！->大佬博客
讲的非常清楚qwq

代码

#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<cctype>
#include<ctime>
#include<iostream>
#include<string>
#include<map>
#include<queue>
#include<stack>
#include<set>
#include<vector>
#include<iomanip>
#include<list>
#include<bitset>
#include<sstream>
#include<fstream>
#include<complex>
#include<algorithm>
#if __cplusplus >= 201103L
#include <unordered_map>
#include <unordered_set>
#endif
#define int long long
using namespace std;
const int INF = 0x3f3f3f3f;
const int mod = 998244353;
int a,b,c,d,x,y;
int ans=1;
int quick_pow(int a,int b){
	int res=1;
	while(b){
		if(b&1){
			res=res*a%mod;
		}
		a=a*a%mod;
		b>>=1;
	}
	return res;
}
void solve(int p){
	int cnt1=0,cnt2=0,res=0;
	while(x%p==0) x/=p,cnt1++;
	while(y%p==0) y/=p,cnt2++;
	if(cnt1==0||cnt2==0) return;
	for(int i=a;i<=b;i++){
		int j=i*cnt1/cnt2,tmp;
		if(j<c){
			tmp=(d-c+1)*i%(mod-1)*cnt1%(mod-1);
		}
		else if(j>=d){
			tmp=(d+c)*(d-c+1)/2%(mod-1)*cnt2%(mod-1);
		}
		else{
			tmp=(j+c)*(j-c+1)/2%(mod-1)*cnt2%(mod-1);
			tmp=((d-j)*i%(mod-1)*cnt1%(mod-1)+tmp)%(mod-1);
		}
		res=(res+tmp)%(mod-1);
	}
	ans=ans*quick_pow(p,res)%mod;
}
signed main(){
	ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
	cin>>a>>b>>c>>d>>x>>y;
	for(int i=2;i<=x;i++){
		if(x%i==0) solve(i);//遇到一个因数进行一次操作 
	}
	cout<<ans<<endl;
    return 0;
}