[bzoj3560] DZY Loves Math V

最新推荐文章于 2019-07-30 19:02:00 发布

原创最新推荐文章于 2019-07-30 19:02:00 发布 · 610 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

数论专栏收录该内容

8 篇文章

订阅专栏

本文介绍了一道数学竞赛题目，通过运用欧拉函数的性质来解决该问题。利用欧拉函数的积性特点，文章详细解释了如何将原问题转化为质数处理，并给出了具体的算法实现过程。

题目大意

给定n个正整数的数组，求 $\sum_{i1|a1} \sum_{i2|a2}...\sum_{in|an} \phi(i1i2i3...in)$ 模 $10^9+7$ 的值。

n≤ $10^5$
ai≤ $10^7$

分析

计算这题要用到欧拉函数的性质。
1. 欧拉函数是积性函数
2. 如果n= $p^k$ （p是质数）， $\phi(n)=(p-1)p^{k-1}=p^k*\frac{p-1}{p}$

首先所有质数分开处理。
设当前质数为p，对于第i个数，假设它分解质因数后p的次数为ai，那么p的答案就是

[(1 + p 1 + . . . + p a 1) (1 + p 1 + . . . + p a 2) . . . (1 + p 1 + . . . + p a n) - 1] p - 1 p + 1

$[(1+p^1+...+p^{a1})(1+p^1+...+p^{a2})...(1+p^1+...+p^{an})-1]\frac{p-1}{p}+1$
式子的意义就是每个数分别取0,1…ai个p，减掉的1、最后加上的1是全部取1的。
最后把每个质数的答案乘起来就好了。

#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;

const int N=1e7,mo=1e9+7;

typedef long long LL;

int T,tot,p[N],f[N+5],n,ans,s[N+5],t[N],id[N+5],a[N];

bool bz[N+5];

char c;

int read()
{
    for (c=getchar();c<'0' || c>'9';c=getchar());
    int x=c-48;
    for (c=getchar();c>='0' && c<='9';c=getchar()) x=x*10+c-48;
    return x;
}

int quick(int x,int y)
{
    if (!y) return 1;
    int ans=quick(x,y>>1);
    ans=(LL)ans*ans%mo;
    if (y&1) ans=(LL)ans*x%mo;
    return ans;
}

int main()
{
    for (int i=2;i<=N;i++)
    {
        if (!bz[i])
        {
            p[id[i]=tot++]=i;
        }
        for (int j=0;j<tot && i*p[j]<=N;j++)
        {
            bz[i*p[j]]=1;
            if (i%p[j]==0) break;
        }
    }
    for (int i=0;i<tot;i++) s[i]=1;
    n=read();
    for (int i=1;i<=n;i++) a[i]=read();
    for (int i=1;i<=n;i++)
    {
        for (int j=0;p[j]*p[j]<=a[i];j++) if (a[i]%p[j]==0)
        {
            int sum=1;
            for (int k=p[j];a[i]%p[j]==0;a[i]/=p[j],k*=p[j]) sum+=k;
            s[j]=(LL)s[j]*sum%mo;
        }
        if (a[i]>1) s[id[a[i]]]=(LL)s[id[a[i]]]*(1+a[i])%mo;
    }
    ans=1;
    for (int i=0;i<tot;i++)
    {
        s[i]=(s[i]+mo-1)%mo;
        ans=(LL)ans*((LL)s[i]*(p[i]-1)%mo*quick(p[i],mo-2)%mo+1)%mo;
    }
    printf("%d\n",ans);
    return 0;
}