HDU -- 2421 Deciphering Password（基础数论）

题目vj链接

题面：

题意：
给定两个数$A,B$，求$A^B$的各因子的因子个数的立方之和。

题解：
我们设 $d$ 为因子个数函数，$d(n)$ 表示 $n$ 的正因子个数。
那么 $d$ 函数为积性函数，若$gcd(n,m)=1$，那么$d(nm)=d(n)\ast d(m)$，那么则有：
$d^3(nm)=d^3(n)\ast d^3(m)$

我们设唯一分解$x=p_1^{e_1}\ast p_2^{e_2}\ast p_3^{e_3}\ast ...\ast p_k^{e_k}$，那么$d(x)=(e_1+1)\ast (e_2+1)\ast (e_3+1)\ast ...\ast (e_k+1)$

现在我们设$A^B=p_1^{e_1}\ast p_2^{e_2}\ast p_3^{e_3}\ast ...\ast p_k^{e_k}$，我们以以下形式表示 $A^B$的因子：

$(p_1^0,p_1^1,p_1^2,...,p_1^{e_1})\ast (p_2^0,p_2^1,p_2^2,...,p_2^{e_2})\ast ...\ast (p_k^0,p_k^1,p_k^2,...,p_k^{e_k})$

我们可以从每一对小括号中选择一个数来组成 $A^B$，且每个小括号选出来的数与其他小括号选出来的数一定都是互质的。

假设我们现在从第一个小括号里面选出来的数是$p_1^0$，其后的数任选均可组成一个 $A^B$的因子，那么现在因子个数的立方的贡献为$((0+1{)}^3)\ast (1^3+2^3+3^3+...+(e_2+1{)}^3)\ast ...\ast (1^3+2^3+3^3+...+(e_k+1{)}^3)$

假设我们现在从第一个小括号里面选出来的数是$p_1^1$，其后的数任选均可组成一个 $A^B$的因子，那么现在因子个数的立方的贡献为$((1+1{)}^3)\ast (1^3+2^3+3^3+...+(e_2+1{)}^3)\ast ...\ast (1^3+2^3+3^3+...+(e_k+1{)}^3)$

那么最终答案$ans=(1^3+2^3+3^3+...+(e_1+1{)}^3)\ast (1^3+2^3+3^3+...+(e_2+1{)}^3)\ast ...\ast (1^3+2^3+3^3+...+(e_k+1{)}^3)$

我们有公式${\sum }_{i=1}^n{i}^3=({\sum }_{i=1}^{n}i{)}^2=(\frac{n\ast (n+1)}{2}{)}^2$

代码：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<string>
#include<queue>
#include<bitset>
#include<map>
#include<unordered_map>
#include<set>
#define ui unsigned int
#define ll long long
#define llu unsigned ll
#define ld long double
#define pr make_pair
#define pb push_back
#define lc (cnt<<1)
#define rc (cnt<<1|1)
#define len(x)  (t[(x)].r-t[(x)].l+1)
#define tmid ((l+r)>>1)
#define fhead(x) for(int i=head[(x)];i;i=nt[i])
#define max(x,y) ((x)>(y)?(x):(y))
#define min(x,y) ((x)>(y)?(y):(x))
using namespace std;

const int inf=0x3f3f3f3f;
const ll lnf=0x3f3f3f3f3f3f3f3f;
const double dnf=1e18;
const int mod=10007;
const double eps=1e-8;
const double pi=acos(-1.0);
const int hp=13331;
const int maxn=5010;
const int maxm=100100;
const int maxp=100100;
const int up=100100;

int cnt=0;
int tot=0;
ll p[108],e[108];
ll a,b;

void only(ll n)
{
    memset(e,0,sizeof(e));
    tot=0;
    for(int i=2;i*i<=n;i++)
    {
        if(n%i) continue;
        p[++tot]=i;
        while(n%i==0)
        {
            e[tot]++;
            n/=i;
        }
    }
    if(n>1) p[++tot]=n,e[tot]=1;
    for(int i=1;i<=tot;i++)
        e[i]*=b;
}

int main(void)
{
    int tt=0;
    while(scanf("%lld%lld",&a,&b)!=EOF)
    {
        only(a);

        ll ans=1;
        ll res=1;
        for(int i=1;i<=tot;i++)
        {
            res=(e[i]+1)*(e[i]+2)/2;
            ans=(ans*res%mod*res)%mod;
        }
        printf("Case %d: %lld\n",++tt,ans%mod);
    }
    return 0;
}