显然,今天(好吧已经算是昨天了)的一天还算是有点意义的。
首先,不晓得应该是感谢还是鄙视军训的时候肚子疼的我~真的不是故意的。。。。
这两天主要搞的是树形DP和状态压缩DP。不过, LG又打算喊 我 和 猪鱼 做一些无聊的搜索题。。。
好吧,还是按顺序说::::::::(冒号是毫无意义的)
第一题,NOIP的加分二叉树
#include<bits/stdc++.h>
using namespace std;
int dp[30][30],a[30],rt[30][30];
int score(int l,int r){
if(dp[l][r])return dp[l][r];
if(l>r)return dp[l][r]=1;
if(l==r){rt[l][l]=l;return dp[l][l]=a[l];}
int maxx=-1;
for(int i=l;i<=r;i++){
int s=score(l,i-1)*score(i+1,r)+a[i];
if(s>maxx){
maxx=s;rt[l][r]=i;
}
}
return dp[l][r]=maxx;
}
void Pre(int l,int r){
if(l>r)return;
printf("%d ",rt[l][r]);
Pre(l,rt[l][r]-1);
Pre(rt[l][r]+1,r);
}
int main(){
int i,j,k,m,n;
scanf("%d",&n);
for(i=1;i<=n;i++)scanf("%d",a+i);
int ans=score(1,n);
printf("%d\n",ans);
Pre(1,n);
return 0;
}
没什么,显然是一道很右脑的题目,没什么多说的。
第二道POJ2486 Apple Tree
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<vector>
using namespace std;
vector <int> old[210],now[210];
int dp[2][210][210],n,m,k,val[210];
bool p[210];
void updata(int x){
p[x]=1;
for(int i=0;i<old[x].size();i++){
if(p[old[x][i]])continue;
now[x].push_back(old[x][i]);
updata(old[x][i]);
}
}
void dfs(int x){
for(int i=0;i<=k;i++)dp[0][x][i]=dp[1][x][i]=val[x];
for(int i=0;i<now[x].size();i++){
int y=now[x][i];
dfs(y);
for(int j=k;j>=0;j--){
for(int jj=0;jj<=j;jj++){
dp[0][x][j+1]=max(dp[0][x][j+1],dp[1][x][jj]+dp[0][y][j-jj]);
dp[0][x][j+2]=max(dp[0][x][j+2],dp[0][x][jj]+dp[1][y][j-jj]);
dp[1][x][j+2]=max(dp[1][x][j+2],dp[1][x][jj]+dp[1][y][j-jj]);
}
}
}
}
int main(){
int x,y,i,j;
while(scanf("%d%d",&n,&k)!=EOF){
memset(dp,0,sizeof(dp));
memset(p,0,sizeof(p));
for(i=0;i<210;i++){
old[i].clear();
now[i].clear();
}
for(i=1;i<=n;i++)scanf("%d",val+i);
for(i=1;i<n;i++){
scanf("%d%d",&x,&y);
old[x].push_back(y);
old[y].push_back(x);
}
updata(1);
dfs(1);
printf("%d\n",dp[0][1][k]);
}
return 0;
}
哇啦啦那些该死的状态转移方程真是搞得我想死了!!
不过,幸好我还没死(—>机智的我)于是继续~
下一个是HUD2196 Computer
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<vector>
#define ll long long
using namespace std;
const ll maxn=10010;
struct Edge{
ll v,w;
};
ll f[maxn][2];
bool p[maxn];
vector <Edge> e[maxn];
ll dfs1(ll u){
//printf("dfs : %lld\n",u);
ll i;p[u]=1;
for(i=0;i<e[u].size();i++){
ll v=e[u][i].v;
ll w=e[u][i].w;
if(p[v])continue;
ll tmp=dfs1(v)+w;
//printf("f[%lld][0] : %lld |||| %lld\n",u,f[u][0],tmp);
f[u][0]=max(f[u][0],tmp);
}
return f[u][0];
}
void dfs2(ll u){
p[u]=1;
ll i,max1=0,max2=0,v1,v2;
for(i=0;i<e[u].size();i++){
ll v=e[u][i].v;
ll w=e[u][i].w;
if(p[v])continue;
ll tmp=w+f[v][0];
if(tmp>max1){
max2=max1;v2=v1;
max1=tmp;v1=v;
}else if(tmp==max1||tmp>max2){
max2=tmp;v2=v;
}
}
if(u!=1){
ll tmp=f[u][1],v=-1;
if(tmp>max1){
max2=max1;v2=v1;
max1=tmp;v1=v;
}else if(tmp==max1||tmp>max2){
max2=tmp;v2=v;
}
}
for(i=0;i<e[u].size();i++){
ll v=e[u][i].v;
ll w=e[u][i].w;
if(p[v])continue;
f[v][1]=v==v1?max2+w:max1+w;
dfs2(v);
}
}
int main(){
ll i,j,u,v,k,m,n;
while(scanf("%lld",&n)!=EOF){
for(i=0;i<=n+1;i++)e[i].clear();
memset(f,0,sizeof(f));
for(i=2;i<=n;i++){
scanf("%lld%lld",&u,&v);
e[i].push_back((Edge){u,v});
e[u].push_back((Edge){i,v});
}
/*for(i=1;i<=n;i++)
for(j=0;j<e[i].size();j++)
printf("%lld -> %lld : %lld\n",i,e[i][j].v,e[i][j].w);
*/
memset(p,0,sizeof(p));
dfs1(1);
//for(i=1;i<=n;i++)cout<<f[i][0]<<" ";puts("");
memset(p,0,sizeof(p));
dfs2(1);
for(i=1;i<=n;i++)
printf("%lld\n",max(f[i][0],f[i][1]));
}
return 0;
}
显然,从这些迷之注释就可以看出来本人莫名其妙的错到了多interesting的地步。。。。。于是,坚强的我又挺了过来。
下一个是POJ2342 Anniversary party
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
using namespace std;
const int maxn=6010;
int a[maxn][2],fa[maxn],n,m;
bool p[maxn];
void DP(int x){
p[x]=1;
for(int i=1;i<=n;i++){
if(fa[i]==x&&!p[i]){
DP(i);
a[x][1]+=a[i][0];
a[x][0]+=max(a[i][1],a[i][0]);
}
}
}
int main(){
int x,y,i,j,k;
while(scanf("%d",&n)!=EOF){
memset(a,0,sizeof(a));
memset(p,0,sizeof(p));
memset(fa,0,sizeof(fa));
for(i=1;i<=n;i++)scanf("%d",&a[i][1]);
int rt=0;
while(scanf("%d%d",&x,&y)&&(x||y)){
fa[x]=y;
rt=x==rt?y:rt;
}
while(fa[rt])rt=fa[rt];
DP(rt);
int ans=max(a[rt][1],a[rt][0]);
printf("%d\n",ans);
}
return 0;
}
搞完这个东西就可以说找到一点树形DP的感觉了—–>于是我就走了;
其实不能怪我,是 无聊的LG 又送了我一下状态压缩的东西:)
于是,状态压缩第一枪,POJ3254 Corn Field
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<stack>
using namespace std;
const int maxn=13;
const int maxm=1<<13;
int st[maxm],map[maxm],dp[maxn][maxm];
bool Judge1(int x){return (x&(x<<1))?0:1;};
bool Judge2(int x,int y){return (map[x]&st[y])?0:1;};
void Clear(){
memset(st,0,sizeof(st));
memset(map,0,sizeof(map));
memset(dp,0,sizeof(dp));
}
stack <bool> Stack;
void Out_put(int x){
do{
Stack.push(x&1);
x>>=1;
}while(x>0);
while(!Stack.empty()){
cout<<Stack.top();
Stack.pop();
}
cout<<" ";
}
int main(){
int i,j,k,m,n;
while(scanf("%d%d",&n,&m)!=EOF){
Clear();
for(i=1;i<=n;i++)
for(j=1;j<=m;j++){
int x;
scanf("%d",&x);
if(x==0)map[i]+=1<<(j-1);
}
int top=0;
for(i=0;i<(1<<m);i++)
if(Judge1(i))
st[++top]=i;
for(i=1;i<=top;i++)
if(Judge2(1,i))
dp[1][i]=1;
for(i=2;i<=n;i++)
for(j=1;j<=top;j++){
if(!Judge2(i,j))continue;
for(k=1;k<=top;k++){
if(!Judge2(i-1,k))continue;
if(st[j]&st[k])continue;
dp[i][j]+=dp[i-1][k];
}
}
int ans=0;
for(i=1;i<=top;i++){
ans+=dp[n][i];
ans%=100000000;
}
printf("%d\n",ans);
}
return 0;
}
这道题目做起来还算顺风顺水,具体也没什么好讲的。但是该死的我的。。。。迷之数组都是[0][0][0]…….
下一个没多大区别,
POJ1185 炮兵阵地
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
using namespace std;
int map[105];
int stk[65], sum[65];
int dp[105][65][65],top;
int get(int x){
int ans=0;
while(x){
ans+=x&1;
x>>=1;
}return ans;
}
bool judge(int x){
if(x&(x<<1))return false;
if(x&(x<<2))return false;
return true;
}
void pre(int n){
}
int main(){
int n, m, r, c, i, j, k;
cin>>n>>m;getchar();
for(i=0;i<n;i++,getchar()){
for(j=0;j<m;j++){
char c;
cin>>c;
if(c=='H')map[i]|=1<<j;
}
}
int top=0;
for(int i=0;i<(1<<m);i++)
if(judge(i)){
stk[top]=i;
sum[top++]=get(i);
}
for(i=0;i<top;i++){
if(!(stk[i]&map[0]))
dp[0][0][i]=sum[i];
}
for(i=1;i<n;i++){//line now
for(j=0;j<top;j++){//the situation about this(i) line
if(stk[j]&map[i])continue;
for(k=0;k<top;k++){//the situation about i-1 line
if(stk[k]&stk[j])continue;
for(int tmp=0;tmp<top;tmp++){// the situation about i-2 line
if(stk[j]&stk[tmp])continue;
if(!dp[i-1][tmp][k])continue;
dp[i][k][j]=max(dp[i][k][j],dp[i-1][tmp][k]+sum[j]);
}
}
}
}
int ans=0;
for(i=0;i<top;i++)
for(j=0;j<top;j++)
ans=max(ans,dp[n-1][i][j]);
cout<<ans<<endl;
return 0;
}
然后,开始为猪鱼的到来做准备(虽然不晓得他会不会游过来){
加了个vjudge的比赛;
第一题=洛谷新手村的function=POJ1579;
第二道=POJ1088滑雪;
//显然是两道毫无意义的弱鸡题目,但是被坑了
}
POJ 1579
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#define ll long long
using namespace std;
ll ans[30][30][30];
ll w(ll a,ll b,ll c){
if(a<=0||b<=0||c<=0)return 1;
if(a>20||b>20||c>20)a=b=c=20;
if(ans[a][b][c])return ans[a][b][c];
if(a>20||b>20||c>20)return ans[20][20][20]=1048576;
if(a<b&&b<c)return ans[a][b][c]=w(a,b,c-1)+w(a,b-1,c-1)-w(a,b-1,c);
else return ans[a][b][c]=w(a-1,b,c)+w(a-1,b-1,c)+w(a-1,b,c-1)-w(a-1,b-1,c-1);
}
int main(){
ll i,j,k,m,n,a,b,c;
while(scanf("%lld%lld%lld",&a,&b,&c)!=EOF){
if(a==-1&&b==-1&&c==-1)return 0;
printf("w(%lld, %lld, %lld) = %lld\n",a,b,c,w(a,b,c));
}
return 0;
}
啊,记得long long !!!!!!!!
POJ1088
#include<bits/stdc++.h>
using namespace std;
int a[111][111],ans[111][111],n,m;
int go[4][2]={ {1,0},{-1,0},{0,1},{0,-1} };
bool pd(int x,int y){
if(x<0||y<0||x>n+1||y>m+1)return false;
return true;
}
int dfs(int x,int y){
if(ans[x][y]!=1)return ans[x][y];
//printf("dfs : %d %d\n",x,y);
for(int i=0;i<4;i++){
int u=x+go[i][0];
int v=y+go[i][1];
//if(p[u][v])continue;
if(!pd(u,v))continue;
if(a[x][y]>a[u][v]){
ans[x][y]=max(ans[x][y],dfs(u,v)+1);
}
}
return ans[x][y];
}
int main(){
int i,j,k;
cin>>n>>m;
for(i=1;i<=n;i++)
for(j=1;j<=m;j++){
cin>>a[i][j];
ans[i][j]=1;
}
int ANS=0;
for(i=1;i<=n;i++)
for(j=1;j<=m;j++){
//printf("Start dfs %d %d : \n",i,j);
ANS=max(ANS,dfs(i,j));
//cout<<endl;
}
/*
for(i=1;i<=n;i++,cout<<endl)
for(j=1;j<=m;j++)cout<<ans[i][j]<<" ";
*/
cout<<ANS<<endl;
return 0;
}
显然我是不会告诉你们我忘记边界条件是x<0啊x>n+1,显然我也不会告诉你们我没有把dfs的结果+1更新了。。。。。