HDU 5656 CA Loves GCD

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5656

题意：给n个数，求出所有非空子集的gcd之和。题目中所有出现的数的范围都不超过1000。

思路：由于范围很小，我们用f[i][j]表示前i个数构成gcd为j的子集的方法数，有递推式f[i][gcd(a[i],j)] += f[i-1][j](选ai)，f[i][j] +=f[i-1][j](不选ai) f[i][a[i]]++ 单独选ai。

注意计算的时候取模，还有需要先预处理出所有的gcd，不然会超时。

#include <cstdio>
#include <cmath>
#include <cstring>
#include <string>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <stack>
#include <map>
#include <set>
#include <vector>
#include <sstream>
#include <queue>
#include <utility>
using namespace std;

#define rep(i,j,k) for (int i=j;i<=k;i++)
#define Rrep(i,j,k) for (int i=j;i>=k;i--)

#define Clean(x,y) memset(x,y,sizeof(x))
#define LL long long
#define ULL unsigned long long
#define inf 0x7fffffff
#define mod 100000007
const int maxn = 1009;

LL f[maxn][maxn];
int num[maxn];
int g[maxn][maxn];
int T;
int n;

int gcd( int a , int b )
{
    return (b==0?a:gcd(b,a%b));
}
void init()
{
    rep(i,1,1000)
        rep(j,i,1000)
            g[i][j] = g[j][i] = gcd(i,j);
}
int main()
{
    init();
    cin>>T;
    while( T-- )
    {
        int uplim = 0;
        scanf("%d",&n);
        rep(i,1,n)
        {
            scanf("%d",&num[i]);
            uplim = max( uplim , num[i] );
        }
        Clean(f,0);
        rep(i,1,n)  // f[i][j]  前i个数 gcd为j
        {
            f[i][num[i]]++;
            rep(j,1,uplim)
                f[i][j]+=f[i-1][j];
            rep(j,1,uplim)
            {
                f[i][ g[j][num[i]] ] += f[i-1][j];
            }
            rep(j,1,uplim) f[i][j]%=mod;
        }
        LL ans = 0;
        rep(i,1,uplim)
        ans = ( ans + f[n][i]*i ) % mod;
        printf("%I64d\n",ans);
    }
    return 0;
}