| ||||||||||
| Online Judge | Problem Set | Authors | Online Contests | User | ||||||
|---|---|---|---|---|---|---|---|---|---|---|
| Web Board Home Page F.A.Qs Statistical Charts | Current Contest Past Contests Scheduled Contests Award Contest | |||||||||
为了这道题贡献了32个WA啊。。。。给大家一点提示不知怎么,这么一道简单的TreeDP弄得我这么狼狈。。。
关于这道题,网上有两种解法,一种是双DP,过程叫清楚,也适合于对于TreeDP不是太了解或者自由树转有根树理解不太透彻的人来做,抑或对背包问题了解的很牛的人,过双DP也很简单。。其实双DP也很具有启发意义,是很精巧的算法。。。但是个人认为第二种算法更体现了TreeDP的精髓。。
这儿介绍一下第二种方法
先把自由树转化成有根树(直接转成孩子兄弟表示法),网上我找到的唯一一个用这种方法做的人,在树的转化这个环节用了100行(Orz,膜拜),其实一个遍历就可以了。孩子兄弟表示法到处都有教程,这儿就不多赘述)
void Zh(int p)
{
f[p]=1;
for(int i=0;i<d[p];i++)
if (!f[next[p][i]])
{
bro[next[p][i]]=son[p];
son[p]=next[p][i];
Zh(next[p][i]);
}
}
然后直接TreeDP
用maxt[0][p][m]表示p的右边兄弟和下面儿子(这个有点XX吧,但是相信大家能理解)消耗m步还要回起点p,最多能吃的个数。。
用maxt[1][p][m]表示p的右边兄弟和下面儿子,消耗m步不需回起点p,最多能吃的个数。。
对于每一个节点p,先考虑无视p的情况(p取不到),直接把“走路的权利”传给自己右边兄弟的情况。
接下来我们考虑的情况,p一定能取道。
分“只往下扩展孩子”和“只往右扩展兄弟”两种情况。
最后,考虑既扩展孩子又扩展兄弟的情况。。。
其中最复杂的就算是“走路的权利”的问题了。。。
我们分maxt[0]和maxt[1]两种情况讨论
maxt[0]时,一定要回到原点,那么很简单,扩展孩子和扩展右边兄弟的时候,“走路的权利”都要-2(想想为什么?)
maxt[0]时,不需要回到原点,那么扩展右边兄弟的时候,“走路的权利”-2,而扩展孩子的时候,只需要-1(想想为什么?)
这两个问题,如果你画一下图,也许能很快明白。。。。
如果你全部明白了,恭喜你,这道题AC将会在30分钟内诞生。。。
祝你成功!
PS:
以下是小弟的代码,为了看的清楚,我把程序写的简洁,但也损失了很多时间,当然对于p==0和m==0的时候的剪枝,其实很可以做的更好。
祝大家AC愉快!
Source Code
Problem: 2486 User: forsona
Memory: 620K Time: 422MS
Language: G++ Result: Accepted
Source Code
#include"iostream"
#define INF -99999999
using namespace std;
int d[110],next[110][110],son[110],bro[110],c[110],f[110],n,k,maxt[2][110][210],r;
void Zh(int p)
{
f[p]=1;
for(int i=0;i<d[p];i++)
if (!f[next[p][i]])
{
bro[next[p][i]]=son[p];
son[p]=next[p][i];
Zh(next[p][i]);
}
memset(maxt,-1,sizeof(maxt));
}
void Inp(void)
{
int i,j,k;
memset(d,0,sizeof(d));
for(i=1;i<=n;i++)
scanf("%d",&c[i]);
for(i=0;i<n-1;i++)
{
scanf("%d%d",&j,&k);
next[j][d[j]++]=k;
next[k][d[k]++]=j;
}
memset(bro,0,sizeof(bro));
memset(son,0,sizeof(son));
memset(f,0,sizeof(f));
Zh(1);
}
int Loop(int b,int p,int m)
//b==0 back to root b==1 no need to root
{
int i,mt;
if (maxt[b][p][m]!=-1)
return maxt[b][p][m];
if (!p)
return maxt[b][p][m]=0;
if (m==0)
return maxt[b][p][m]=max(Loop(b,bro[p],m),c[p]);
if (b==0)
{
mt=Loop(0,bro[p],m);
if (m>=2)
{
mt=max(mt,c[p]+Loop(0,bro[p],m-2));
mt=max(mt,c[p]+Loop(0,son[p],m-2));
}
for(i=0;i<=m-4;i++)
mt=max(mt,Loop(0,bro[p],i)+Loop(0,son[p],m-4-i)+c[p]);
return maxt[b][p][m]=mt;
}
else
{
mt=Loop(1,bro[p],m);
mt=max(mt,c[p]+Loop(1,son[p],m-1));
if (m>=2)
mt=max(mt,c[p]+Loop(1,bro[p],m-2));
for(i=0;i<=m-4;i++)
mt=max(mt,Loop(1,bro[p],i)+Loop(0,son[p],m-4-i)+c[p]);
for(i=0;i<=m-3;i++)
mt=max(mt,Loop(0,bro[p],i)+Loop(1,son[p],m-3-i)+c[p]);
return maxt[b][p][m]=mt;
}
}
int main()
{
while(scanf("%d%d",&n,&k)!=EOF)
{
Inp();
r=Loop(1,1,k);
printf("%d\n",r);
}
return 0;
}
Followed by:
Post your reply here: |
All Rights Reserved 2003-2013 Ying Fuchen,Xu Pengcheng,Xie Di
Any problem, Please Contact Administrator