[hdu5351]找规律,大整数模板
题意:f(1)="a",f(2)="b",f(i)=f(i-1)+f(i-2),"+"表示连接符。给定n,m,求f(n)的前m个字符的“next值”。
思路:并不知道如何一步步推出结论,只能打个表找找规律了:找到最小的i使得f(i)>m+1,则答案就是m-f(i-2)。然后就是大整数模板了
|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
|
#include <map>#include <set>#include <cmath>#include <ctime>#include <deque>#include <queue>#include <vector>#include <cstdio>#include <string>#include <cstdlib>#include <cstring>#include <iostream>#include <algorithm>using namespace std;#define X first#define Y second#define pb push_back#define mp make_pair#define all(a) (a).begin(), (a).end()#define fillchar(a, x) memset(a, x, sizeof(a))typedef long long ll;typedef pair<int, int> pii;typedef unsigned long long ull;#ifndef ONLINE_JUDGEvoid RI(vector<int>&a,int n){a.resize(n);for(int i=0;i<n;i++)scanf("%d",&a[i]);}void RI(){}void RI(int&X){scanf("%d",&X);}template<typename...R>void RI(int&f,R&...r){RI(f);RI(r...);}void RI(int*p,int*q){int d=p<q?1:-1;while(p!=q){scanf("%d",p);p+=d;}}void print(){cout<<endl;}template<typename T>void print(const T t){cout<<t<<endl;}template<typename F,typename...R>void print(const F f,const R...r){cout<<f<<", ";print(r...);}template<typename T>void print(T*p, T*q){int d=p<q?1:-1;while(p!=q){cout<<*p<<", ";p+=d;}cout<<endl;}#endiftemplate<typename T>bool umax(T&a, const T&b){return b<=a?false:(a=b,true);}template<typename T>bool umin(T&a, const T&b){return b>=a?false:(a=b,true);}template<typename T>void V2A(T a[],const vector<T>&b){for(int i=0;i<b.size();i++)a[i]=b[i];}template<typename T>void A2V(vector<T>&a,const T b[]){for(int i=0;i<a.size();i++)a[i]=b[i];}const double PI = acos(-1.0);const int INF = 1e9 + 7;/* -------------------------------------------------------------------------------- */struct BigInt { const static int maxI = 1e8; const static int Len = 8; typedef vector<int> vi; typedef long long LL; vi num; bool symbol; BigInt() { num.clear(); symbol = 0; } BigInt(int x) { symbol = 0; if (x < 0) { symbol = 1; x = -x; } num.push_back(x % maxI); if (x >= maxI) num.push_back(x / maxI); } BigInt(bool s, vi x) { symbol = s; num = x; } BigInt(char s[]) { int len = strlen(s), x = 1, sum = 0, p = s[0] == '-'; symbol = p; for (int i = len - 1; i >= p; i--) { sum += (s[i] - '0') * x; x *= 10; if (x == 1e8 || i == p) { num.push_back(sum); sum = 0; x = 1; } } while (num.back() == 0 && num.size() > 1) num.pop_back(); } void push(int x) { num.push_back(x); } BigInt abs() const { return BigInt(false, num); } bool smaller(const vi &a, const vi &b) const { if (a.size() != b.size()) return a.size() < b.size(); for (int i = a.size() - 1; i >= 0; i--) { if (a[i] != b[i]) return a[i] < b[i]; } return 0; } bool operator < (const BigInt &p) const { if (symbol && !p.symbol) return true; if (!symbol && p.symbol) return false; if (symbol && p.symbol) return smaller(p.num, num); return smaller(num, p.num); } bool operator > (const BigInt &p) const { return p < *this; } bool operator == (const BigInt &p) const { return !(p < *this) && !(*this < p); } bool operator >= (const BigInt &p) const { return !(*this < p); } bool operator <= (const BigInt &p) const { return !(p < *this); } vi add(const vi &a, const vi &b) const { vi c; c.clear(); int x = 0; for (int i = 0; i < a.size(); i++) { x += a[i]; if (i < b.size()) x += b[i]; c.push_back(x % maxI); x /= maxI; } for (int i = a.size(); i < b.size(); i++) { x += b[i]; c.push_back(x % maxI); x /= maxI; } if (x) c.push_back(x); while (c.back() == 0 && c.size() > 1) c.pop_back(); return c; } vi sub(const vi &a, const vi &b) const { vi c; c.clear(); int x = 1; for (int i = 0; i < b.size(); i++) { x += maxI + a[i] - b[i] - 1; c.push_back(x % maxI); x /= maxI; } for (int i = b.size(); i < a.size(); i++) { x += maxI + a[i] - 1; c.push_back(x % maxI); x /= maxI; } while (c.back() == 0 && c.size() > 1) c.pop_back(); return c; } vi mul(const vi &a, const vi &b) const { vi c; c.resize(a.size() + b.size()); for (int i = 0; i < a.size(); i++) { for (int j = 0; j < b.size(); j++) { LL tmp = (LL)a[i] * b[j] + c[i + j]; c[i + j + 1] += tmp / maxI; c[i + j] = tmp % maxI; } } while (c.back() == 0 && c.size() > 1) c.pop_back(); return c; } vi div(const vi &a, const vi &b) const { vi c(a.size()), x(1, 0), y(1, 0), z(1, 0), t(1, 0); y.push_back(1); for (int i = a.size() - 1; i >= 0; i--) { z[0] = a[i]; x = add(mul(x, y), z); if (smaller(x, b)) continue; int l = 1, r = maxI - 1; while (l < r) { int m = (l + r + 1) >> 1; t[0] = m; if (smaller(x, mul(b, t))) r = m - 1; else l = m; } c[i] = l; t[0] = l; x = sub(x, mul(b, t)); } while (c.back() == 0 && c.size() > 1) c.pop_back(); return c; } BigInt operator + (const BigInt &p) const { if (!symbol && !p.symbol) return BigInt(false, add(num, p.num)); if (!symbol && p.symbol) { return *this >= p.abs() ? BigInt(false, sub(num, p.num)) : BigInt(true, sub(p.num, num)); } if (symbol && !p.symbol) { return (*this).abs() > p ? BigInt(true, sub(num, p.num)) : BigInt(false, sub(p.num, num)); } return BigInt(true, add(num, p.num)); } BigInt operator - (const BigInt &p) const { return *this + BigInt(!p.symbol, p.num); } BigInt operator * (const BigInt &p) const { BigInt res(symbol ^ p.symbol, mul(num, p.num)); if (res.symbol && res.num.size() == 1 && res.num[0] == 0) res.symbol = false; return res; } BigInt operator / (const BigInt &p) const { if (p == BigInt(0)) return p; BigInt res(symbol ^ p.symbol, div(num, p.num)); if (res.symbol && res.num.size() == 1 && res.num[0] == 0) res.symbol = false; return res; } BigInt operator % (const BigInt &p) const { return *this - *this / p * p; } void show() const { if (symbol) putchar('-'); printf("%d", num[num.size() - 1]); for (int i = num.size() - 2; i >= 0; i--) { printf("%08d", num[i]); } putchar('\n'); } int TotalDigit() const { int x = num[num.size() - 1] / 10, t = 1; while (x) { x /= 10; t++; } return t + (num.size() - 1) * Len; }};const int md = 258280327;BigInt f[1234];char s[1000];void pre_init() { f[0] = 0; f[1] = f[2] = 1; for (int i = 3; i <= 1111; i ++) { f[i] = f[i - 1] + f[i - 2]; }}int main() {#ifndef ONLINE_JUDGE freopen("in.txt", "r", stdin); //freopen("out.txt", "w", stdout);#endif // ONLINE_JUDGE pre_init(); int T; cin >> T; while (T --) { int n; cin >> n; scanf("%s", s); BigInt buf(s); int L = 1, R = 1111; while (L < R) { int M = (L + R) >> 1; if (f[M] > buf + 1) R = M; else L = M + 1; } ((buf - f[L - 2]) % md).show(); } return 0;} |
[hdu5351]找规律,大整数模板的更多相关文章
- HDU 4910 Problem about GCD 找规律+大素数判断+分解因子
Problem about GCD Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others ...
- [带符号大整数模板]vector版
#include <iostream> #include <cstdio> #include <vector> #include <cstring> u ...
- OpenJudge 2737 大整数除法
链接地址:http://bailian.openjudge.cn/practice/2737/ 题目: 总时间限制: 1000ms 内存限制: 65536kB 描述 求2个大的正整数相除的商 输入 第 ...
- vijos - P1447开关灯泡 (大数模板 + 找规律 + 全然数 + python)
P1447开关灯泡 Accepted 标签:CSC WorkGroup III[显示标签] 描写叙述 一个房间里有n盏灯泡.一開始都是熄着的,有1到n个时刻.每一个时刻i,我们会将i的倍数的灯泡改变状 ...
- Pollard-Rho大整数拆分模板
随机拆分,简直机智. 关于过程可以看http://wenku.baidu.com/link?url=JPlP8watmyGVDdjgiLpcytC0lazh4Leg3s53WIx1_Pp_Y6DJTC ...
- HDU - 4722 Good Numbers 【找规律 or 数位dp模板】
If we sum up every digit of a number and the result can be exactly divided by 10, we say this number ...
- hdu 5051 找规律?+大trick
http://acm.hdu.edu.cn/showproblem.php?pid=5051 打表找规律 据说是http://zh.wikipedia.org/wiki/%E6%9C%AC%E7%A6 ...
- 大整数类BIGN的设计与实现 C++高精度模板
首先感谢刘汝佳所著的<算法竞赛入门经典>. 众所周知,C++中储存能力最大的unsigned long long 也是有着一个上限,如果我们想计算非常大的整数时,就不知所措了,所以,我写了 ...
- [国家集训队]整数的lqp拆分 数学推导 打表找规律
题解: 考场上靠打表找规律切的题,不过严谨的数学推导才是本题精妙所在:求:$\sum\prod_{i=1}^{m}F_{a{i}}$ 设 $f(i)$ 为 $N=i$ 时的答案,$F_{i}$ 为斐波 ...
随机推荐
- PHP函数:array_rand
array_rand() - 从数组中随机取出一个或多个单元. 说明: array_rand ( array $array [, int $num = 1 ] ) : mixed 参数: array ...
- MySQL 50题练习
表名和字段 –1.学生表 Student(s_id,s_name,s_birth,s_sex) –学生编号,学生姓名, 出生年月,学生性别 –2.课程表 Course(c_id,c_name,t_id ...
- ubuntu16.04-交叉编译opencv3.4.6
0.前言 在要移植opecv和SeetaFaceEngine-master到ARM板子上运行的所有步骤之前,有几点需要注意的: 查看板子运行的Kernel版本 交叉编译工具链的gcc版本,关键就是工具 ...
- Java IO 流--FileUtils 工具类封装
IO流的操作写多了,会发现都已一样的套路,为了使用方便我们可以模拟commosIo 封装一下自己的FileUtils 工具类: 1.封装文件拷贝: 文件拷贝需要输入输出流对接,通过输入流读取数据,然后 ...
- Mac安装多版本JDK
0. 配置JDK环境 安装完成之后,配置.bash_profile文件 使用source ./bash_profile激活 2. 如何切换默认的jdk? 使用java -version就可以看默认版本 ...
- Python模块---制作新冠疫情世界地图()
目录 pyecharts模块 简介 安装pyecharts 测试pyecharts模块 pyecharts实战:绘制新冠肺炎疫情地图 需求分析 请求数据 提取数据 处理数据 制作可视化地图 设置可视化 ...
- http请求返回的数字代表的含义
一些常见的状态码为: 200 - 服务器成功返回网页 404 - 请求的网页不存在 503 - 服务器超时 下面提供 HTTP 状态码的完整列表.点击链接可了解详情.您也可以访问 HTTP 状态码上的 ...
- Docker虚拟化管理:30分钟教你学会用Docker
关于Docker的官方介绍网上太多了我就不贴了,就实际体验来说Docker可以极大的简化环境搭建及服务部署的操作流程,大大降低部署的时间成本,解放你的双手. 本文不会深入讲解Docker底层架构及运行 ...
- 【React踩坑记五】React项目中引入并使用react-ace代码编辑插件(自定义列表提示)
最近有一个引入sql编辑器插件的需求,要求代码高亮显示,代码智能提示,以及支持自定义代码提示列表等功能.中途在自定义代码提示列表中由于没有相关demo,所以踩了一些坑,遂将其整理如下,以便日后查看. ...
- mysql查询连接数
最近公司的测试服务器数据库经常是连接爆满,几次加大了依旧满了. 明明只有几个人在用这个数据库,但是连接数到了三四百.于是就想是谁一直开着连接不释放,于是写了个SQL查了下连接使用情况. SQL如下: ...