BZOJ 3601 一个人的数论

第一次被搞。

这个套路深啊。。。。。。
$$\begin{gathered} \sum _{i=1}^{n-1}[\gcd (i,n)==1]i^d \\ =\sum _{i=1}^n\sum _{p|i,p|n}\mu (p)i^d \\ =\sum _{p|n}^n\mu (p)\sum _{i=1}^{\frac{n}{p}}(ip{)}^d \\ =\sum _{p|n}^n\mu (p)p^d\sum _{i=1}^{\frac{n}{p}}{i}^d \end{gathered}$$
然后对于自然数幂和我们可以用拉格朗日插值将它改写为关于上界的$d+1$次多项式：设${\sum }_{i=1}^x{i}^d={\sum }_{i=0}^{d+1}{a}_i{x}^i$
那么原式$$\begin{gathered} =\sum _{p|n}^n\mu (p)p^d\sum _{i=0}^{d+1}{a}_i(\frac{n}{p}{)}^i \\ =\sum _{i=0}^{d+1}{a}_i\sum _{p|n}^n\mu (p)p^d(\frac{n}{p}{)}^i \end{gathered}$$
那么如果可以算出后面一个和式的值久可以算出答案了。
后面一个和式发现是三个积性函数的点积，是可以在令$n=p_i^{a_i}$的每个$i=1....w$处算出答案然后乘起来得到的。

$\mathrm{A}\mathrm{C}\mathrm{C}\mathrm{o}\mathrm{d}\mathrm{e}$

#include<bits/stdc++.h>
#define maxn 1005
#define mod 1000000007
using namespace std;

int d,w,s[maxn],inv[maxn]={1,1},invf[maxn]={1,1},rinv[maxn]={1,-1},dp[maxn]={1},x[maxn],p[maxn],a[maxn];
int Pow(int base,int k){
	int ret = 1;
	for(;k;k>>=1,base=1ll*base*base%mod)
		if(k&1)
			ret=1ll*ret*base%mod;
	return ret;
}

int main(){
	scanf("%d%d",&d,&w);
	for(int i=1;i<=d+2;i++){
		s[i] = (s[i-1] + Pow(i,d)) % mod;
		for(int j=d+2;j>=0;j--)
			dp[j] = ((j ? dp[j-1] : 0) - 1ll * i * dp[j]) % mod;
	}
	for(int i=2;i<=d+2;i++) 
		inv[i] = 1ll * (mod - mod / i ) * inv[mod % i] % mod,
		invf[i] = 1ll * invf[i-1] * inv[i] % mod,
		rinv[i] = 1ll * rinv[i-1] * (-inv[i]) % mod;
	for(int i=1;i<=d+2;i++){
		for(int j=0;j<=d+2;j++){
			dp[j] = ((j ? dp[j-1] : 0) - dp[j]) * 1ll * inv[i] % mod,
			x[j] = (x[j] + dp[j] * 1ll * invf[i-1] % mod * rinv[d+2-i] % mod * s[i]) % mod;
		}
		for(int j=d+2;j>=0;j--)
			dp[j] = ((j ? dp[j-1] : 0) - 1ll * i * dp[j]) % mod;
	}
	for(int i=1;i<=w;i++) scanf("%d%d",&p[i],&a[i]);
	int ans = 0;
	for(int i=0;i<=d+1;i++){
		int prd=1;
		for(int j=1;j<=w;j++)
			prd = 1ll * prd * (Pow(Pow(p[j],a[j]),i) - 1ll * Pow(p[j],d) * Pow(Pow(p[j],a[j]-1),i) % mod) % mod;
		ans = (ans + 1ll * prd * x[i]) % mod;
	}
	printf("%d\n",(ans+mod)%mod);
}