DP+预处理优化——玲珑学院OJ 1108

• 题目连接：http://www.ifrog.cc/acm/problem/1108

• 分析：用DP[i][j]表示前i个数达到和j的u的取值方法数，那么我们每次更新：
  $DP[i][j] = \sum ^{w[i]}_{u = 0} DP[i-1][j-gcd(u,A[i])]$ 即可，但是N的范围是1E3，加上gcd的时间，算法复杂度为$O(N^3log(N))$ 妥妥的超时。所以我们还需要优化一下。

• 优化：我们发现在更新 $DP[i][j]$ 的时候，有很多$DP[i-1][j-gcd(u,A[i])]$的值是相同的，我们是可以一次性加进去的。即我们可以预处理出每一个 $gcd$ 的值和这个值的次数：

• 预处理方法：计算出 $A[i]$ 的每一个约数，然后用 $W[i]$ 除以这个约数，得到它可能出现的次数，如果约数之间可以互相除尽，那么就有重复，所以我们从大往小计算约数出现的次数，如果大的约数能够除尽小的约数，那么小的约数里需要减去大的约数出现的次数（这个应该是比较好理解的，写2个数算算就出来了）。

• 新的DP方程：

for(int j=G;j>=0;j--)
{
    DP[j] = 0;
    for( int k=0;k<GCD.size() && GCD[k]<=j;k++)
    {
        DP[j] += DP[j-GCD[k]]*Times[k]; //GCD为第几个约数，Times为出现的次数
        if(DP[j]>mod) DP[j]%=mod;
    }
}

• AC代码：

/*************************************************************************
    > File Name: C.cpp
    > Author: Akira
    > Mail: qaq.febr2.qaq@gmail.com
    > Created Time: 2017年03月18日 星期六 16时35分51秒
 ************************************************************************/

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<cstdlib>
#include<algorithm>
#include<queue>
#include<stack>
#include<map>
#include<cmath>
#include<vector>
#include<set>
#include<list>
#include<ctime>
typedef long long LL;
typedef unsigned long long ULL;
typedef long double LD;
#define MST(a,b) memset(a,b,sizeof(a))
#define CLR(a) MST(a,0)
#define Sqr(a) ((a)*(a))

using namespace std;

#define MaxN 1001
#define MaxM MaxN*10
#define INF 0x3f3f3f3f
#define PI 3.1415926535897932384626
const int mod = 1234321237;
const int eps = 1e-8;
#define bug cout << 88888888 << endl;

int N,G;
int A[MaxN];
int W[MaxN];

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

LL DP[MaxN*100];
vector<int>::iterator it;
void solve()
{
    DP[0] = 1; 
    for(int i=1;i<=N;i++)
    {
        vector<int> GCD;
        for(int k=1;k*k<=A[i];k++)
        {
            if(A[i]%k==0)
            {
                GCD.push_back(k);
                if(k*k<A[i]) GCD.push_back(A[i]/k);
            }
        }
        sort(GCD.begin(),GCD.end());
        vector<int> Times (GCD.size(), 0);
        for(int k=GCD.size()-1;k>=0;k--)
        {
            Times[k] = W[i]/GCD[k];
            for(int l=k+1;l<GCD.size();l++)
            {
                if(GCD[l] % GCD[k] == 0) Times[k] -= Times[l];
            }
        }
        for(int j=G;j>=0;j--)
        {
            DP[j] = 0;
            for( int k=0;k<GCD.size() && GCD[k]<=j;k++)
            {
                DP[j] += DP[j-GCD[k]]*Times[k];
                if(DP[j]>mod) DP[j]%=mod;
            }
        }
    }
    printf("%lld\n", DP[G]);
}

int main()
{
    scanf("%d%d", &N, &G);
    for(int i=1;i<=N;i++) scanf("%d", &A[i]);
    for(int i=1;i<=N;i++) scanf("%d", &W[i]);
    solve();
}