リーインカーネイション

気づけば I came too far 此博客无限期停用

BZOJ3621: 我想那还真是令人高兴啊 复数运算

题意:给出一对三角形,已知其中一个是由另一个绕某一点旋转放缩而得到,求这个点的坐标。每个测试点有T(<=10)组数据。
一开始以为是计算几何,然而向量的旋转和放缩是两种操作,复数却可以一步完成。
对于复数a和b,a*b的几何意义为a与b长度相乘,极角相加。
那么本题就可以看成从旋转点出发,指向其中一个三角形顶点的三个复数向量乘以同一个复数后变成了从旋转点出发,指向另一个三角形顶点的三个复数向量。设出旋转点坐标和同乘的复数,枚举顶点的六种可能对应关系。列出三个方程,用前两个方程解出解,代到第三个方程中检验是否正确即可。
解方程需要用到复数除法。对于(a+bi)/(c+di),分子分母同乘(c-di)可进行分母有理化,
整理得结果为((ac+bd)/(cc+dd))+((bc-ad)/(cc+dd))i。
然而复数除法非常掉精度,取eps=1e-5都会WA,要1e-4才能过。。。

#include<cstdio>
#include<cmath>
using namespace std;
#define y1 __y__
struct cpx
{
    double a,b;
    cpx(const double &a=0,const double &b=0):a(a),b(b){}
    cpx(FILE* f){fscanf(f,"%lf%lf",&a,&b);}
    #define c ano.a
    #define d ano.b
    inline cpx operator + (const cpx &ano) const
    {return cpx(a+c,b+d);}
    inline cpx operator - (const cpx &ano) const
    {return cpx(a-c,b-d);}
    inline cpx operator * (const cpx &ano) const
    {return cpx(a*c-b*d,a*d+b*c);}
    inline cpx operator / (const cpx &ano) const
    {
        double temp=c*c+d*d;
        return cpx((a*c+b*d)/temp,(b*c-a*d)/temp);
    }
    inline bool operator == (const cpx &ano) const
    {return fabs(a-c)<1e-4&&fabs(b-d)<1e-4;}
    inline void operator() (FILE* f){fprintf(f,"%lf %lf\n",a,b);}
}ans,a1,a2,a3;
int t;
inline bool check(const cpx &b1,const cpx &b2,const cpx &b3)
{
    cpx __t=(b1-b2)/(a1-a2);
    ans=(a1*__t-b1)/(__t-cpx(1.0,0.0));
    return (a3-ans)*__t==(b3-ans);
}
int main()
{
    scanf("%d",&t);
    while(t--)
    {
        a1=stdin,a2=stdin,a3=stdin;
        cpx b1(stdin),b2(stdin),b3(stdin);
        if(check(b1,b2,b3)) goto successed_to_match;
        if(check(b1,b3,b2)) goto successed_to_match;
        if(check(b2,b1,b3)) goto successed_to_match;
        if(check(b2,b3,b1)) goto successed_to_match;
        if(check(b3,b1,b2)) goto successed_to_match;
        if(check(b3,b2,b1)) goto successed_to_match;
        successed_to_match:
        ans(stdout);
    }
    return 0;
}
阅读更多
版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/Mima_Reincarnation/article/details/53995565
所属专栏: BZOJ做题纪录
想对作者说点什么? 我来说一句

android应用添加广告测试

2011年07月18日 82KB 下载

没有更多推荐了,返回首页

不良信息举报

BZOJ3621: 我想那还真是令人高兴啊 复数运算

最多只允许输入30个字

加入CSDN,享受更精准的内容推荐,与500万程序员共同成长!
关闭
关闭