HDU4341-Gold miner-分组DP

模拟黄金矿工这个游戏，给出每一个金子的位置和所需时间，计算在给定时间内最大收益。

刚看这道题以为金子的位置没什么用，直接DP就行，WA了一发终于明白如果金子和人共线的话只能按顺序抓。

这就是需要考虑先后关系问题。看了背包⑨讲之后以为是“有依赖关系的背包”，感觉解决方案很不明显，想不出来做法。

后来想到，可以把共线的金子按1，1+2，1+2+3。。。变成若干个，然后共线的金子组成一组。

显然这个问题就变成了在组内互斥的情况下的分组DP，背包⑨讲已经给出了伪代码。

预处理的时候使用了优先队列，map，乱搞了一下就好了。。。//好不容易自己写出来了一道不太水的题。。。

 1 #include <cstdio>
 2 #include <algorithm>
 3 #include <cstring>
 4 #include <queue>
 5 #include <stack>
 6 #include <map>
 7
 8 using namespace std;
 9
10 const int maxn = 200+10;
11
12 struct node
13 {
14     int x,y;
15     int w,v;
16     node(){x=0;y=0;w=0;v=0;}
17     node(int x,int y) {w=x;v=y;}
18     bool operator < (const node &b) const
19     {
20         return x*x+y*y > b.x*x+b.y*y;
21     }
22
23 }gold[maxn];
24
25 struct item
26 {
27     int w,v;
28     item(){w=0;v=0;}
29     item(int a,int b){w=a;v=b;}
30 };
31
32 int x[maxn],y[maxn],w[maxn],v[maxn],dp[40010];
33 int N,T;
34 vector <item> itm[maxn];
35
36 int main()
37 {
38     int cas = 1;
39     while(~scanf("%d%d",&N,&T))
40     {
41         map<pair<int,int> , int > mp;
42         priority_queue<node> pq[maxn];
43         int cnt_group = 0;
44         for(int i=0;i<maxn;i++) itm[i].clear();
45
46         for(int i=0,x,y;i<N;i++)
47         {
48             scanf("%d%d%d%d",&gold[i].x,&gold[i].y,&gold[i].w,&gold[i].v);
49             x = gold[i].x;y = gold[i].y;
50
51             if(mp.count(pair<int,int>( x/__gcd(x,y) , y/__gcd(x,y)) ))
52             {
53                 pq[mp[pair<int,int>( x/__gcd(x,y) , y/__gcd(x,y))]].push(gold[i]);
54             }
55             else
56             {
57                 mp[pair<int,int>( x/__gcd(x,y) , y/__gcd(x,y))]=cnt_group++;
58                 pq[mp[pair<int,int>( x/__gcd(x,y) , y/__gcd(x,y))]].push(gold[i]);
59             }
60         }
61
62         map<pair<int,int>,int >::iterator it;
63         node cur;
64         for(int i=0;i<cnt_group;i++)
65         {
66             cur = node(0,0);
67             while(!pq[i].empty())
68             {
69                 cur.w += pq[i].top().w;
70                 cur.v += pq[i].top().v;
71                 pq[i].pop();
72                 itm[i].push_back(item(cur.w,cur.v));
73                 //printf("grp:%d w:%d v:%d\n",i,cur.w,cur.v);
74             }
75         }
76         memset(dp,0,sizeof dp);
77
78         for(int gp=0;gp<cnt_group;gp++)
79         {
80             int n = itm[gp].size();
81
82             for(int t=T;t>=0;t--)
83             {
84                 for(int i=0;i<n;i++) if(t>=itm[gp][i].w)
85                 {
86                     //printf("t:%d w:%d\n",t,itm[gp][i].w);
87                     dp[t] = max(dp[t],dp[t-itm[gp][i].w] + itm[gp][i].v);
88                 }
89             }
90         //for(int i=0;i<T;i++)
91             //printf("%d ",dp[i]);
92         //printf("\n");
93         }
94
95         printf("Case %d: %d\n",cas++,dp[T]);
96     }
97 }