CodeForces Round 428 div2 839D Winter is here(数学)(详解)

唉,要是哪天我能不看题解就做出D题来就好了.
感觉自己数学底子好差啊,我怎么都想不出来的问题,他们一句就带过了,应该是很显然的吧.
这题一开始我想到枚举gcd,但是我发现会有重复的情况,然后我就不知道怎么处理了.
这题我觉得最关键的就是这里了,假如很多数,他们都是2的倍数,但是显然,他们中的一些数的gcd不会还是2,但是(!!!!!!)他们的gcd,也仍然是2的倍数.如果你读到这里毫无感觉,那这就是我菜的地方了呀.
那么我们就先枚举gcd,然后假如有n个数是gcd的倍数,那么他们对答案的贡献假设为 gcd*x,那这个x是什么呢?
x = (1*C(n,1)+2*C(n,2)+….+n*C(n,n))(连我都很容易想到)
那问题来了,这个x要怎么样快速算出来呢?只要你知道组合数的定义(连我都知道),就能简单的把他化为 x = n * (2^(n-1)).
那我们现在来处理一下这n个里面某些组合在一起,他们的gcd不再是我们枚举的gcd的情况.咋办呢,从后往前就好了,只要减掉他们这种情况的x(没乘gcd呢)们,就行了.
举个例子把,因为我自己也不是很清楚.
有100个数他们都是2的倍数,但是其中有三个的gcd是18,那我们在计算gcd等于2时,我们的x中不能有刚好是他们三个的情况了(3*C(100,3)).
具体看代码把.

/*  xzppp  */
#include <iostream>
#include <vector>
#include <set>
#include <queue>
#include <map>
#include <algorithm>
#include <stdio.h>
#include <string.h>
#include <list>
#include <iomanip>
using namespace std;
#define FFF freopen("in.txt","r",stdin);freopen("out.txt","w",stdout);
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define mp make_pair
typedef long long  LL;
typedef unsigned long long ULL;
const int MAXN = 1e6;
const int  INF = 0x7fffffff;
const int MOD = 1e9+7;
LL a[MAXN+17],num[MAXN+17];
LL pow[MAXN/5+17];
int main()
{
    //FFF
    int n,mx = 0;
    cin>>n;
    for (int i = pow[0]= 1; i < 2*1e5+17; ++i)
        pow[i] = (pow[i-1]*2)%MOD;
    for (int i = 0; i < n; ++i)
    {
        int temp;
        scanf("%d",&temp);
        mx = max(temp,mx);
        a[temp]++;
    }
    LL ans = 0;
    for (int i = mx; i > 1 ; --i)
    {
        LL cnt = 0,tosub = 0;
        for (int j = i; j <= mx; j+=i)
        {
            cnt += a[j];
            tosub += num[j];
        }
        num[i] = (cnt*pow[cnt-1])%MOD-tosub;
        ans += (num[i]*i)%MOD;
    }
    cout<<ans%MOD<<endl;
}