题意:定义 f(x) = x 的因子个数 ^ 3,g(x) = \sum_{i|x} f(x) 。求 g(A^B) mod 10007 。 A,B<=1e6 。
solution: 积性函数好题。
考点:唯一质因数分解 + 积性函数。
首先 f(x) 是积性函数。假设 f(x) = f(u*v) ,(u,v)=1 。设 d(x) = x 的质因子个数,由定义可知 d(x) = d(u) + d(v) ,又因为 f(x) = d(x)^3 = (d(u) + d(v))^3 = d(u)^3 d(v)^3 = f(u) * f(v) ,证毕。
容易发现 g(x) 也是积性函数。 假设 g(x) = g(u*v),(u,v)=1 。那么对于 i|x , i 的一部分 x1 属于 u, 另一部分 x2 则属于 v ,又因为 f(i) = f(x1)*f(x2) ,所以 g(x) = g(u) * g(v) ,证毕。
再利用唯一质因数分解定理,对 A 进行因式分解,只需要解决 g ( p i a i ) g(p_i^{a_i}) g(piai) 。又因为 g ( p i a i ) = ∑ i = 1 a i + 1 i 3 = ( a i + 1 ) 2 ( a i + 2 ) 2 / 4 g(p_i^{a_i}) = ∑_{i=1}^{a_i+1}i^3 = (a_i+1)^2(a_i+2)^2 / 4 g(piai)=∑i=1ai+1i3=(ai+1)2(ai+2)2/4 ,可以直接计算,答案就是 ∏ i g ( p i a i ∗ B ) \prod_{i}{g(p_i^{a_i*B})} ∏ig(piai∗B) 。
#include <bits/stdc++.h>
#define ll long long
#define INF 0x3f3f3f3f
#define PII pair<int,int>
using namespace std;
const int mx=1e6+5;
const int mod=10007;
int n,m,prime[mx],cnt,vis[mx];
signed main() {
// freopen("data.in","r",stdin);
for(int i=2;i<=1e6;i++) {
if(!vis[i]) prime[++cnt]=i,vis[i]=i;
for(int j=1;j<=cnt&&1ll*i*prime[j]<=1e6;j++) {
vis[i*prime[j]]=prime[j]; //用来分离 i 的最小质因数
if(i%prime[j]==0) break;
}
}
int Case(0);
while(~scanf("%d%d",&n,&m)) {
ll ans(1),res;
while(n>1) {
ll tmp=vis[n],cnt=0;
while(n%tmp==0) n/=tmp,cnt++;
cnt=cnt*m%mod;
res=(cnt+1)*(cnt+2)/2%mod;
res=res*res%mod;
ans=ans*res%mod;
}
printf("Case %d: %lld\n",++Case,ans);
}
}