Aggregated Counting-----hdu5439(2015 长春网络赛 找规律)
2015-09-23 10:20
519 查看
#include<stdio.h> #include<string.h> #include<iostream> #include<math.h> #include<stdlib.h> #include<algorithm> using namespace std; #define N 500000 #define mod 1000000007 __int64 a , b , ans ; void Init() { int i; a[1] = b[1] = 1; a[2] = 2; b[2] = 3; for(i=1; i<N; i++) { a[i] = lower_bound(b+1, b+i, i) - b; b[i] = b[i-1] + a[i]; } ans[1] = 1; for(i=2; i<N; i++) { ans[i] = ans[i-1] + (b[i]-b[i-1])*(b[i-1]+1 + b[i]) / 2 % mod * i % mod; } } int main() { int T, n, i; Init(); scanf("%d", &T); while(T--) { scanf("%d", &n); int pos = lower_bound(b+1, b+N+1, n) - b; __int64 Ans = ans[pos-1]; for(i=b[pos-1]+1; i<=n; i++) Ans = (Ans + (__int64)i*pos) % mod; printf("%I64d\n", Ans); } return 0; }
View Code
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5439
函数lower_bound()在first和last中的[b]前闭后开区间进行二分查找,返回大于或等于val的第一个元素位置。如果所有元素都小于val,则返回last的位置[/b]
pos = lower_bound(a, a+n, k) - a; 是在a中>=k的第一个数的位置;
题意:刚开始给一个1,序列a是由a[i]个i组成,最后就变成了1,2,2,3,3,4,4,4,5,5,5.......
就是求Ans: 即n在a中最后出现的位置m,m在a中的最后出现的位置Ans;
我们用b【i】表示i在a中出现的最后位置;
用ans【i】表示b【i】在a中出现的最后位置;
#include<stdio.h> #include<string.h> #include<iostream> #include<math.h> #include<stdlib.h> #include<algorithm> using namespace std; #define N 500000 #define mod 1000000007 __int64 a , b , ans ; void Init() { int i; a[1] = b[1] = 1; a[2] = 2; b[2] = 3; for(i=1; i<N; i++) { a[i] = lower_bound(b+1, b+i, i) - b; b[i] = b[i-1] + a[i]; } ans[1] = 1; for(i=2; i<N; i++) { ans[i] = ans[i-1] + (b[i]-b[i-1])*(b[i-1]+1 + b[i]) / 2 % mod * i % mod; } } int main() { int T, n, i; Init(); scanf("%d", &T); while(T--) { scanf("%d", &n); int pos = lower_bound(b+1, b+N+1, n) - b; __int64 Ans = ans[pos-1]; for(i=b[pos-1]+1; i<=n; i++) Ans = (Ans + (__int64)i*pos) % mod; printf("%I64d\n", Ans); } return 0; }
#include<stdio.h> #include<string.h> #include<iostream> #include<math.h> #include<stdlib.h> #include<algorithm> using namespace std; #define N 500000 #define mod 1000000007 __int64 a , b , ans ; void Init() { int i; a[1] = b[1] = 1; a[2] = 2; b[2] = 3; for(i=1; i<N; i++) { a[i] = lower_bound(b+1, b+i, i) - b; b[i] = b[i-1] + a[i]; } ans[1] = 1; for(i=2; i<N; i++) { ans[i] = ans[i-1] + (b[i]-b[i-1])*(b[i-1]+1 + b[i]) / 2 % mod * i % mod; } } int main() { int T, n; Init(); scanf("%d", &T); while(T--) { scanf("%d", &n); int pos = upper_bound(b+1, b+N+1, n) - b - 1; __int64 Ans = (ans[pos] + (b[pos] + 1 + n)*(n - b[pos])/2 % mod * (pos+1)) % mod; printf("%I64d\n", Ans); } return 0; }
相关文章推荐
- HTTP返回的状态码
- linux下面测试网络带宽 bmon iftop nload
- 从数据长征和网络协议对TCP/IP,http,socket 进行概念梳理
- 使用nginx搭建http代理服务器
- UI:网络请求
- TCP-socket发送结构体类型数据
- ASIHTTPRequest详解 [经典3]
- 浅谈HTTP中Get与Post的区别
- AFNetworking判断当前手机的网络状态
- ASIHTTPRequest详解 [经典]
- ASIHTTPRequest-插件的使用
- C# 程序员最常犯的 10 个错误http://www.oschina.net/translate/top-10-mistakes-that-c-sharp-programmers-make
- Linux管道编程实例 http://blog.csdn.net/pcliuguangtao/article/details/6453743
- 通过QQ查看对方的ip地址
- RHCSA 系列(八): 加固 SSH,设定主机名及启用网络服务
- 关于HTTPS://不能访问的问题
- TCP/IP协议详解
- TCP/IP 寻址(二)
- hiho 1228 : Mission Impossible 6,北京网络赛B题
- 【计算机网络】:经典计算题整理