To kill all the bugs is always easier than to capture their brains. A map is drawn for you, with all the rooms marked by the amount of bugs inside, and the possibility of containing a brain. The cavern's structure is like a tree in such a way that there is one unique path leading to each room from the entrance. To finish the battle as soon as possible, you do not want to wait for the troopers to clear a room before advancing to the next one, instead you have to leave some troopers at each room passed to fight all the bugs inside. The troopers never re-enter a room where they have visited before.
A starship trooper can fight against 20 bugs. Since you do not have enough troopers, you can only take some of the rooms and let the nerve gas do the rest of the job. At the mean time, you should maximize the possibility of capturing a brain. To simplify the problem, just maximize the sum of all the possibilities of containing brains for the taken rooms. Making such a plan is a difficult job. You need the help of a computer.
The last test case is followed by two -1's.
5 10 50 10 40 10 40 20 65 30 70 30 1 2 1 3 2 4 2 5 1 1 20 7 -1 -1
50 7
答案:
#include <stdio.h>
#include <string>
#define max(a, b) a>b?a:b;
int n,m;
int path[101][101]; // 路径
int len[101],ls[101]; // 房间连接的其他房间个数
int room[101][2]; // 房间--0代表虫子数;1代表可能性
int ass[101][101]; // 在第i个房间放置j个士兵所带来的可能性
bool vist[101]; // 第i个房间是否已经被访问处理
int tmp[101]; // 放置第i个士兵所得到的可能性
void DP(int now)
{
int i,j,k,t;
if (ls[now] == 1)
{
vist[now] = true;
t = room[now][0] / 20;
if (room[now][0] % 20) t ++;
for (j=m-t;j>=0;j--) ass[now][j+t] = ass[now][j] + room[now][1];
for (t--;t>=0;t--) ass[now][t] = 0; // 重新置零
for (i=0;i<len[now];i++)
{
int next = path[now][i];
if (vist[next]) continue;
ls[next] --;
ass[now][0] = ass[next][0] = 0;
for (j=0;j<=m;j++) tmp[j] = max(ass[now][j],ass[next][j]);
for (j=0;j<=m;j++)
{
for (k=0;k<=m;k++)
{
if (j+k > m) break;
tmp[j+k] = max(tmp[j+k],ass[now][j]+ass[next][k]);
}
}
for (j=0;j<=m;j++) ass[next][j] = tmp[j];
if (ls[next] == 1) DP(next);
}
}
}
int solve() {
int i;
if (m == 0) return 0;
memset(ass,0,sizeof(ass));
memset(vist,0,sizeof(vist));
memcpy(ls,len,sizeof(len));
for (i = 1; i <= n; ++i)
DP(i);
int ret = 0;
for (i=1;i<=m;i++) ret = max(ret,ass[0][i]);
return ret;
}
int main() {
int i;
while (scanf("%d %d",&n,&m)==2)
{
if (n==-1 && m==-1) break;
for (i=1;i<=n;i++)
scanf("%d %d",&room[i][0],&room[i][1]);
for (i=0;i<=n;i++)
len[i] = 0;
for (i=0;i<n-1;i++)
{
int a,b;
scanf("%d %d",&a,&b);
path[a][len[a]++] = b;
path[b][len[b]++] = a;
}
path[1][ len[1]++ ] = 0;
printf("%d/n",solve());
}
return 0;
}
参考:
http://hi.baidu.com/%D3%FB%BB%B6%D0%A1%CF%B7/blog/item/b06b7d98dfd5df046f068c56.html