[bzoj3512]DZY Loves Math IV

最新推荐文章于 2020-09-23 16:08:09 发布

YZH__12345

最新推荐文章于 2020-09-23 16:08:09 发布

阅读量549

点赞数

分类专栏：杜教筛

本文链接：https://blog.csdn.net/YZHYZHYZH3/article/details/78743235

版权

杜教筛专栏收录该内容

3 篇文章 0 订阅

订阅专栏

给定n,m，求这里写图片描述模10^9+7的值。
听说是杜教筛裸题…
可能是我太鶸了.感觉好难啊..
Orz w_yqts
Orz xudyh

(i, j) = g c d (i, j)

$(i,j)=gcd(i,j)$

s (i, j) = \sum i = 1 m φ (n i)

$s(i,j)=\sum_{i=1}^mφ(ni)$

a n s = \sum i = 1 n s (i, m)

$ans=\sum_{i=1}^n s(i,m)$
对于n%p=0,有φ(n*p)=φ(n)*p
所以只需要求

x=Πp1i $x=\Pi p_i^1$ 的s(x,m).

令d=(n,i)(从某处开始d为变量)
于是有i%d=0,φ(d*i)=φ(i)*d
∵n中最高次数为一次
∴(n/d,i*d)=1,即n/d与i*d互质
∴

φ (n * i) = φ (n d * i * d) = φ (n d) * φ (i) * d

$φ(n*i)=φ(\frac n d*i*d)=φ(\frac n d)*φ(i)*d$
∵

\sum d | n φ (d) = n

$\sum_{d|n} φ(d)=n$
∴

φ (n * i) = φ (n d) * φ (i) * \sum e | d φ (e) = φ (n d) * φ (i) * \sum e | d φ (d e) = φ (i) * \sum e | d φ (d e) * φ (n d) = φ (i) * \sum e | d φ (n e)

$φ(n*i)=φ(\frac n d)*φ(i)*\sum_{e|d} φ(e)=φ(\frac n d)*φ(i)*\sum_{e|d} φ(\frac d e)=φ(i)*\sum_{e|d} φ(\frac d e)* φ(\frac n d)= φ(i)*\sum_{e|d} φ(\frac n e)$

s (n, m) = \sum i = 1 n \sum d | n φ (n d) * \sum i = 1 ⌊ m d ⌋ φ (d i) = \sum d | n φ (n d) * s (d, ⌊ m d ⌋)

$s(n,m)=\sum_{i=1}^n\sum_{d|n} φ(\frac n d)*\sum_{i=1}^{\lfloor\frac m d\rfloor} φ(di)=\sum_{d|n}φ(\frac n d)*s(d,\lfloor\frac m d\rfloor)$
s(1,m)用杜教筛求
然后记忆化一下就好了.

#include <bits/stdc++.h>
using namespace std;
#define p 1000000007
#define ll long long
#define N 2500005
map <ll,ll> f,s;
int pn,pr[N],pri[N],flag[N];
ll phi[N],Ans[N],ph[N],ttt;
inline void init()
{
    pri[1]=flag[1]=phi[1]=ph[1]=1LL;
    for (int i=2;i<N;++i)
    {
        if (!flag[i]) pr[++pn]=pri[i]=i,phi[i]=i-1;
        ph[i]=phi[i];
        for (int j=1;j<=pn && pr[j]*i<N;++j)
        {
            flag[i*pr[j]]=1;
            if (i%pr[j]==0) {phi[i*pr[j]]=phi[i]*pr[j];pri[i*pr[j]]=pri[i];break;}
            phi[i*pr[j]]=phi[i]*(pr[j]-1);
            pri[i*pr[j]]=pri[i]*pr[j];
        }
    }
    for (int i=1;i<N;++i) phi[i]+=phi[i-1];
}
inline ll calc(int n)
{
    if (n<N) return phi[n];
    if (f[n]) return f[n];
    ll res=(ll)(n+1)*n>>1;
    for (int i=2,pos;i<=n;i=pos+1)
    {
        pos=n/(n/i);
        res-=calc(n/i)*(pos-i+1);
    }
    return f[n]=res;
}
inline int sum(int n,int m)
{
    if (!m) return 0;
    if (s[ttt*(n-1)+m]) return s[ttt*(n-1)+m];
    if (n==1) return calc(m)%p;
    ll ans=0LL;
    for (int i=1;i*i<=n;++i)
    if (n%i==0)
    {
        ans+=ph[n/i]*sum(i,m/i)%p;
        if (i*i!=n) ans+=ph[i]*sum(n/i,m/(n/i))%p;
    }
    return s[ttt*(n-1)+m]=ans%p;
}
int main()
{
    init();
    ll n,m;
    cin>>n>>m;ttt=(ll)m;
    ll ans=0LL;
    for (int i=1;i<=n;++i)
    {
        if (pri[i]!=i) Ans[i]=Ans[pri[i]]*(i/pri[i])%p;else
        Ans[i]=sum(i,m);
        ans=(ans+Ans[i])%p;
    }
    cout<<ans<<endl;
}