hdu 5445 Food Problem 分组背包，二进制拆分

题意：要得到至少p的能量，有n种甜品，每种甜品有能量值t，体积u，数量v，运送这些甜品的卡车有m种，每种有容积x，花费y，数量z。求运送的甜品至少能量达到p的最小花费。p<=50000，最小花费<=50000，n,m<=200，其他<=100，每块甜品可以分开装在不同的车上，但必须把这块甜品全部装载。

分析：首先可以算出的到p能量所需的最小体积，在算出得到这个最小体积的最小花费。最小体积可能很大，由于最小花费<=50000，可以dp出每个花费得到的最大容积，去满足的最小花费。物品总数过多，可以把每种物品按二进制拆分为log种物品。

#include<iostream>
#include<cstring>
#include<cstdio>
#include<ostream>
#include<istream>
#include<algorithm>
#include<queue>
#include<string>
#include<cmath>
#include<set>
#include<map>
#include<stack>
#include<vector>
#define fi first
#define se second
#define ll long long
#define pii pair<int,int>
#define inf (1<<30)
#define eps 1e-8
#define pb push_back
using namespace std;
const int maxn=110005;
int n,m,p;
int t,u,v;
int x,y,z;
struct Cot
{
    int w;
    int o;
};
Cot cot1[20000];
Cot cot2[20000];
int cnt1;
int cnt2;
int dp1[maxn];
int dp2[maxn];
int main()
{
    int T;
    scanf("%d",&T);
    while(T--) {
        scanf("%d%d%d",&n,&m,&p);
        cnt1=cnt2=0;
        for(int i=0;i<n;i++) {
            scanf("%d%d%d",&t,&u,&v);
            int k=1;
            while(k<v) {
                cot1[cnt1].w=k*t;
                cot1[cnt1++].o=u*k;
                v-=k;
                k*=2;
            }
            if(v) {
                cot1[cnt1].w=v*t;
                cot1[cnt1++].o=u*v;
            }
        }
        for(int i=0;i<m;i++) {
            scanf("%d%d%d",&x,&y,&z);
            int k=1;
            while(k<z) {
                cot2[cnt2].w=y*k;
                cot2[cnt2++].o=x*k;
                z-=k;
                k*=2;
            }
            if(z) {
                cot2[cnt2].w=y*z;
                cot2[cnt2++].o=x*z;
            }
        }
        dp1[0]=0;
        for(int i=1;i<=p+100;i++)
            dp1[i]=inf;
        int minv=inf;
        for(int i=0;i<cnt1;i++) {
            for(int j=p+100;j>=0;j--) {
                if(j-cot1[i].w>=0)
                    dp1[j]=min(dp1[j],dp1[j-cot1[i].w]+cot1[i].o);
                if(j>=p) {
                    minv=min(minv,dp1[j]);
                }
            }
        }
        memset(dp2,0,sizeof(dp2));
        int ans=inf;
        for(int i=0;i<cnt2;i++) {
            for(int j=50000;j>=0;j--) {
                if(j>=cot2[i].w) {
                    dp2[j]=max(dp2[j],dp2[j-cot2[i].w]+cot2[i].o);
                }
                if(dp2[j]>=minv) {
                    ans=min(ans,j);
                }
            }
        }
        if(ans!=inf)
            printf("%d\n",ans);
        else
            printf("TAT\n");
    }
    return 0;
}