【USACO】2014 Feb Cow Decathlon 全能奶牛

Cow Decathlon 全能奶牛

---

• Description

约翰有N头奶牛，组成了一直队伍参加全能比赛。比赛一共有N项，每头奶牛必须参加一项比赛，每项比赛也必须有一头奶牛参加。任何一头奶牛可以胜任任何一项比赛，但得分不一样。如果第i头奶牛参加第j项比赛，在比赛结束的时候，可以为团体总分增加$S_{i,j}$ 。
比赛是按照顺序依次进行的。除了上述获得分数的方法之外，还有B种奖励分。获得奖励的方法是在前几项比赛里获得足够的分数。具体来说，第i项奖励会在第$K_i$项比赛结束的时候检查，如果 当时的总分大于或等于$P_i$，奶牛们就可以立即获得额外的$A_i$ 分。如果有多项奖励在同一时刻检查，奶牛可以自由安排检查和加分的顺序。请问约翰应该如何安排奶牛参加比赛，才能让它们获得最高的 分数？

• Input Format

第一行：两个整数N和B，1≤N≤20 , 1≤B≤20
第二行到第B+1行：第i+1行有三个整数$K_i$，$P_i$和$A _i$，1≤$K_i$≤N ,1≤$P_i$≤40000 , 1≤$A _i$≤1000
第B+2行到第B+N+1行：第i+B+1 行有N个整数，代表$S_{i,1}$到$S_{i,N}$，对每个1≤j≤N , 1≤$S_{i,j}$≤1000

• Output Format

单个整数：表示奶牛们可以获得的最大得分

• Sample Input

3 1
2 7 6
5 1 7
2 2 4
4 2 1

• Sample Output

17

• Hint

第一项比赛由第一头奶牛参加，第二项比赛由第三头奶牛参加，第三项比赛由第二头奶牛参加

---

• 分析

这道题第一眼看上去就像是动规的状压，用$F[i][S]$ 表示前i场比赛奶牛的参加情况。虽然时间好像很卡，但还是很高兴的打出来。
结果发现只能过三个点，太坑了！！！
仔细一想发现这个状态的第一维是冗余的，因为S里有几个1，i就是几。所以我们可以去掉第一维。然后用一个队列进行递推。

---

30%

#include <queue>
#include <stack>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int B[21][21],F[21][(1<<20)+2],start[22],n,b;
struct Data{int k,p,a;}A[21];
inline bool cmp(const Data&A,const Data&B){return A.k==B.k?A.p<B.p:A.k<B.k;}
int main(){
    freopen("in.txt","r",stdin);
    freopen("out.txt","w",stdout);
    scanf("%d%d",&n,&b);
    for (int i=1;i<=b;i++) scanf("%d%d%d",&A[i].k,&A[i].p,&A[i].a);
    for (int i=1;i<=n;i++) 
        for (int j=1;j<=n;j++) scanf("%d",&B[i][j]);
    sort(A+1,A+1+b,cmp);
    for (int i=1;i<=b;i++)
        if (!start[A[i].k]) start[A[i].k]=i;
    A[0].k=1<<30; A[b+1].p=1<<30;
    memset(F,128,sizeof(F)); int Inf=F[0][0];
    for (int i=0;i<n;i++) F[1][1<<i]=B[i+1][1];
    int State=1<<n;
    for (int i=1;i<n;i++)
        for (int j=0;j<State;j++){
            if (F[i][j]==Inf) continue;
            for (int k=start[i];k!=start[i+1] && F[i][j]>=A[k].p;k++) F[i][j]+=A[k].a;
            for (int k=0;k<n;k++)
                if (j^(1<<k)) F[i+1][j|(1<<k)]=max(F[i][j]+B[k+1][i+1],F[i+1][j|(1<<k)]);
        }
    for (int i=start[n];F[i][State-1]>=A[i].p;i++) F[n][State-1]+=A[i].a;
    printf("%d",F[n][State-1]);
    fclose(stdin); fclose(stdout);
    return 0;
}

100%

#include <queue>
#include <stack>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
queue <int> Q;
int B[21][21],F[(1<<20)+2],start[22],n,b;
struct Data{int k,p,a;}A[21];
inline bool cmp(const Data&A,const Data&B){return A.k==B.k?A.p<B.p:A.k<B.k;}
int main(){
    freopen("in.txt","r",stdin);
    freopen("out.txt","w",stdout);
    scanf("%d%d",&n,&b);
    for (int i=1;i<=b;i++) scanf("%d%d%d",&A[i].k,&A[i].p,&A[i].a);
    for (int i=1;i<=n;i++) 
        for (int j=1;j<=n;j++) scanf("%d",&B[i][j]);
    sort(A+1,A+1+b,cmp);
    for (int i=1;i<=b;i++)
        if (!start[A[i].k]) start[A[i].k]=i;
    A[0].k=1<<30; A[b+1].p=1<<30;
    memset(F,128,sizeof(F)); int Inf=F[0];
    for (int i=0;i<n;i++) F[1<<i]=B[i+1][1],Q.push(1<<i);
    int State=1<<n;
    for (int i=0,j;!Q.empty();i=0){
        j=Q.front(); Q.pop();
        for (int k=j;k;k^=(k&(-k))) i++;
        for (int k=start[i];i==A[k].k && F[j]>=A[k].p;k++) F[j]+=A[k].a;
        for (int k=0;k<n;k++)
            if (!(j&(1<<k))){
                if (F[j|(1<<k)]==Inf) Q.push(j|(1<<k));
                F[j|(1<<k)]=max(F[j]+B[k+1][i+1],F[j|(1<<k)]);
            }
    }
    printf("%d",F[State-1]);
    fclose(stdin); fclose(stdout);
    return 0;
}