JZOJ3736【NOI2014模拟7.11】数学题(math)

题目就不给了。

题解：干啊！！！！！老子从早思考到晚终于下午3：20切了这道要人命的题目。。
我第一次认识到c++的误差到底是有多烂，sqrt过去sqr回来居然结果不一样！！坑了我半个小时。
首先，我建议你先看完论文再来看我的题解，因为论文中有些地方很明显作者笔误&怠惰。
首先这道题目，你要把题目模型转换好，不然就会像我上午比赛一样啥都没搞出来。
前言：论文中关于当向量a,b同向时结果为gcd(|a|,|b|)明显错误。并不是结果，而是x和y的分母。这里可能有很多人就不懂了（其实应该没人懂。。）
我们知道，向量同向，那么只要反向的话其实就是答案。。那么我们可以通过a/gcd,-b/gcd然后求出lcm，这样不就刚好圆上了？
然后我们来化简模型。。
其实就是论文里的那个式子|ax+by|^2。。
这个一看就容易联想到exgcd，上午那个想了半天gcd根本没有任何触及。。
然后，这跟gcd有啥关系呢。。
看了论文都知道，，有个结论，就是a,b之间夹角>3/pi的时候，答案为min(|a|,|b|)。
那么其实接下来就很明显了，就是如何变换使得(a,b)之间夹角>3/pi，那么又有一个结论出来了。。(a,b)=(a,b+ka)；
那么。。
用二次函数解极值的方法求出k，每次逼近就可以了，cos的话用点积算。
判断是否为0不要用斜率！！c++误差大的惊人，所以用叉积就好了。。
算向量的模长的时候不要直接开方，用的时候在开方，不然最后你统计答案用的是平方，sqrt之后再sqr回去答案会有+-1的误差！
code(短小精悍)

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<math.h>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define sqr(x) ((x)*(x))
#define eps 0.5
using namespace std;
int n,m;
const int N=1e5+5;
typedef long long ll;
struct node
{
    ll x,y;
}a,b;
inline int dianji(ll x1,ll y1,ll x2,ll y2)
{
    return x1*x2+y1*y2;
}
inline ll dist(ll x1,ll y1)
{
    return sqr(x1)+sqr(y1);
}
inline int chaji(ll x1,ll y1,ll x2,ll y2)
{
    return x1*y2-x2*y1;
}
int main()
{
    freopen("math.in","r",stdin);
    freopen("math.out","w",stdout);
    while (scanf("%lld%lld%lld%lld",&a.x,&a.y,&b.x,&b.y)!=EOF)
    {
        if (!chaji(a.x,a.y,b.x,b.y))printf("0\n");
        else
        {
            while (1)
            {
                if (dist(a.x,a.y)>dist(b.x,b.y))
                {
                    node t=a;
                    a=b;
                    b=t;
                }
                double cos=dianji(a.x,a.y,b.x,b.y)/sqrt(dist(a.x,a.y)*dist(b.x,b.y));
                if (cos<0)a.x=-a.x,a.y=-a.y;
                else
                {
                    if (cos<eps)break;
                    int k=sqrt(dist(b.x,b.y))*cos/sqrt(dist(a.x,a.y))+eps;
                    b.x-=k*a.x,b.y-=k*a.y;
                }
            }
            ll ans=min(dist(a.x,a.y),dist(b.x,b.y));
            printf("%lld\n",ans);
        }
    }
    return 0;
}