bzoj3560 DZY Loves Math V

最新推荐文章于 2018-08-25 08:40:02 发布

sdfzyhx

最新推荐文章于 2018-08-25 08:40:02 发布

阅读量362

点赞数

分类专栏：数学 bzoj

本文链接：https://blog.csdn.net/sdfzyhx/article/details/72967831

版权

bzoj 同时被 2 个专栏收录

251 篇文章 0 订阅

订阅专栏

数学

229 篇文章 0 订阅

订阅专栏

根据欧拉函数的积性，我们可以把每个质因子单独考虑。可以得到

= \sum i 1 = 0 q 1 \dots \sum i n = 0 q n φ (p \sum n j - 1 i j) p - 1 p [(\prod i = 1 n \sum j = 0 q i p j) - 1] + 1

$\begin{align} & \sum_{i_1=0}^{q_1}\cdots\sum_{i_n=0}^{q_n}\varphi(p^{\sum_{j-1}^ni_j}) \\ = & \frac{p-1}p[(\prod_{i=1}^n\sum_{j=0}^{q_i}p^j)-1]+1 \end{align} \\$
因为

qi $q_i$ 很小，直接暴力计算就可以了。

#include<cstdio>
#include<vector>
#include<algorithm>
#include<cmath>
using namespace std;
#define LL long long
const int p=1000000007,maxn=100010,maxx=3200;
int prm[maxx+10],vis[maxx+10],a[maxn],f[maxn],s[maxx],pw[maxx],n,tot,num;
vector<int> v[maxx+10];
int pow(int b,int k)
{
    int ret=1;
    for (;k;k>>=1,b=(LL)b*b%p)
        if (k&1) ret=(LL)ret*b%p;
    return ret;
}
int inc(int x,int y)
{
    x+=y;
    return x>=p?x-p:x;
}
int dec(int x,int y)
{
    x-=y;
    return x<0?x+p:x;
}
int main()
{
    //freopen("123.in","r",stdin);
    int x,c,mx,tem,cnt,ans;
    for (int i=2;i<=maxx;i++)
    {
        if (!vis[i]) prm[++tot]=i;
        for (int j=1;i*prm[j]<=maxx;j++)
        {
            vis[i*prm[j]]=1;
            if (i%prm[j]==0) break;
        }
    }
    scanf("%d",&n);
    for (int i=1;i<=n;i++)
    {
        scanf("%d",&x);
        for (int j=1;x>1&&j<=tot;j++)
            if (x%prm[j]==0)
            {
                c=0;
                while (x%prm[j]==0) x/=prm[j],c++;
                v[j].push_back(c);
            }
        if (x>1) f[++num]=x;
    }
    pw[0]=s[0]=ans=1;
    for (int i=1;i<=tot;i++) if (v[i].size())
    {
        mx=0;
        tem=1;
        for (vector<int>::iterator it=v[i].begin();it!=v[i].end();++it)
        {
            for (int j=mx+1;j<=*it;j++)
            {
                pw[j]=(LL)pw[j-1]*prm[i]%p;
                s[j]=inc(s[j-1],pw[j]);
            }
            mx=max(mx,*it);
            tem=(LL)tem*s[*it]%p;
        }
        tem=inc((LL)dec(tem,1)*(prm[i]-1)%p*pow(prm[i],p-2)%p,1);
        ans=(LL)ans*tem%p;
    }
    sort(f+1,f+num+1);
    cnt=1;
    for (int i=2;i<=num;i++)
        if (f[i]==f[i-1]) cnt++;
        else ans=(LL)ans*inc((LL)dec(pow(f[i-1]+1,cnt),1)*(f[i-1]-1)%p*pow(f[i-1],p-2)%p,1)%p,cnt=1;
    ans=(LL)ans*inc((LL)dec(pow(f[num]+1,cnt),1)*(f[num]-1)%p*pow(f[num],p-2)%p,1)%p;
    printf("%d\n",ans);
}