839D

839D - Winter is here

题意：

题解：

那个式子我都没想着去化简~太弱了=_=

参考的dalao的代码~

#include <bits/stdc++.h>
using namespace std;
#define ll long long
const int maxn=200010;
const int mod=1e9+7;

int cnt1[maxn*5],cnt2[maxn*5];
ll dp[maxn*5];
int base[maxn];

int main(){
    base[0]=1;
    for(int i=1;i<maxn;i++) base[i]=base[i-1]*2%mod;  //预处理
    int n;
    scanf("%d",&n);
    int x;
    for(int i=0;i<n;i++){
        scanf("%d",&x);
        cnt1[x]++;
    }
    for(int i=2;i<=1000000;i++){
        for(int j=i;j<=1000000;j+=i){
            cnt2[i]+=cnt1[j];
        }
    }
    for(int i=1000000;i>=2;i--){
        if(cnt2[i]==0) continue;
        dp[i]=(1LL*base[cnt2[i]-1]*cnt2[i])%mod;
        for(int j=i*2;j<=1000000;j+=i) dp[i]=(dp[i]+mod-dp[j])%mod;
    }
    int ans=0;
    for(int i=2;i<=1000000;i++) ans=(ans+i*dp[i])%mod;
    printf("%d
",ans);
    return 0;
}

下面这个更优一些（不完整

int N;
ll A[202020];
ll mo=1000000007;

ll num[1010100];
ll p[1010100];

ll modpow(ll a, ll n = mo-2) {
    ll r=1;
    while(n) r=r*((n%2)?a:1)%mo,a=a*a%mo,n>>=1;
    return r;
}

void solve() {
    int i,j,k,l,r,x,y; string s;
    
    cin>>N;
    FOR(i,N) {
        cin>>A[i];
        for(x=1;x*x<=A[i];x++) if(A[i]%x==0) {
            num[x]++;
            if(x*x!=A[i]) num[A[i]/x]++;
        }
    }
    
    ll tot=0;
    for(i=1000000;i>=2;i--) {
        if(num[i]) {
            p[i]=num[i]*modpow(2,num[i]-1)%mo;
            for(j=i*2;j<=1000000;j+=i) p[i]-=p[j];
            p[i]=(p[i]%mo+mo)%mo;
            
            (tot+=p[i]*i)%=mo;
        }
    }
    cout<<tot<<endl;
}

另一种解法，还没仔细看~