“玲珑杯”ACM比赛 Round #12【dp】

题目：http://www.ifrog.cc/acm/contest/1014
题解：http://www.ifrog.cc/acm/solution/17

A签到题。
B。
题意：
给一个长度是n的字符串，有’0’,’1’,’?’，现在要在把?’变成0或1，变0的代价是c0,变1的代价是c1，而且不能有连续的n0个’0’或者n1个1，输入保证有解，问最小代价？
分析：
刚开始以为dp，想了想贪心即可。如果c0< c1，那么？尽量变成0。要注意的是判断当前？是否变0，要判断{ 前面已经有了多少个连续的0+后面的连续的0+1}< n0。
c0>=c1的情况也是这样。

C
题意：
有n种药，每种药有Wi单位，和一个特征值Ai，现在要从每种药取出Ui（Ui>0）单位，满足

问方案数

Hint：
You can verify the sample by brute force.
1≤N≤10^3
1≤Ai,Wi≤10^5
1≤G≤10^3

分析：
官方题解：
设(F[i][j])表示已经考虑前(i)种药，并且前 i 种药的(∑ix=1gcd(Ax,Ux)=j) 的不同的, .U1…i.的种数。
转移显然是(F[i][j]=∑Wiu=1F[i−1][j−gcd(u,Ai)])，然而这不够优。
对于每一对((Ai,Wi))
，预处理出对于所有可能的Ui，所有可能出现的(gcd)的值以及各个值对应的(Ui)个数，然后再开始在背包上转移，整道题的时间复杂度是(O(∑i=1Nd2(Ai)+G∗∑i=1Nd(Ai)))。

dp方程还是比较容易想的，比赛的时候自己只想到O(n^3)的转移。预处理gcd差不多也是O（1e9），时间给了5s，冲一冲没准能过 果然过不了QAQ。
超时代码：

#include<bits/stdc++.h>
using namespace std;
int gcd(int a,int b){return b==0?a:gcd(b,a%b);}
typedef pair<int,int>pii;
typedef long long ll;
const int mod=1234321237;
const int N=1e3+7;
int n,g,a[N],w[N],f[N][N];
ll d[N][N];
int main() {
   freopen("in.txt","r",stdin);
    scanf("%d%d",&n,&g);
    for(int i=0;i<n;i++)scanf("%d",a+i);
    for(int i=0;i<n;i++)scanf("%d",w+i);
    if(g<n){
        puts("0");return 0;
    }

    memset(f,0,sizeof(f));
    for(int i=0;i<n;i++){
        for(int j=1;j<=w[i];j++){
            int tt=gcd(a[i],j);
            if(tt>g)continue;
            f[i+1][tt]++;
        }
    }
    d[0][0]=1;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=g;j++){
            for(int k=1;k<=j;k++){
                if(f[i][k])
                d[i][j]=(d[i][j]+d[i-1][j-k]*f[i][k]%mod)%mod;
            }
        }
    }
    printf("%d\n",d[n][g]);
    return 0;
}

赛后优化了一下，勉强过了。

#include<bits/stdc++.h>
using namespace std;
int gcd(int a,int b){return b==0?a:gcd(b,a%b);}
typedef pair<int,int>pii;
typedef long long ll;
const int mod=1234321237;
const int N=1e3+2;
int n,g,a[N],w[N],f[N][N];
ll d[N][N];
int main() {
   //freopen("in.txt","r",stdin);
    scanf("%d%d",&n,&g);
    for(int i=0;i<n;i++)scanf("%d",a+i);
    for(int i=0;i<n;i++)scanf("%d",w+i);
    if(g<n){
        puts("0");return 0;
    }
    g-=n; //每一种药肯定会拿一单位
    memset(f,0,sizeof(f));
    for(int i=0;i<n;i++){
        for(int j=1;j<=w[i];j++){
            int tt=gcd(a[i],j);
            if(tt>g+1)continue;
            f[i+1][tt-1]++;
        }
    }
    d[0][0]=1;
    for(int i=1;i<=n;i++){
        for(int j=0;j<=g;j++){
            for(int k=0;k<=j;k++){
                if(f[i][k])
                d[i][j]=(d[i][j]+d[i-1][j-k]*f[i][k]%mod)%mod;
            }
        }
    }
    printf("%d\n",d[n][g]);
    return 0;
}

看到网上有人把标程贴出来了，不知道有没有经过出题人同意惹
http://www.cnblogs.com/ECJTUACM-873284962/p/6574397.html

继续优化就要优化求gcd了，学习一下~~

#include <bits/stdc++.h>
const int MOD = 1234321237;
int F[100001], N, G, a[1000], w[1000];

void DP(int x, int y)
{
    std::vector < int > Div;
    for (int i = 1; i * i <= x; i++)
        if (x % i == 0)
        {
            Div.push_back(i);
            if (i * i < x)
                Div.push_back(x / i);
        }
    std::sort(Div.begin(), Div.end());
    int L = Div.size();
    std::vector < int > Use(L, 0);
    for (int i = L - 1; ~i; i--)
    {
        Use[i] = y / Div[i];
        for (int j = i + 1; j < L; j++)
            if (Div[j] % Div[i] == 0)
                Use[i] -= Use[j];
    }
    for (int i = G; ~i; i--)
    {
        F[i] = 0;
        for (int j = 0; j < L && Div[j] <= i; j++)
            F[i] = (F[i] + (long long)F[i - Div[j]] * Use[j]) % MOD;
    }
}
int main()
{
    scanf("%d%d", &N, &G);
    for (int i = 0; i < N; i++)
        scanf("%d", a + i);
    for (int i = 0; i < N; i++)
        scanf("%d", w + i);
    F[0] = 1;
    for (int i = 0; i < N; i++)
        DP(a[i], w[i]);
    printf("%d\n", F[G]);
    return 0;
}

【吐槽】
这场过了两题，因为有事比赛很晚才开始打的，第三题差一点QAQ，然后我却随机到抱枕了，难道命运使然？还是挺开心的。