poj 1739 Tony's Tour

最新推荐文章于 2012-10-10 17:08:09 发布
最新推荐文章于 2012-10-10 17:08:09 发布
阅读量1.6k 收藏 1
点赞数
分类专栏： 插头dp 文章标签： up struct 2010
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/xymscau/article/details/8035484
版权

楼教主的男人八题之一，这是我切的第一题，昨天搞定了括号表示法，今天趁热打铁，这题很容易想到把左下和右下两个格子连起来，构成一条哈密顿回路，就和上一题一样，不过这条路有一条必须经过的路径，我们需要特殊处理，做法有两种，我是添加了一层，然后特殊处理最后一层，后来发现小hh的解法说不用，又写了另一种解法，就是将最后一层的状态变成初始状态，然后从底往上dp，这种写法比较好，简单，而且状态比第一种写法少。

Run IDUserProblemResultMemoryTimeLanguageCode LengthSubmit Time
108732392010307204251739Accepted1456K110MSG++2941B2012-10-01 17:28:50
108731932010307204251739Accepted1444K125MSG++3493B2012-10-01 17:15:55
108730992010307204251739Accepted1456K94MSG++2943B2012-10-01 16:48:33
第二种写法 

#include<cstdio>
#include<cstring>
#include<string>
#include<iostream>
#include<algorithm>
#include<cmath>
#define LL long long
using namespace std;
const int maxn=30001;
int n,m,lx,ly,rx,ry,mov[9]={0,2,4,6,8,10,12,14,16};
LL getbit(int st,int k){return (st>>mov[k])&3;}
LL pybit(int st,int k){return (st<<mov[k]);}
LL clrbit(int st,int a,int b){return st&(~(3<<mov[a]))&(~(3<<mov[b]));}
struct node
{
	int head[maxn],next[maxn],size;
	LL sum[maxn],sta[maxn];
	void clear()
	{
		memset(head,-1,sizeof(head));
		memset(sum,0,sizeof(sum));
		size=0;
	}
	void push(LL st,const LL v)
	{
		LL hash=st%maxn;
		for(int i=head[hash];i>=0;i=next[i])
		{
			if(sta[i]==st)
			{
				sum[i]+=v;
				return;
			}
		}
		sum[size]=v,sta[size]=st;
		next[size]=head[hash],head[hash]=size++;
	}
}dp[2];
int fl(LL st,int k)
{
	int cnt=1;
	for(int j=k+1;j<=m;j++)
	{
		int e=getbit(st,j);
		if(e==2)cnt--;
		else if(e==1)cnt++;
		if(cnt==0) return j;
	}
}
int fr(LL st,int k)
{
	int cnt=1;
	for(int j=k-1;j>=0;j--)
	{
		int e=getbit(st,j);
		if(e==2)cnt++;
		else if(e==1)cnt--;
		if(cnt==0) return j;
	}
}
char gp[20][20];
LL DP()
{
	LL ans=0;
	dp[0].clear();dp[0].push(pybit(1,0)|pybit(2,m-1),1);
	int now=0,pre=1;
	for(int i=n;i>=1;i--)
	{
		pre=now,now^=1,dp[now].clear();
		for(int j=0;j<dp[pre].size;j++)
			dp[now].push(dp[pre].sta[j]<<2,dp[pre].sum[j]);
		for(int j=1,jj=m;j<=m;j++,jj--)
		{
			pre=now,now^=1,dp[now].clear();
			for(int k=0;k<dp[pre].size;k++)
			{
				LL l=getbit(dp[pre].sta[k],j-1);
				LL up=getbit(dp[pre].sta[k],j);
				LL st=clrbit(dp[pre].sta[k],j-1,j);
			   // cout<<i<<' '<<jj<<' '<<l<<' '<<up<<' '<<st<<' '<<dp[pre].sta[k]<<' '<<dp[pre].sum[k]<<endl;
				if(!l&&!up)
				{
					if(gp[i][jj]=='#')
						dp[now].push(dp[pre].sta[k],dp[pre].sum[k]);
					else if(i>1&&jj>1&&gp[i-1][jj]=='.'&&gp[i][jj-1]=='.')
						dp[now].push(st|pybit(1,j-1)|pybit(2,j),dp[pre].sum[k]);
				}
				else if(!l||!up)
				{
					LL e=l==0?up:l;
					if(i>1&&gp[i-1][jj]=='.')
					dp[now].push(st|pybit(e,j-1),dp[pre].sum[k]);
					if(jj>1&&gp[i][jj-1]=='.')
					dp[now].push(st|pybit(e,j),dp[pre].sum[k]);
				}
				else if(l==1&&up==1)
					dp[now].push(st^pybit(3,fl(st,j)),dp[pre].sum[k]);
				else if(l==2&&up==2)
					dp[now].push(st^pybit(3,fr(st,j-1)),dp[pre].sum[k]);
				else if(l==2&&up==1)
					dp[now].push(st,dp[pre].sum[k]);
				else if(lx==i&&ly==jj)
					dp[now].push(st,dp[pre].sum[k]);
			}
		}
	}
	for(int i=0;i<dp[now].size;i++)
		if(dp[now].sta[i]==0)
			ans+=dp[now].sum[i];
	return ans;
}
int main()
{
	while(~scanf("%d%d",&n,&m))
	{
		if(n==0&&m==0) break;
		for(int i=1;i<=n;i++)
			scanf("%s",&gp[i][1]);
		lx=ly=0;
		for(int i=1;i<=n&&!lx;i++)
		{
			for(int j=1;j<=m&&!lx;j++)
				if(gp[i][j]=='.')
					lx=i,ly=j;
		}
	//	cout<<lx<<' '<<ly<<endl;
		cout<<DP()<<endl;
	}
	return 0;
}

