Codeforces Gym
来源:互联网 发布:刚毕业的程序员工资 编辑:程序博客网 时间:2024/05/18 01:33
题目
题意
给你两个数X和Y,要你构造一个三角形或一个四边形(凹凸任意),满足以下条件:
1.至少一个点的x坐标为0
2.至少一个点的x坐标为X
3.至少一个点的y坐标为0
4.至少一个点的y坐标为Y
5.图形面积为<=25000
6.任意点坐标(x,y),x和y必须是整数,且范围:0<=x<=X,0<=y<=Y
数据范围:n个case,1<=n<=1e5,2<=X<=1e9,2<=Y<=1e9
分析
首先考虑下极限数据,X = Y = 1e9
数据很大,面积又规定2.5e4以内,那考虑下凹四边形,一个方案是这样的
四个点坐标分别是O(0,0),B(1e9-1,1e9),A(1,1),C(1e9,1e9-1),这样面积是OA * BC / 2 = 1
这个策略是这样的:B的位置是(X,Y)这个点往左一格,即B(X-1,Y),C的位置是(X,Y)这个点往下一格,即C(X,Y-1)
由于BC段长度固定为1,面积由A点坐标决定,设A(a,b),考虑A坐标为整数时最接近O点的坐标取值
由于A点必须夹在OB和OC中间,由斜率有Y/(X-1) < b/a < (Y-1)/X,考虑b=Y且a=X的情况下不等式恒成立,也就是A点在右上角那里
但这样面积太大,由于分子分母同时除以一个数后,分式值不变,那么要在坐标为整数的情况下a和b取值最小
则令b = Y/g,a = X/g,g为X和Y的最大公约数
接下来计算此时面积S,我们先推导一个小公式,因为待会三角形也要用到,当然你硬着算也行。
推个OA一般式再代个点到直线距离,得此时面积S = |bc - ad|/2 ,因为B点在下方,由斜率,S = (bc-ad)/2
回到上面面积计算,分为OAB和OAC两个三角形,利用推导的公式,B(X-1,Y),A(a,b),C(X,Y-1)
面积S = (Y*a - b*(Y-1))/2 + (X*b - (Y-1)*a)/2 = (a+b)/2 = (X/g + Y/g)/2
问题在于g很小且X和Y很大的时候,面积会超过2.5e4的限制
考虑三角形,就再拿第二张图来用下好了,令A点为(X,Y),即a = X,b = Y,已经推导S = (bc-ad)/2
观察bc-ad,c和d未知,a和b已知,很像exgcd的形式,ax-by,令g = gcd(a,b),解方程bx-ay = g,此时面积S = g/2,这样,在g很小情况下,就可以构造一个满足要求的三角形
为什么题目是一定有解的?
我们看下两个面积,一个是S1 =(X/g + Y/g)/2,一个是S2 =g/2,试着将其相乘再开根,sqrt(S1*S2) = (sqrt(X+Y))/2<=25000
也就是说S1和S2不会同时超过25000,得证
代码
#include<iostream>#include<cstdio>#define test(a) cout<<(a)<<endl#define test2(a,b) cout<<(a)<<" "<<(b)<<endl#define test3(a,b,c) cout<<(a)<<" "<<(b)<<" "<<(c)<<endl#define IOS ios::sync_with_stdio(false);cin.tie(0)using namespace std;void gcd(int a,int b,int &d,int &x,int &y) { if(!b) {d = a;x = 1;y = 0;} else {gcd(b,a%b,d,y,x);y-=x*(a/b);}}int myabs(int a) { if (a>0) return a; return -a;}void work() { int a,b,g,x,y; cin>>a>>b; gcd(a,-b,g,y,x); g = myabs(g); if (g<50000) { test(3); test2(0,0); test2(a,b); test2(x,y); } else { test(4); test2(0,0); test2(a-1,b); test2(a/g,b/g); test2(a,b-1); }}int main() { #ifdef local freopen("in.txt","r",stdin); #endif // local int t; cin>>t; while(t--) work();}
本文参考了官方题解后写出
- 【Codeforces Gym
- 【codeforces Gym
- codeforces Gym
- codeforces Gym
- codeforces Gym
- codeforces Gym
- Codeforces Gym
- Codeforces Gym
- Codeforces Gym
- codeforces Gym
- Codeforces Gym
- Codeforces Gym
- Codeforces Gym
- [codeforces] Gym
- [codeforces] Gym
- [codeforces] Gym
- [codeforces] Gym
- [codeforces] Gym
- 莫队(bzoj 2038: [2009国家集训队]小Z的袜子(hose))
- 上海家长给孩子择校的福利 虎妈家,一站式教育规划网站
- FPGA学习网站推荐
- 原生ajax上传多图然后显示
- 浅谈zk(3).领导者与追随者如何协作?
- Codeforces Gym
- 我的博客地址
- Vue+NodeJS的跨域请求如何保存session
- [编程之美-15]输入一颗二元查找树将该树转换为它的镜像
- 实习两周来的体验
- Python::判断语句流程控制
- rtop 通过 SSH 监控远程主机
- [算法与数据结构]
- python数据分析与展示(三)