BestCoder Round #78 1002/hdu 5656 CA Loves GCD 动态规划

题意：求n个数的全部子集的gcd之和。

我们令dp[i][j]表示在前i个数中，选出若干个数使得它们的gcd为j的方案数，于是只需要枚举第i+1个数是否被选中来转移就可以了。

令第i+1个数为v，当考虑dp[i][j]的时候，我们令dp[i+1][j] += dp[i][j],dp[i+1][gcd(j,v)] += dp[i][j]。

复杂度O(N*MaxV) MaxV 为出现过的数的最大值。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#define MOD 100000007
#define N 1100
#define ll long long
using namespace std;

int gcd(int a,int b)
{
    return b==0?a:gcd(b,a%b);
}

int d[N][N],a[N];

int main()
{
    int T;
    cin>>T;
    while(T--)
    {
        int n;
        scanf("%d",&n);
        memset(d,0,sizeof(d));
        for(int i=0;i<n;i++)
        {
            scanf("%d",&a[i]);
            d[i][a[i]]++;
        }
        for(int i=0;i<n-1;i++)
            for(int j=1;j<=1000;j++)
                if(d[i][j])
                {
                    (d[i+1][j]+=d[i][j])%=MOD;
                    (d[i+1][gcd(j,a[i+1])]+=d[i][j])%=MOD;
                }
        ll ans=0;
        for(int i=1;i<=1000;i++)   (ans+=(ll)i*d[n-1][i])%=MOD;
        cout<<ans<<endl;
    }
}