第一种写法 

#include<cstdio>
#include<cstring>
#include<string>
#include<iostream>
#include<algorithm>
#include<cmath>
#define LL long long
using namespace std;
const int maxn=30001;
int n,m,lx,ly,rx,ry,mov[9]={0,2,4,6,8,10,12,14,16};
LL getbit(int st,int k){return (st>>mov[k])&3;}
LL pybit(int st,int k){return (st<<mov[k]);}
LL clrbit(int st,int a,int b){return st&(~(3<<mov[a]))&(~(3<<mov[b]));}
struct node
{
	int head[maxn],next[maxn],size;
	LL sum[maxn],sta[maxn];
	void clear()
	{
		memset(head,-1,sizeof(head));
		memset(sum,0,sizeof(sum));
		size=0;
	}
	void push(LL st,const LL v)
	{
		LL hash=st%maxn;
		for(int i=head[hash];i>=0;i=next[i])
		{
			if(sta[i]==st)
			{
				sum[i]+=v;
				return;
			}
		}
		sum[size]=v,sta[size]=st;
		next[size]=head[hash],head[hash]=size++;
	}
}dp[2];
int fl(LL st,int k)
{
	int cnt=1;
	for(int j=k+1;j<=m;j++)
	{
		int e=getbit(st,j);
		if(e==2)cnt--;
		else if(e==1)cnt++;
		if(cnt==0) return j;
	}
}
int fr(LL st,int k)
{
	int cnt=1;
	for(int j=k-1;j>=0;j--)
	{
		int e=getbit(st,j);
		if(e==2)cnt++;
		else if(e==1)cnt--;
		if(cnt==0) return j;
	}
}
char gp[20][20];
LL DP()
{
	LL ans=0;
	dp[0].clear();dp[0].push(0,1);
	int now=0,pre=1;
	for(int i=1;i<=n;i++)
	{
		pre=now,now^=1,dp[now].clear();
		for(int j=0;j<dp[pre].size;j++)
			dp[now].push(dp[pre].sta[j]<<2,dp[pre].sum[j]);
		for(int j=1;j<=m;j++)
		{
			pre=now,now^=1,dp[now].clear();
			for(int k=0;k<dp[pre].size;k++)
			{
				LL l=getbit(dp[pre].sta[k],j-1);
				LL up=getbit(dp[pre].sta[k],j);
				LL st=clrbit(dp[pre].sta[k],j-1,j);
		//	    cout<<i<<' '<<j<<' '<<l<<' '<<up<<' '<<st<<' '<<dp[pre].sta[k]<<' '<<dp[pre].sum[k]<<endl;
				if(!l&&!up)
				{
					if(gp[i][j]=='#')
						dp[now].push(dp[pre].sta[k],dp[pre].sum[k]);
					else if((i<n&&j<m&&gp[i+1][j]=='.'&&gp[i][j+1]=='.')||i==n)
						dp[now].push(st|pybit(1,j-1)|pybit(2,j),dp[pre].sum[k]);
				}
				else if(!l||!up)
				{
					LL e=l==0?up:l;
					if((i<n&&gp[i+1][j]=='.')||i==n)
					dp[now].push(st|pybit(e,j-1),dp[pre].sum[k]);
					if(j<m&&gp[i][j+1]=='.')
					dp[now].push(st|pybit(e,j),dp[pre].sum[k]);
				}
				else if(l==1&&up==1)
					dp[now].push(st^pybit(3,fl(st,j)),dp[pre].sum[k]);
				else if(l==2&&up==2)
					dp[now].push(st^pybit(3,fr(st,j-1)),dp[pre].sum[k]);
				else if(l==2&&up==1)
					dp[now].push(st,dp[pre].sum[k]);
			}
		}
	}
	pre=now,now^=1,dp[now].clear();
	for(int j=0;j<dp[pre].size;j++)
		dp[now].push(dp[pre].sta[j]<<2,dp[pre].sum[j]);
	for(int j=1;j<=m;j++)
	{
		pre=now,now^=1,dp[now].clear();
		for(int k=0;k<dp[pre].size;k++)
		{
			LL l=getbit(dp[pre].sta[k],j-1);
			LL up=getbit(dp[pre].sta[k],j);
			LL st=clrbit(dp[pre].sta[k],j-1,j);
		//	cout<<j<<' '<<l<<' '<<up<<' '<<st<<' '<<dp[pre].sta[k]<<' '<<dp[pre].sum[k]<<endl;
			if(j==1&&up==1)
				dp[now].push(st|pybit(1,j),dp[pre].sum[k]);
			if(j!=m&&l==1&&up==0)
				dp[now].push(st|pybit(1,j),dp[pre].sum[k]);
			if(l==1&&up==2&&j==m)
				dp[now].push(st,dp[pre].sum[k]);
		}
	}
	for(int i=0;i<dp[now].size;i++)
		if(dp[now].sta[i]==0)
			ans+=dp[now].sum[i];
	return ans;
}
int main()
{
	while(~scanf("%d%d",&n,&m))
	{
		if(n==0&&m==0) break;
		for(int i=1;i<=n;i++)
			scanf("%s",&gp[i][1]);
		lx=ly=0;
		for(int i=1;i<=n&&!lx;i++)
		{
			for(int j=1;j<=m&&!lx;j++)
				if(gp[i][j]=='.')
					lx=i,ly=j;
		}
	//	cout<<lx<<' '<<ly<<endl;
		cout<<DP()<<endl;
	}
	return 0;
}

