[2017CCPC杭州] HDU 6265 数论

以${\prod }_{i=1}^k{p}_i^{q_i}$的形式给出一个数$n$，$p_i\le 1e8,q_i\le 1e8$，并且不同的$p_i$不会超过$20$个。求${\sum }_{d|n}\varphi (d)\frac{n}{d}$。
所求的式子是一个与$n$有关的狄利克雷卷积，而$\varphi$和$id$都是积性函数，所求的是式子也是积性函数。对于积性函数，可以单独考虑每个质因数的答案，最后乘在一起。
也就是说对于第$i$个质因数$p_i^{q_i}$，记为$x_i$
$\sum _{d|x_i}\varphi (d)\frac{x_i}{d}$
$=\varphi (1)\ast p_i^{q_i}+\sum _{j=1}^{q_i}\varphi (p_i^j)\frac{p_i^{q_i}}{p_i^j}$
$=p_i^{q_i}+\sum _{j=1}^{q_i}(p_i^j-p_i^{j-1})\frac{p_i^{q_i}}{p_i^j}$
$=p_i^{q_i-1}(p_i-q_i+p_i{q}_i)$
答案就是$\prod _{i=1}^k{p}_i^{q_i-1}(p_i-q_i+p_i{q}_i)$

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int mod=998244353;
const int N=21;
ll p[N],q[N];
void mul(ll &x,ll y) {
	if(x*y>=mod) x=(x*y)%mod;
	else x*=y; 
}
ll fpow(ll x,ll y) {
	ll ans=1;
	while(y) {
		if(y&1) ans=(ans*x)%mod;
		x=(x*x)%mod;
		y>>=1;
	}
	return ans;
}
int main() {
	int T;
	scanf("%d",&T);
	while(T--) {
		int n;
		scanf("%d",&n);
		for(int i=1;i<=n;i++) 
			scanf("%lld%lld",&p[i],&q[i]);
		ll ans=1;
		for(int i=1;i<=n;i++) {
			mul(ans,fpow(p[i],q[i]-1));
			mul(ans,((p[i]-q[i]+mod)%mod+p[i]*q[i]%mod)%mod);
		}
		printf("%lld\n",ans); 
	}
	return 0;
}