BZOJ 1100 POI2007 對稱軸osi 計算幾何+KMP演算法
阿新 • • 發佈:2019-01-02
題目大意:給定一個多邊形,求對稱軸數量
我X 這究竟是怎麼想到KMP的……
首先 將邊字元化 即找到這個多邊形的中心 然後用與中心構成的三角形的邊-角-邊的方式表示這條邊
將邊順時針掃一遍 然後倍增至長度為2n-1 再逆時針掃一遍 逆時針掃的那遍在順時針那遍中出現的次數就是對稱軸數目
用KMP演算法就能搞出來 證明自己YY吧
出題人卡精度喪心病狂。。。
#include <cmath> #include <cstdio> #include <cstring> #include <iostream> #include <algorithm> #define M 100100 #define EPS 1e-6 using namespace std; struct point{ double x,y; point(){} point(double _,double __): x(_),y(__){} void Read() { scanf("%lf%lf",&x,&y); } void operator += (const point &Y) { x+=Y.x;y+=Y.y; } point operator - (const point &Y) const { return point(x-Y.x,y-Y.y); } double operator * (const point &Y) const { return x*Y.y-Y.x*y; } point operator / (double a) const { return point(x/a,y/a); } }points[M],centre; struct line{ double d1,d2,cross; line(){} line(double _,double __,double ___): d1(_),d2(__),cross(___/_/__){} bool operator == (const line &Y) const; }a[M],b[M<<1]; int n; bool line :: operator == (const line &Y) const { if(fabs(d1-Y.d1)>EPS) return false; if(fabs(d2-Y.d2)>EPS) return false; if(fabs(cross-Y.cross)>EPS) return false; return true; } double Distance(const point &p1,const point &p2) { return sqrt( (p1.x-p2.x)*(p1.x-p2.x) + (p1.y-p2.y)*(p1.y-p2.y) ); } int KMP(int len) { static int next[M]; int i,fix=0,re=0; for(i=2;i<=n;i++) { while( fix && !(a[fix+1]==a[i]) ) fix=next[fix]; if( a[fix+1]==a[i] ) ++fix; next[i]=fix; } fix=0; for(i=1;i<=len;i++) { while( fix && !(a[fix+1]==b[i]) ) fix=next[fix]; if( a[fix+1]==b[i] ) ++fix; if(fix==n) ++re,fix=next[fix]; } return re; } int main() { //freopen("osi.in","r",stdin); //freopen("osi.out","w",stdout); int T; for(cin>>T;T;T--) { int i; cin>>n;centre=point(0.0,0.0); for(i=1;i<=n;i++) points[i].Read(),centre+=points[i]/static_cast<double>(n); for(i=1;i<=n;i++) b[i]=b[i+n]=line(Distance(points[i],centre),Distance(points[i%n+1],centre), (points[i]-centre)*(points[i%n+1]-centre) ); for(i=n;i;i--) a[n-i+1]=line(b[i].d2,b[i].d1,0),a[n-i+1].cross=b[i].cross; cout<<KMP(n+n-1)<<endl; } }