[USACO 5.4.2 tour]解题报告

这道题一开始写的费用流,分析时间复杂度的话，差不多O(m)吧，觉得数据这么小很奇怪。

后来看了题解，原来数据小是为DP准备的；想了想DP的话觉得没什么思路，但是看了题解以后却恍然大悟！像这种两个人走的DP，往往是要把两个状态集乘起来。最经典的莫过于传纸条，然后寒假里做了一道高斯消元的题亦然如此。可见这种思想是非常重要的，以后做题的时候一定要有这个意识。

代码（费用流）：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
#include<cstring>
#include<iostream>
#include<string>
string getstr(){
	char c=getchar();
	string ans="";
	while((c<'A'||c>'Z')&&(c<'a'||c>'z')&&(c<'0'||c>'9'))c=getchar();
	for(;(c>='A'&&c<='Z')||(c>='a'&&c<='z')||(c>='0'&&c<='9');c=getchar())ans+=c;
	return ans;
}
#include<map>
struct ES{
	int from,to,cost,c;
}e[20000];
int tot=2,p[1602],next[20000];
int addedge(int from,int to,int cost,int c){
	e[tot]=(ES){from,to,cost,c};
	next[tot]=p[from],p[from]=tot++;
}
int main(){
	freopen("tour.in","r",stdin);
	//freopen("tour.out","w",stdout);
	map<string,int> Map;
	int n,i,m;
	scanf("%d%d",&n,&m);
	if(n==1){
		printf("1\n");
		return 0;
	}
	for(i=1;i<=n;++i)Map[getstr()]=i;
	int u,v;
	for(i=0;i<m;++i){
		u=Map[getstr()],v=Map[getstr()];
		if(u==v)continue;
		if(u>v)swap(u,v);
		u<<=1,v=(v<<1)-1;
		addedge(u,v,0,0x7fffffff),addedge(v,u,0,0);
	}
	for(i=2;i<n;++i)addedge((i<<1)-1,i<<1,1,1),addedge(i<<1,(i<<1)-1,-1,0);
	addedge(1,2,0,2),addedge(2,1,0,0);
	addedge((n<<1)-1,n<<1,0,2),addedge(n<<1,(n<<1)-1,0,0);
	
	int ans=0,dis[1602],q[2001],h,t,path[1602];//最长路抵达此点的边。 
	bool flag[1602]={0};
	for(int k=0;k<2;++k){
		//SPFA
		memset(dis,0xd0,sizeof(dis));
		q[0]=1,h=0,t=1,flag[1]=1,dis[1]=0;
		for(;h!=t;){
			flag[q[h]]=0,i=p[q[h]];
			h=h==2000?0:h+1;
			for(;i;i=next[i])
				if(e[i].c&&dis[e[i].to]<dis[e[i].from]+e[i].cost){
					dis[e[i].to]=dis[e[i].from]+e[i].cost;
					path[e[i].to]=i;
					if(!flag[e[i].to]){
						if(h==t||dis[e[i].to]>dis[q[h]]){
							h=h?h-1:2000;
							q[h]=e[i].to;
						}
						else{
							q[t]=e[i].to;
							t=t==2000?0:t+1;
						}
						flag[e[i].to]=1;
					}
				}
		}
		if(dis[n<<1]<-1E7){
			printf("1\n");
			return 0;
		}
		for(i=n<<1;i;i=e[path[i]].from){
			--e[path[i]].c,++e[path[i]^1].c;
			ans+=e[path[i]].cost;
		}
	}
	printf("%d\n",ans+2);
}

代码（DP）：

#include<cstdio>
using namespace std;
int f[500][500];
#include<iostream>
#include<string>
string getstr(){
	char c=getchar();
	string ans="";
	while((c<'0'||c>'9')&&(c<'a'||c>'z')&&(c<'A'||c>'Z'))c=getchar();
	for(;c>='0'&&c<='9'||c>='a'&&c<='z'||c>='A'&&c<='Z';c=getchar())ans+=c;
	return ans;
}
int p[500],next[5001],succ[5001];
#include<map>
int main(){
	freopen("tour.in","r",stdin);
	//freopen("tour.out","w",stdout);
	int n,m;
	bool flag[500]={0};
	scanf("%d%d",&n,&m);
	int u,v,i,j;
	map<string,int> Map;
	for(i=0;i<n;++i)Map[getstr()]=i;
	int tot=1;
	for(i=0;i<m;++i){
		u=Map[getstr()],v=Map[getstr()];
		if(u==v)continue;
		if(u>v)swap(u,v);
		next[tot]=p[u],p[u]=tot,succ[tot++]=v;
		if(v==n-1)flag[u]=1;
	}
	f[0][0]=1,f[n-1][n-1]=2;
	int k;
	for(i=0;i<n;++i)
		for(j=0;j<n;++j)
			if(f[i][j]){
				for(k=p[i];k;k=next[k])
					if(succ[k]>j)
						f[succ[k]][j]=max(f[succ[k]][j],f[i][j]+1);
				for(k=p[j];k;k=next[k])
					if(succ[k]>i)
						f[i][succ[k]]=max(f[i][succ[k]],f[i][j]+1);
			}
	for(i=0;i<n;++i)
		if(flag[i])
			f[n-1][n-1]=max(f[n-1][n-1],f[i][n-1]+1);
	printf("%d\n",f[n-1][n-1]-1);
}