xymscau
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
poj.rar_MáS
09-21
题意：给一个无向图，求从起点s到终点e尽力n条边的路径最小值。 边的数量2<=t<=100(每条边可以重复遍历) 算法：利用dp求出 path[ l ] [ i ] [ k ]=MIN(path[ l ][ i ] [ k ] , path[ l ][ i ] [ j ] + path...
poj.rar_poj
09-23
标题中的"poj.rar_poj"暗示了这是一个与POJ（Programming Online Judge）相关的压缩文件，POJ是一个在线编程竞赛平台，主要供程序员们练习和提交算法解决方案。在这个压缩包中，很可能包含了用户在POJ上参与编程挑战...
POJ.rar_poj
09-19
【标题】"POJ.rar_poj" 指的是一份与编程竞赛相关的压缩文件，主要涉及POJ（Problem Overloading Judge）平台上的问题解答。POJ是一个在线的编程竞赛平台，它为参赛者提供了大量的算法题目进行练习和比赛，以提升...
POJ2262-Goldbach's Conjecture
07-29
【标题】"POJ2262-Goldbach's Conjecture"是一个编程竞赛题目，源自北京大学的在线判题系统POJ（Problem Online Judge）。这个题目是基于著名的数学猜想——哥德巴赫猜想（Goldbach's Conjecture）设计的。 【描述...
hdu 1693 插头dp
XYM-AC之路
09-07 6068
题意：在n*m的矩阵中，有些格子有树，没有树的格子不能到达，找一条或多条回路，吃完所有的树，求有多少中方法。 这题是插头dp，刚刚学习，不是很熟悉，研究了好几天才明白插头dp的方法，他们老是讲一些什么轮廓线啊，插头啊什么的，刚开始完全不知道这些事干什么的，看完cdq的大作后也是
插头dp的几个模板
XYM-AC之路
10-10 4537
/* ural1519 求经过所有可行点的哈密顿回路的个数 括号匹配法,转移有点复杂,但是时间空间比较小 */ #include #include #include #include #include #include #include #include #define LL long long using namespace std; const int maxn=30001; int n,m,
uva10572 black&white 插头dp
XYM-AC之路
10-10 1499
终于切了第一道广义路径的题目，也确实感到比较复杂，最小表示法的优势也可以体现出来，因为这题，昨天的澡到现在才洗完，按照小hh的思路，用一条轮廓线记录联通状态，注意这里的轮廓线已经不是m+1啦，因为我们不需要左插头了，每个格子可以插至多4个插头，另外一条轮廓线表示当前的染色状态，这个要m+1因为还要保留左上格子的颜色，然后我们对于4个格子的颜色讨论，左，上，左上，当前，一共16中情况，8种是对称的，
poj 3133 Manhattan Wiring 插头dp
XYM-AC之路
10-08 1497
这题的和上一份报告那题有点像，不过是确定了起点和终点，而且分别有两对，按照小hh的说法就是确定了起点和终点就相当于回路问题，另一边相当于虫洞连接起来，不过他说起点为左括号，终点为右括号，这还是有点问题的，昨晚思考了一下，如果，起点和终点的位置在其他位置，而不是边界，想想楼教男人8题的那题，如果位置在中间该怎么处理，左括号和右括号的设置方法肯定不可行，而是应该像上篇报告那题要独立插头，但是这个独立插
zoj 3213 Beautiful Meadow 插头dp
XYM-AC之路
10-07 1058
求任意路径可以得到的最大权值和，这题和前面的题目不一样，插头的起点和结束点都不是确定的，所以要加一个插头表示单插头，表示起点或终点，接着就是恶心的转移，一直漏一个条件，后来看小hh的代码发现问题，终于过了。 Run ID Submit Time Judge Status Problem ID Language Run Time(ms) Run Memory(KB)
zoj 3466 The Hive II
XYM-AC之路
09-28 922
小hh插头dp http://www.notonlysuccess.com/index.php/plug-dp-complete/ 的第二题。 这题比第一题难处理一点，但是方法是一样的，用轮廓线扫一遍就好了，每个可行格子有且仅有两个插头，要分奇偶格子处理 两条轮廓线就像这样，然后层层扫描即可 Run ID Submit Time Judge Status
ural 1519. Formula 1
XYM-AC之路
10-01 890
插头dp第二题，debug了一个晚上，导致今天中秋都没有看到月亮，不过总算搞定了，思维习惯上打错个东西，找的好久，思路其实和前面轮廓线扫描差不多，但是要保证在最后一个合并最后两个插头，这样才能保证只有一条回路，单纯的二进制扫描法不能完成这个任务，所以要用括号表示法，cqd的ppt说的很明白，我就不说明啦，要用三进制表示，我们用二进制的两位来表示三进制，这样状态有2^(26)，太多啦，就算空间允许，
zoj 3256 Tour in the Castle
XYM-AC之路
10-03 863
这题是poj1739的改编版，在那份题解里我用了两种写法，这题的m很大，显然快速幂，BFS所有状态，不过矩阵也相当的大，我算了最大状态是114这么多,所以我优化了快速幂，预存幂矩阵，然后n^2快速幂，跑了600+,那个第一的200+ms是有减少状态么，还是有其他优化方法，我被他时间空间完胜了，Orz。 Run ID Submit Time Judge Status Prob
hdu 3377 Plan
XYM-AC之路
10-04 851
昨晚熬夜把它给切了，1y，这题可以相当于poj 1739+hdu 1964的加强版，我的做法是在最左边和最下边加上一层外围，权值全为0,然后就是插头dp了，不过要注意的是，题目不要求全部格子都经过，所以有一边插头的格子一定要有另一边，而没有插头的格子可以选择放两个插头，或者不放，另外要再注意下加的那两层的插头状态是确定的。 Run ID Submit Time Judge S
fzu 1977 Pandora adventure
XYM-AC之路
10-05 841
这题比上题更简单，没有确定的起始点，是求回路，不过有些点是必须走的，那么这些点上必须有插头，否则不能转移，最后两个插头合并只能在所有的必须点扫描后。 RunID Submit Time PID Language Time Mem Len 466516 2012-10-05 12:52:39 1977 GNU
hdu 4285 circuits
XYM-AC之路
10-05 827
被这题卡了一整天，测试无数的数据，和对拍数据无误，几乎崩溃，后来重新敲时，发现自身写法的漏洞，果不其然，就是这个潜伏多年的bug，在ural，poj，以及我的vim下得出的数据都是一样的，我现在有些明白，为什么以前出现过在G++AC的题目，到C++是wa，应该是写法上的漏洞，编译器不同结果是不同的，改完这个bug，就AC了，rank排到第一去了。 Run ID Submit T
NET语言程序设计课件-第2章 Visual C#.NET语法基础.ppt
最新发布
08-08
Visual C# .NET是一门面向对象的程序设计语言，是当前流行的.NET系列的语言。 本课程以Visual Studio为程序设计环境，对C#进行了全面阐述。 本电子课件适用于课程讲解或者课程分析参考使用。
[毕业设计]PHP图书借阅与归还管理系统(源代码+论文).zip
08-08
[毕业设计]PHP图书借阅与归还管理系统(源代码+论文)
ssm_012_mysql宝康药房销售管理系统.zip
08-08
随着我国市场经济的蓬勃发展和人们对医药产品需求的迅速增加，医药销售行业正处于一个高速发展的时期。行业的快速发展必然导致竞争的加剧，面对药品销售业日益严酷的竟争现实，加强管理、提高工作效率和改善服务质量成了急待解决的问题。而解决这些问题的关键措施之一就是利用计算机等现代信息技术，建立实用、先进、高效的药房销售管理系统，引进创新的经营机制，实现药品销售过程的全面信息化管理，以适应企业生存和发展的新形势。因此，开发一个宝康药房销售管理系统是十分必要的。 根据需求，确定系统采用JSP技术，SSM框架，JAVA作为编程语言，MySQL作为数据库。整个系统要操作方便、易于维护、灵活实用。主要实现了频道维护、新闻维护、药品维护、订单维护、系统管理等功能。
POJ2135 Farm Tour —— 最小费用最大流
05-25
题目描述 农夫约翰在他的农场里建立了一些水沟，以帮助水从一个地方流到另一个地方。他需要为这些水沟安排一些管道，以便水可以通过它们流动。然而，由于约翰的农场是非常大的，因此在不同位置的水沟之间建立管道的...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值