Description
FGD小朋友——一个闻名遐迩的年轻数学家——有一个小MM,yours。
FGD小朋友非常喜欢他的MM,所以他很乐意帮助他的MM做数学作业。
但是,就像所有科学的容器一样,FGD的大脑拒绝不停地重复思考同样的问题。
不幸的是,yours是一个十分用功的学生,所以她不停地让FGD帮助她检查她的作业。
一个阳光明媚的周末,yours的数学老师布置了非常多的寻找多边形的对称轴的题,足够她做相当长的一段时间了。
在此之前FGD已经决定去海边度过这个难得的假期,不过他还是觉得应该帮助他的MM对付可爱的数学作业。
很快地,他找到了解决方案,最好写一个程序来帮助yours检查她的数学作业。
因为FGD并非一个计算机科学家,所以他找到了他的好朋友你,请你帮助他完成这个任务。
请写一个程序:读入多边形的描述计算出每个多边形的对称轴数将计算的结果输出
Input
输入的第一行包含一个正整数t(1<=t<=10),为多边形的边数。
接下来,为t个多边形的描述,每个描述的第一行为一个正整数n(3<=n<=100000),表示了多边形的点数。
然后在后面n行每行两个整数x和y(?100000000<=x, y<=100000000),依次表示多边形的顶点坐标。
多边形不一定是凸的,但是不自交——任何两条边都只有最多一个公共点——他们的公共端点。
此外,没有两条连续的边平行。
Output
你的程序应该输出正好t行,第k行包含了一个整数nk——表示第k个多边形有多少个对称轴。
Sample Input
12
1 -1
2 -1
2 1
1 1
1 2
-1 2
-1 1
-2 1
-2 -1
-1 -1
-1 -2
1 -2
6
-1 1
-2 0
-1 -1
1 -1
2 0
1 1
Sample Output
2
HINT
题解Here!
这是一道神题!题目要求:求多边形的对称轴条数。
乍一看,一道有思维难度的计算几何——先求出凸包,然后开始旋转卡壳。
但是这题比较好玩,不用这么麻烦:
将每一个边和角(叉积)$hash$一下,然后就会得到一个字串,对这个字串跑一遍最长回文子串,其长度即为对称轴数量。
我当时看到这个做法时第一反应:吼啊!妙啊!神奇啊!
至于最长回文子串怎么求,就丢给了$O(n)$的$Manacher$。。。
注:附加多组询问。
附代码:
#include<iostream> #include<algorithm> #include<cstdio> #include<cmath> #include<cstring> #define MAXN 100010 using namespace std; int n,m,p[MAXN<<3],s[MAXN<<3]; struct Point{ int x,y; }a[MAXN]; inline int read(){ int date=0,w=1;char c=0; while(c<'0'||c>'9'){if(c=='-')w=-1;c=getchar();} while(c>='0'&&c<='9'){date=date*10+c-'0';c=getchar();} return date*w; } inline int dis(int i,int j){ return (a[i].x-a[j].x)*(a[i].x-a[j].x)+(a[i].y-a[j].y)*(a[i].y-a[j].y); } inline int cross(int i,int j,int k){ return (a[i].x-a[k].x)*(a[j].y-a[k].y)-(a[j].x-a[k].x)*(a[i].y-a[k].y); } void buildstring(){ for(int i=0;i<n;i++){ s[i<<1]=cross((i-1+n)%n,(i+1)%n,i); s[i<<1|1]=dis(i,(i+1)%n); } m=n<<1; for(int i=0;i<n;i++)s[i+m]=s[i]; m<<=1; } void manacher(){ int ans=0,maxn=0,id; for(int i=0;i<m;i++){ if(maxn>i)p[i]=min(p[id*2-i],p[id]+id-i); else p[i]=1; for(;i-p[i]>=0&&i+p[i]<=m&&s[i+p[i]]==s[i-p[i]];p[i]++); if(p[i]+i>maxn){ maxn=p[i]+i; id=i; } if(p[i]>=n+1)ans++; } printf("%d ",ans); } void work(){ buildstring(); manacher(); } void init(){ n=read(); memset(a,0,sizeof(a)); memset(p,0,sizeof(p)); memset(s,0,sizeof(s)); for(int i=0;i<n;i++){a[i].x=read();a[i].y=read();} } int main(){ int t=read(); while(t--){ init(); work(); } return 0; }