【BZOJ2144】Throw 数论
题目大意
给你三个数\(a,b,c\),每次你可以选择一个数\(s_1\),再选择一个数\(s_2\),把\(s_1\)变成\(2s_2-s_1\),但要求\(s_3\)不在\(s_1\)到\(2s_2-s_1\)之间。
再给你三个数\(x,y,z\),问你是否能把\(a,b,c\)变成\(x,y,z\)。
\(|a|,|b|,|c|,|x|,|y|,|z|\leq {10}^9\)
题解
首先三个数\(a,b,c(a<b<c)\)只能有三种转移:
\begin{cases}
(2a-b,a,c)\\
(a,c,2c-b)\\
(b,2b-a,c)~~~~~(c-b>b-a)\\
(a,2b-c,b)~~~~~(c-b<b-a)
\end{cases}
\]
可以发现,这些状态构成了一棵树。
转移可以用辗转相除法加速。
直接暴力往上跳找LCA即可。
时间复杂度:\(O(\log (\max_a-\min_a))\)
代码
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<ctime>
#include<utility>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
void sort(ll &a,ll &b,ll &c)
{
if(b<a)
swap(a,b);
if(c<a)
swap(a,c);
if(c<b)
swap(b,c);
}
void open(const char *s)
{
#ifdef DEBUG
char str[100];
sprintf(str,"%s.in",s);
freopen(str,"r",stdin);
sprintf(str,"%s.out",s);
freopen(str,"w",stdout);
#endif
}
ll gcd(ll a,ll b)
{
return b?gcd(b,a%b):a;
}
ll getmid(ll a,ll b,ll c)
{
while(b-a!=c-b)
if(b-a>c-b)
{
ll j=(b-a)/(c-b);
ll k=c-b;
c-=j*k;
b-=j*k;
if(b==a)
{
b+=k;
c+=k;
}
}
else
{
ll j=(c-b)/(b-a);
ll k=b-a;
a+=j*k;
b+=j*k;
if(b==c)
{
b-=k;
a-=k;
}
}
return b;
}
pll a1[110];
int d1[110];
pll a2[110];
int d2[110];
int cnt1,cnt2;
int main()
{
open("c");
ll a,b,c,x,y,z;
scanf("%lld%lld%lld%lld%lld%lld",&a,&b,&c,&x,&y,&z);
sort(a,b,c);
sort(x,y,z);
if(gcd(b-a,c-b)!=gcd(y-x,z-y))
{
printf("NO\n");
return 0;
}
if(getmid(a,b,c)!=getmid(x,y,z))
{
printf("NO\n");
return 0;
}
printf("YES\n");
pll s1(b-a,c-b);
pll s2(y-x,z-y);
cnt1=cnt2=0;
d1[0]=d2[0]=0;
while(s1.first!=s1.second)
{
a1[++cnt1]=s1;
if(s1.first>s1.second)
{
d1[cnt1]=s1.first/s1.second;
s1.first%=s1.second;
if(!s1.first)
{
s1.first=s1.second;
d1[cnt1]--;
}
}
else
{
d1[cnt1]=s1.second/s1.first;
s1.second%=s1.first;
if(!s1.second)
{
s1.second=s1.first;
d1[cnt1]--;
}
}
}
a1[++cnt1]=s1;
d1[cnt1]=0;
reverse(a1+1,a1+cnt1+1);
reverse(d1+1,d1+cnt1+1);
while(s2.first!=s2.second)
{
a2[++cnt2]=s2;
if(s2.first>s2.second)
{
d2[cnt2]=s2.first/s2.second;
s2.first%=s2.second;
if(!s2.first)
{
s2.first=s2.second;
d2[cnt2]--;
}
}
else
{
d2[cnt2]=s2.second/s2.first;
s2.second%=s2.first;
if(!s2.second)
{
s2.second=s2.first;
d2[cnt2]--;
}
}
}
a2[++cnt2]=s2;
d2[cnt2]=0;
reverse(a2+1,a2+cnt2+1);
reverse(d2+1,d2+cnt2+1);
ll ans=0;
int i,j;
for(i=1;i<=cnt1;i++)
d1[i]+=d1[i-1];
for(i=1;i<=cnt2;i++)
d2[i]+=d2[i-1];
i=cnt1,j=cnt2;
while(a1[i-1]!=a2[j-1])
if(d1[i-1]>d2[j-1])
{
ans+=d1[i]-d1[i-1];
i--;
}
else
{
ans+=d2[j]-d2[j-1];
j--;
}
if(a1[i].first==a2[j].first||a1[i].second==a2[j].second)
{
if(d1[i]>d2[j])
ans+=d1[i]-d2[j];
else
ans+=d2[j]-d1[i];
}
else
ans+=d1[i]-d1[i-1]+d2[j]-d2[j-1];
printf("%lld\n",ans);
return 0;
}
【BZOJ2144】Throw 数论的更多相关文章
- 数学概念——I - 数论,线性方程
I - 数论,线性方程 Time Limit:1000MS Memory Limit:32768KB 64bit IO Format:%I64d & %I64u Submit ...
- 浅谈Java的throw与throws
转载:http://blog.csdn.net/luoweifu/article/details/10721543 我进行了一些加工,不是本人原创但比原博主要更完善~ 浅谈Java异常 以前虽然知道一 ...
- C++异常处理:try,catch,throw,finally的用法
写在前面 所谓异常处理,即让一个程序运行时遇到自己无法处理的错误时抛出一个异常,希望调用者可以发现处理问题. 异常处理的基本思想是简化程序的错误代码,为程序键壮性提供一个标准检测机制. 也许我们已经使 ...
- Codeforces Round #382 Div. 2【数论】
C. Tennis Championship(递推,斐波那契) 题意:n个人比赛,淘汰制,要求进行比赛双方的胜场数之差小于等于1.问冠军最多能打多少场比赛.题解:因为n太大,感觉是个构造.写写小数据, ...
- java中的throw与throws的区别
什么时运行时异常?什么是非运行时异常? 通俗的讲: 运行时异常:就是编译通过,运行时就崩了,比如数组越界. 非运行时异常:就是编译不通过,这时就得必须去处理了.不然就没法运行了. 全面的讲: Thro ...
- js 利用throw 写的一个小程序
在下边的小程序中比较特殊的是使用isNaN()函数判断一个参数是不是数字, <!DOCTYPE html> <!DOCTYPE html> <html> <h ...
- C++异常处理: try,catch,throw,finally的用法
写在前面 所谓异常处理,即让一个程序运行时遇到自己无法处理的错误时抛出一个异常,希望调用者可以发现处理问题. 异常处理的基本思想是简化程序的错误代码,为程序键壮性提供一个标准检测机制. 也许我们已经使 ...
- Java 中 手动抛出异常: throw new Exception("错误信息") 错误信息的获得
当然需要先用try catch捕获,但注意new Exception("")括号里的字符串其实是异常原因,所以获取是要用ex.getCause().getMessage() int ...
- NOIP2014 uoj20解方程 数论(同余)
又是数论题 Q&A Q:你TM做数论上瘾了吗 A:没办法我数论太差了,得多练(shui)啊 题意 题目描述 已知多项式方程: a0+a1x+a2x^2+..+anx^n=0 求这个方程在[1, ...
随机推荐
- android linux 传文件
EStrongs File Explorer 即: Es文件浏览器 网络 -> 远程管理器 设置 ->设置根目录 linux 使用浏览器访问即可.
- 01 Django REST Framework 介绍
01-Django REST Framework的介绍 Django REST框架是一个用于构建Web API的强大而灵活的工具包. 您可能希望使用REST框架的一些原因: 1. Web可浏览API对 ...
- 周末时间学习Linux
大家都是如何度过周末时光的呢?好多人都认为一周的工作后要好好休息下,于是在家疯狂的补觉,刷剧,打游戏,自我觉得很是正常,工作几天了,休息下不是当然嘛.是的,休息下很正常,但是把周末的时光都用到这些东西 ...
- vue: WebStorm设置快速编译运行
WebSorm是一款优秀的前端开发工具,而Vue项目可以使用Node进行编译运行,平常我们可以通过命令行部署项目进行调试. 本文介绍设置Webstorm进行快速部署Vue项目. 第一步 点击启动快捷按 ...
- CMMI摘要
CMMI_百度百科https://baike.baidu.com/item/CMMI CMMI分为哪几个等级?CMMI等级介绍_百度经验https://jingyan.baidu.com/articl ...
- apach ab 安装时的错误
1.cmd进入bin下:执行httpd -k install 需要管理员 登陆安装 我们打开conf文件夹,找到httpd.conf,修改如下内容,让serverroot指向你的安装位置: Defin ...
- [转帖]SAP一句话入门:Plant Maintenance
SAP一句话入门:Plant Maintenance http://blog.vsharing.com/MilesForce/A618273.html PM就是Plant Maintenance(本文 ...
- macbookpro 以及 surface 的技术规格
macbookpro 13.3 英寸 (对角线) LED 背光显示屏 (采用 IPS 技术):初始分辨率 x ( ppi),支持数百万色彩 15.4 英寸 (对角线) LED 背光显示屏 (采用 IP ...
- Day 5-3 多态与多态性
多态与多态性 鸭子类型 多态与多态性 多态:一类事物有多种形态.比如,动物有多种形态,人,狗,猪,豹子.水也有多种形态,冰,雪,水蒸气. #多态:同一类事物的多种形态 import abc class ...
- MyBatis源码分析1 参数映射分析
首先我们拿出之前的代码,在如图位置打上断点,开始调试 我们规定了一个mapper接口,而调用了mapper接口的getEmpByIdAndLastName,我们并没有实现这个接口,这是因为Mybati ...