ZOJ 3375 Imperishable Night 【状压DP】

题目链接：http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3375

题意：有个人要去采集一些东西，有N(1 ≤ N ≤ 15)个采集点，每个采集点有两种东西（Item和Time），采集一份Item就会收获之前已经采集的Item（IV）的点数，且增加Time（TV），采集一份Time也类似，问能收获的最大的点数

分析：一看发现N很小，基本可以肯定这是需要状压，而且要使得收获最大就必须要采集完所有的采集点，由此可以类比TSP问题，所以可以确定一个框架。但是这个题中每个采集点所能收获的价值和之前的收获到的TV、IV和采集策略有关，问题似乎变得很复杂，因为采集策略又好像对TV和IV有影响，但是仔细研究发现其实不论是什么策略，没经过一个采集点TV和IV增加的量是一定的，TV+=ax，IV+=by。所以在采集一个新的点之前 ，IV和TV是一定的，可以先计算出来。所以只需要想一个获得最大点数的策略就好了。很明显的有两种策略，①先采集完一种再采集完另一种（具体更据a和b的大小确定先后顺序）②交替采集两种（具体更据x和y的大小确定先后顺序）。这样就直接推公式求最大值即可（没仔细推，是不是有一种策略明显优于另一种,不过我分别交了两次。。，发现策略①是最优的）

代码：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<vector>
#include<cmath>
#define N (1<<15)
using namespace std;
typedef long long ll;

ll dp[N],val[N][2];
struct node
{
    int a,b,x,y;
}s[20];
int n;

void init()
{
    int lim=(1<<n)-1;
    for(int i=0;i<=lim;i++)
    {
        int k=i,p=0;
        val[i][0]=val[i][1]=dp[i]=0;
        while(k)
        {
            if(k&1)
            {
                val[i][0]+=s[p].b*s[p].y;
                val[i][1]+=s[p].a*s[p].x;
            }
            k>>=1;p++;
        }
    }
}

ll getmax(int q,int p)
{
    ll r1,r2;
    if(s[p].a<=s[p].b)r1=val[q][1]*s[p].y+(val[q][0]+s[p].y*s[p].b)*s[p].x;
    else r1=val[q][0]*s[p].x+(val[q][1]+s[p].x*s[p].a)*s[p].y;
    if(s[p].x>=s[p].y)
    {
        int t=s[p].y;
        r2=(val[q][0]+val[q][1])*t+((t-1)*s[p].a*t+(t+1)*s[p].b*t)/2+(s[p].x-t)*(val[q][0]+s[p].y*s[p].b);
    }
    else
    {
        int t=s[p].x;
        r2=(val[q][0]+val[q][1])*t+((t-1)*s[p].a*t+(t+1)*s[p].b*t)/2+(t-s[p].x)*(val[q][1]+s[p].x*s[p].a);
    }
    return max(r1,r2);
}

int main()
{
    while(scanf("%d",&n)!=EOF)
    {
        for(int i=0;i<n;i++)scanf("%d%d%d%d",&s[i].a,&s[i].b,&s[i].x,&s[i].y);
        init();
        int lim=(1<<n)-1;
        for(int i=0;i<=lim;i++)
        {
            int k=i,p=0;
            while(k)
            {
                if(k&1)
                {
                    dp[i]=max(dp[i],dp[i^(1<<p)]+getmax(i^(1<<p),p));
                }
                k>>=1;p++;
            }
        }
        cout<<dp[lim]<<endl;
    }
    return 0;
}