计算几何专项：UVa 10556

来源：互联网发布：人工智能专业就业方向编辑：程序博客网时间：2024/05/22 15:39

刚开始用double类型做，就是先把第一个多边形上一个点位移到另一个多边形的对应点处，然后再旋转和位移。然后果断因为精度问题WA了。。。然后改成整形，推一遍公式就AC了。恩，以后不管什么类型的题，能化整形就尽量化成整形，以避免不必要的精度问题。

#include <iostream>#include <cstdio>#include <cmath>using namespace std;struct point{    int x,y;    point(int x=0,int y=0):x(x),y(y){}};point operator-(point a,point b){return point(a.x-b.x,a.y-b.y);}point operator+(point a,point b){return point(a.x+b.x,a.y+b.y);}int length(point a){return a.x*a.x+a.y*a.y;}int cross(point a,point b){return a.x*b.y-a.y*b.x;}int n;point p1[20],p2[20],p3[20];int main(){    freopen("in.txt","r",stdin);    while(cin>>n)    {        if(!n) break;        for(int i=0;i<n;i++)            cin>>p1[i].x>>p1[i].y;        for(int i=0;i<n;i++)            cin>>p2[i].x>>p2[i].y;        for(int i=0;i<n;i++)            p3[i]=p1[0]+(p2[i]-p2[0]);        int flag=1;        point a0=p3[1]-p3[0],b0=p1[1]-p1[0];        for(int i=0;i<n;i++)        {            point a1=p3[i]-p3[0],b1=p1[i]-p1[0];            if(length(a1)*length(b0)==length(b1)*length(a0)&&               cross(a0,b0)*cross(a0,b0)*length(a1)*length(b1)==cross(a1,b1)*cross(a1,b1)*length(a0)*length(b0))                continue;            flag=0;            break;        }        if(flag) cout<<"similar"<<endl;        else cout<<"dissimilar"<<endl;    }    return 0;}