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

原创 2017年01月03日 16:23:56

题意:给出一对三角形,已知其中一个是由另一个绕某一点旋转放缩而得到,求这个点的坐标。每个测试点有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;
}
版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

BZOJ 3621 我想那还真是令人高兴啊 复数

题目大意:给定两个三角形,其中一个可以通过以某点为中心旋转并放缩的方式得到另一个,求这个中心 それはとっても嬉しいなって。 首先两个复数相乘的几何意义是【极角相加,长度相乘】 这两种变换正好对应旋...

自己对ajax的第一次上手被人说菜死还真是菜

今天算第一次上手ajax,为了再当前页面验证用户名和

poj 2976 && poj 2728 && poj3621 && bzoj 1486 01分数规划

关于01分数规划的具体讲解和例题分析,推荐相当好的博客:http://www.cnblogs.com/perseawe/archive/2012/05/03/01fsgh.html 这里不再赘述 所以...

bzoj1690/poj3621[Usaco2007 Dec]奶牛的旅行

0-1分数规划+Bellman-Ford(最优比率环)

函数重载——复数运算

运算符的重载(复数的相关运算)

运算符的重载实际是一种特殊的函数重载,必须定义一个函数,并告诉C++编译器,当遇到该重载的运算符时调用此函数。这个函数叫做运算符重载函数,通常为类的成员函数。      定义运算符重载函数的一般格式:...

C#复数运算

  • 2017-03-25 16:46
  • 48KB
  • 下载

运算符重载——重载+和-对复数类对象进行加减运算

1.题目: Problem Description 定义一个复数类,该类包含两个double型的数据成员代表复数的实部和虚部,包含构造函数(默认值为0,0),和显示函数,现重载运算符+和-,使其能对...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)