hdu 3902 Swordsman

http://acm.hdu.edu.cn/showproblem.php?pid=3902

题意:给出一个简单多边形,判断这个多边形是否是轴对称多边形.

把原来的n个点扩展,加上各边的中点,扩展为2n个点。

然后对称轴一定是一个点和它对面的一个点组成。(i,i+n)

再验证在这条对称边的两端的对应两点组成的线段是否垂直平分这条对称轴.

理论上的复杂度是 O(n*n),但实际上应该达不到的吧=!

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <string.h>
#include <vector>
using namespace std;

const int maxn = 20005;
const double eps = 1e-8;
struct node
{
	double x,y;
}a[2*maxn];
int n;

bool Can(int s,int e,int st,int et)
{
     if( fabs((a[s].x-a[e].x)*(a[st].x-a[et].x)+(a[s].y-a[e].y)*(a[st].y-a[et].y))>eps)
		 return false;
	 double x=(a[st].x+a[et].x)/2;
	 double y=(a[st].y+a[et].y)/2;
	 if( fabs((x-a[s].x)*(a[e].y-a[s].y)- (a[e].x-a[s].x )*(y-a[s].y) )>eps)
		 return false;
	 return true;
}

bool Is_true(int s,int e)
{
    for(int i=1;i<=n-1;i++)
	{
		int st = s-i;
		if(st<=0)
			st=st+2*n;
		int et = s+i;

		if( Can(s,e,st,et)==false)
			return false;
	}
	return true;
}

int main()
{
   while(scanf("%d",&n)!=EOF)
   {
      for(int i=1;i<=2*n;i+=2)
	      scanf("%lf %lf",&a[i].x,&a[i].y);
	  a[2*n+1]=a[1];
      for(int i=2;i<=2*n;i+=2)
	  {
		  a[i].x=(a[i-1].x+a[i+1].x)/2;
		  a[i].y=(a[i-1].y+a[i+1].y)/2;
	  }
	  int ok=0;
	  for(int i=1;i<=n;i++)
	  {
          if(Is_true(i,i+n)==true)
		  {
             ok=1;
			 break;
		  }	  
	  }
	  if(ok==1)
		  printf("YES\n");
	  else
		  printf("NO\n");
   }
   return 0;
}


下面是另一个版本的代码,复杂度是O(nlogn)

http://hi.baidu.com/nplusnplusnplu/blog/item/d260baef2e9e9c5879f055cb.html

 转化为字符串处理的。

 修改后的版本=

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <string.h>
#include <cmath>
#include <vector>
using namespace std;

const int maxn=20010;
int N,M,next[maxn];
__int64 x[maxn],y[maxn];
__int64 ans;
struct Point 
{
	__int64 ed,ag;
}P[2*maxn],Q[maxn];

bool operator == (Point A,Point B)
{
	return A.ed==B.ed && A.ag==B.ag;
}

void Findnext ()
{
	next[0]=-1;
	next[1]=0;
	for (int i=2 ; i<N ; ++i)
	{
		int p=next[i-1];
		for (; !(Q[p]==Q[i-1])&&p;)
            p=next[p];
		if(Q[p]==Q[i-1])
			next[i]=p+1;
		else 
			next [i]=0;
	}
}

bool comp(__int64 a,__int64 b,__int64 c,__int64 d)
{
	return a*d-b*c>0;
}

inline __int64 sqr(__int64 x)
{
	return x*x;
}
void Getinf()
{
	for (int i=0 ; i<N ; ++i)//正序 
	{
		P[i].ed=sqr(x[i+1]-x[i])+sqr(y[i+1]-y[i]);
		P[i].ag=sqr(x[i+2]-x[i])+sqr(y[i+2]-y[i]);
		if(comp(x[i+2]-x[i],y[i+2]-y[i],x[i+1]-x[i],y[i+1]-y[i]))
			P[i].ag=-P[i].ag;
	}
	for (int i=2 ; i<=N ; ++i )//反序 
	{
		Q[N-i].ed=sqr(x[i]-x[i-1])+sqr(y[i]-y[i-1]);
		Q[N-i].ag=sqr(x[i]-x[i-2])+sqr(y[i]-y[i-2]);
		if(!comp(x[i-2]-x[i],y[i-2]-y[i],x[i-1]-x[i],y[i-1]-y[i]))
			Q[N-i].ag=-Q[N-i].ag; 
	}
	Q[N-1].ed=sqr(x[1]-x[0])+sqr(y[1]-y[0]);
	Q[N-1].ag=sqr(x[1]-x[N-1])+sqr(y[1]-y[N-1]);
	if(!comp(x[N-1]-x[1],y[N-1]-y[1],x[0]-x[1],y[0]-y[1]))
		Q[N-1].ag=-Q[N-1].ag;
	memcpy(P+N,P,N*sizeof(Point));

}

int main ()
{
	while (scanf("%d",&N)!=EOF)
	{
		for (int i=0 ; i<N ; i++)
			scanf("%I64d %I64d",&x[i],&y[i]);
		
		x[N]=x[0];
		y[N]=y[0];
		x[N+1]=x[1];
		y[N+1]=y[1];
		Getinf();

		int p=0;
		ans=0;
		Findnext();
		for(int i=0;i<(N<<1)-1;++i)
		{
			for(;!(P[i]==Q[p])&&p;)p=next[p];
			if(P[i]==Q[p])
			{
				if(p==N-1)
				{
					ans++;//ans记录对称轴的个数 
					p=next[p];
					if(P[i]==Q[p])p++;
					continue;
				}
				p++;
			}
		}
		if(ans)printf("YES\n");
		else printf("NO\n");
	}
	return 0;
}

转载于:https://www.cnblogs.com/lwbaptx/archive/2011/08/04/2128008.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值