JZOJ3736数学题
来源:互联网 发布:淘宝直通车手机端位置 编辑:程序博客网 时间:2024/05/21 20:23
题目
Description
Input
输入有多组测例,每组测例有一行,为4 个整数x1,y1, x2, y2,含义见题目描述。输入文件以EOF 结束。
Output
Sample Input
3 0 1 2
6 0 4 0
Sample Output
5
0
Data Constraint
题解
首先我们分析在什么情况下我们容易得解
观察题目,发现题目是要求二元二次函数的极值
设两个向量的模长为p,q;p>q
移项,处理方程后我们可以发现在cosa>1/2(我们设其为结论1)时最优解是p=1,q=0
那么我们现在就要不断的扩大这两个向量的夹角
怎么样才可以扩大这两个向量的夹角呢?
类似于欧几里得算法的思想,我们有
由于前面结论1,我们可以用类似exgcd的方法进行倒推
然后在倒推的时候我们可以先画一条垂线,然后根据±1的长短选择一种加的方法
如果在某一时刻两个向量的夹角大于九十度那么我们就把其中一个向量反向
贴代码
#include<iostream>#include<algorithm>#include<cstdio>#include<cstring>#include<cmath>#define ll long longusing namespace std;ll a,b,c,d,x,y,x11,y11,x2,y2;ll sqr(ll x) { return x*x;}void make(ll x1,ll y1,ll x2,ll y2,ll &x,ll &y){ ll t1; t1=x1*x2+y1*y2; if (t1<0){ make(x1,y1,-x2,-y2,x,y); y=-y; return; } ll t2,t3; t2=x1*x1+y1*y1; t3=y2*y2+x2*x2; if (t2>t3){ make(x2,y2,x1,y1,y,x); return; } if (sqr(t1)*4<=t2*t3){ x=1; y=0; return; } ll k; k=t1/t2; if (2*t1<=t2*(k+k+1)){ make(x1,y1,x2-x1*k,y2-y1*k,x,y); x-=y*k; } else{ make(x1,y1,x2-x1*(k+1),y2-y1*(k+1),x,y); x-=y*(k+1); }}int main(){ freopen("math.in","r",stdin); freopen("math.out","w", stdout); while (scanf("%lld%lld%lld%lld",&x11,&y11,&x2,&y2)!=EOF){ x=1; y=1; make(x11,y11,x2,y2,x,y); printf("%lld\n",sqr(x11*x+x2*y)+sqr(y11*x+y2*y)); }}
0 0
- [JZOJ3736]数学题
- JZOJ3736数学题
- 【jzoj3736】【数学题】【类欧几里德】
- JZOJ3736【NOI2014模拟7.11】数学题(math)
- 【JZOJ3736】【NOI2014模拟7.11】数学题(math)
- 数学题
- 数学题
- 数学题
- 数学题
- 数学题
- 数学题
- 数学题
- 数学题
- 数学题
- 数学题
- 数学题
- 数学题
- 数学题
- TCP/IP报文头部结构整理
- C语言解决对给定的任一不超过1000的正整数n,简单地数一下,需要多少步(砍几下)才能得到n=1?
- 关于政府网站的设计心得
- 聊聊Linux 五种IO模型
- 网络安全思维导图(NS Mind Map)
- JZOJ3736数学题
- Java线程安全和非线程安全
- QT 传输图片出现数据缺失
- jQuery.map() 函数详解
- CUDA之编程中线程分配的数组在register中还是local memory中?
- 扫描台操作以及集合框架
- Java并发编程-并发工具包java.util.concurrent使用指南
- 进制转换
- jqueryMobile手机开发