XDOJ-1201-Dogs of Qwordance Senior Backend R&D Engineers（数论）

题目链接：XDOJ-1201-Dogs of Qwordance Senior Backend R&D Engineers

设 $f(x)$ 为 $x$ 的因子数量。显然：

$$ans=\sum_{i|x}f(i)f(\frac{x}{i})$$
但是这样似乎不可搞。。
将 $x$ 分解为质数乘积的形式
$$x=\prod_{i=1}^k p_i^{a_i}$$
那么 $x$ 的任意一个因子 $x'=\prod_{i=1}^k p_i^{b_i}$，$b_i\leq a_i$ 。因此 $a_i$ 划分成 $b_i$ 和 $a_i-b_i$ 后就得到 $x$ 的两个因子。
当我们划分 $a_i$ 时，共有 $a_i+1$ 种方案。
例如 $2^3$ 可以划分成 $1*8$ ，$2*4$，$4*2$，$8*1$ 。
则方案数为 $\sum_{j=0}^{a_i}f(p_i^j)f(p_i^{a_i-j})=\sum_{j=0}^{a_i}(j+1)(a_i-j+1)$
对于每个 $a_i$ 求方案数，根据乘法原理：

$$ans=\prod_{i=1}^k\sum_{j=0}^{a_i}(j+1)(a_i-j+1)$$

然后可以优化一下，首先是用线性筛$O(n)$筛出素数。
然后有了素数后，对于每个 $x$ 找出其素因子个数就行了。
复杂度为 $O(n+\frac{x}{lnx})$

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e6+7;
bool vis[maxn];
ll p[80007];
int main()
{
    ll x;
    int k=0;
    for(int i=2;i<maxn;i++)
    {
        if(!vis[i]) p[k++]=i;
        for(int j=0;j<k&&p[j]*i<maxn;j++)
        {
            vis[i*p[j]]=true;
            if(i%p[j]==0) break;
        }
    }
    while(~scanf("%lld",&x))
    {
        ll ans=1;
        for(ll i=0;p[i]*p[i]<=x;i++)
        {
            ll cnt=0;
            while(x%p[i]==0) cnt++,x/=p[i];
            ll res=0;
            for(ll i=0;i<=cnt;i++)
                res+=(i+1)*(cnt-i+1);
            ans*=res;
        }
        if(x!=1) ans*=4;
        printf("%lld\n",ans);
    }
    return 0;
}