洛谷 P2152 [SDOI2009]SuperGCD (高精度)
这道题直接写了我两个多小时……
主要是写高精度的时候还存在着一些小毛病,调了很久
在输入这一块卡了很久。
然后注意这里用while的形式写,不然会炸
最后即使我已经是用的万进制了,但是交上去还是有两个点超时
然后就开始漫长的改进,一直过不了那两个点。
然后突然发现,貌似这道题没有涉及到乘法。
那不就可以直接开10的九次方为一位了(10的9次方是int的最大数量级)
我交上去之后就AC了,全部测试点交上去总和4.6秒
然后原来我一个数开的是2500这么大,因为题目给的10000位,之前用万进制除以4就是2500
然后我现在改成了1200,大约是10000除以9
然后……
交上去总和3秒
直接快了一秒多
而且单个测试点最多也就0.7秒,而这道题的最大限制是2s
看来省空间的时候时间也省了很多
#include<cstdio>
#include<cstring>
#include<algorithm>
#define REP(i, a, b) for(int i = (a); i < (b); i++)
#define _for(i, a, b) for(int i = (a); i <= (b); i++)
using namespace std;
const int MAXN = 1200;
const int base = 1e9;
struct bignum
{
int len, s[MAXN];
bignum() { len = 0; memset(s, 0, sizeof(s)); }
};
bignum operator - (bignum a, const bignum& b)
{
for(int i = a.len; i >= 1; i--)
{
a.s[i] -= b.s[i];
if(a.s[i] < 0) a.s[i+1]--, a.s[i] += base;
}
while(!a.s[a.len] && a.len > 0) a.len--;
return a;
}
bool judge(bignum a, bignum b) //这种写法很简略
{
if(a.len != b.len) return a.len > b.len;
for(int i = a.len; i >= 1; i--)
if(a.s[i] != b.s[i]) return a.s[i] > b.s[i];
return true;
}
bignum operator % (bignum a, bignum b)
{
while(judge(a, b)) a = a - b;
return a;
}
char str[10000 + 5];
void read(bignum& a) //这个输入代码写了好久
{
scanf("%s", str);
reverse(str, str + strlen(str)); //先翻转在说
int& len = a.len = 0;
for(int i = 0, w; i < strlen(str); i++, w *= 10)
{
if(i % 9 == 0) len++, w = 1;
a.s[len] += w * (str[i] - '0');
}
}
void print(bignum a)
{
printf("%d", a.s[a.len]);
for(int i = a.len - 1; i >= 1; i--)
printf("%09d", a.s[i]);
puts("");
}
int main()
{
bignum a, b, c;
read(a); read(b);
while(b.len) //规定len = 0时值为0
{
c = a % b;
a = b;
b = c;
}
print(a);
return 0;
}
洛谷 P2152 [SDOI2009]SuperGCD (高精度)的更多相关文章
- 洛谷 P2152 [SDOI2009]SuperGCD
题意简述 求两个整数a,b的最大公约数0 < a , b ≤ 10 ^ 10000. 题解思路 如果 a % 2 == 0 && b % 2 == 0 gcd(a,b) = gc ...
- 洛谷 2152 [SDOI2009]SuperGCD
Description Sheng bill有着惊人的心算能力,甚至能用大脑计算出两个巨大的数的GCD(最大公约 数)!因此他经常和别人比赛计算GCD.有一天Sheng bill很嚣张地找到了你,并要 ...
- 洛谷P1972 [SDOI2009]HH的项链 题解
[SDOI2009]HH的项链 题目背景 无 题目描述 HH 有一串由各种漂亮的贝壳组成的项链.HH 相信不同的贝壳会带来好运,所以每次散步完后,他都会随意取出一段贝壳,思考它们所表达的含义.HH 不 ...
- BZOJ1880或洛谷2149 [SDOI2009]Elaxia的路线
BZOJ原题链接 洛谷原题链接 显然最长公共路径是最短路上的一条链. 我们可以把最短路经过的边看成有向边,那么组成的图就是一张\(DAG\),这样题目要求的即是两张\(DAG\)重合部分中的最长链. ...
- 洛谷P2148 [SDOI2009]E&D(博弈论)
洛谷题目传送门 先安利蒟蒻仍在施工的博弈论总结 首先根据题目,石子被两两分组了,于是根据SG定理,我们只要求出每一组的SG值再全部异或起来就好啦. 把每一对数看成一个ICG,首先,我们尝试构造游戏的状 ...
- BZOJ1227或洛谷2154 [SDOI2009]虔诚的墓主人
BZOJ原题链接 洛谷原题链接 又是扫描线,题解可看大佬的博客(太懒了不想打) #include<cstdio> #include<algorithm> using names ...
- BZOJ1228或洛谷2148 [SDOI2009]E&D
BZOJ原题链接 洛谷原题链接 完全不会呀.. 写了这题才知道\(SG\)函数原来也能打表找规律... 题解请看大佬的博客 #include<cstdio> using namespace ...
- BZOJ1226或洛谷2157 [SDOI2009]学校食堂
BZOJ原题链接 洛谷原题链接 注意到\(B[i]\)很小,考虑状压\(DP\). 设\(f[i][j][k]\)表示前\(i - 1\)个人已经拿到菜,第\(i\)个人及其后面\(7\)个人是否拿到 ...
- [洛谷P1972][SDOI2009]HH的项链
题目大意:给你一串数字,多次询问区间内数字的种类数 题解:莫队 卡点:洛谷数据加强,开了个$O(2)$ C++ Code: #include <cstdio> #include <a ...
随机推荐
- 深入了解Python--元组
1. 对原元组进行插入 2. 元组的嵌套使用 3. for循环使用嵌套元组实例 4. 命名元组避免对分片混淆
- C语言与java的区别
相同的地方: 1.语法类似:由于Java可以算是从C++发展而来的,因此Java与C语言的语法比较类似 2.编程的熟练程度就是对语言程序库的掌握程度: 从某种程度上来说,编程语言都是由语法和相应的程序 ...
- CodeForcesGym 100548G The Problem to Slow Down You
The Problem to Slow Down You Time Limit: 20000ms Memory Limit: 524288KB This problem will be judged ...
- BZOJ 3674 可持久化并查集加强版(按秩合并版本)
/* bzoj 3674: 可持久化并查集加强版 http://www.lydsy.com/JudgeOnline/problem.php?id=3674 用可持久化线段树维护可持久化数组从而实现可持 ...
- H3C防火墙——回环流量问题(内网终端通过外网IP访问内部服务器)
http://www.bubuko.com/infodetail-1533703.html
- J - Borg Maze
J - Borg Maze 思路:bfs+最小生成树. #include<queue> #include<cstdio> #include<cstring> #in ...
- POJ 1284
想了很久,只想到枚举的方法,估计会超时吧. 原来有这样一条性质:p为素数,则p有phi(p-1)个原根 Orz... #include <iostream> #include <cs ...
- HBase读取代码
HBase读取代码 需要的jar包: activation-1.1.jar aopalliance-1.0.jar apacheds-i18n-2.0.0-M15.jar apacheds-kerbe ...
- stylesheet_link_tag,javascript_include_tag无效?
stylesheet_link_tag,javascript_include_tag无效? http://stackoverflow.com/questions/28241981/rails-4-ex ...
- unity3d 中动画的帧事件
代码事件监听 using UnityEngine; using System.Collections; public class BoxEventScript : MonoBehaviour { vo ...