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;
}