bzoj 2813: 奇妙的Fibonacci 数学+线性筛
来源:互联网 发布:网页游戏传奇源码下载 编辑:程序博客网 时间:2024/05/16 07:25
题意
Fibonacci数列是这样一个数列:
F1 = 1, F2 = 1, F3 = 2 …
Fi = Fi-1 + Fi-2 (当 i >= 3)
pty忽然对这个古老的数列产生了浓厚的兴趣,他想知道:对于某一个Fibonacci数Fi,
有多少个Fj能够整除Fi (i可以等于j),他还想知道所有j的平方之和是多少。答案模1e9+7。
第一行一个整数Q,表示Q个询问。
第二行四个整数:Q1, A, B, C
第i个询问Qi = (Qi-1 * A + B) mod C + 1(当i >= 2)
Q <= 3*10^6,C <= 10^7,A <= 10^7,B <= 10^7,1 <= Q1<= C
分析
这题利用到了斐波那契数列一些比较神奇的性质。
定理一:
证明:
…
定理二:
证明:
设
当
当
…
当
所以
定理三:
证明:
然后通过归纳不难发现
要想
然后线性筛爱怎么筛怎么筛。
因为
代码
#include<iostream>#include<cstdio>#include<cstdlib>#include<cstring>#include<algorithm>using namespace std;typedef long long LL;const int MOD=1000000007;const int N=10000005;int pow[N],num[N],prime[N],a[N],low[N],tot;LL sum[N];bool not_prime[N];void get_prime(int n){ for (int i=2;i<=n;i++) { if (!not_prime[i]) prime[++tot]=i,low[i]=i,a[i]=1,pow[i]=i; for (int j=1;j<=tot&&i*prime[j]<=n;j++) { not_prime[i*prime[j]]=1; low[i*prime[j]]=prime[j]; if (i%prime[j]==0) { a[i*prime[j]]=a[i]+1; pow[i*prime[j]]=pow[i]*prime[j]; break; } a[i*prime[j]]=1; pow[i*prime[j]]=prime[j]; } } sum[1]=num[1]=1; for (int i=2;i<=n;i++) num[i]=num[i/pow[i]]*(a[i]+1),sum[i]=(LL)sum[i/low[i]]*low[i]*low[i]+sum[i/pow[i]];}int main(){ get_prime(10000000); LL T,q,A,B,C; scanf("%lld%lld%lld%lld%lld",&T,&q,&A,&B,&C); LL ans1=0,ans2=0; while (T--) { (ans1+=num[q])%=MOD;(ans2+=sum[q])%=MOD; if (q&1) ans1++,ans2+=4; q=(q*A+B)%C+1; } printf("%lld\n%lld",ans1%MOD,ans2%MOD);}
- bzoj 2813: 奇妙的Fibonacci 数学+线性筛
- BZOJ 2813 奇妙的Fibonacci 线性筛
- BZOJ 2813: 奇妙的Fibonacci 线性筛
- BZOJ 2813: 奇妙的Fibonacci 线性筛
- BZOJ 2813 奇妙的Fibonacci
- bzoj 2813 奇妙的Fibonacci (线性筛求约束平方和)
- 【bzoj2813】奇妙的Fibonacci 线性筛法
- [BZOJ2813]奇妙的Fibonacci(线性筛)
- 【欧拉筛/线性筛】BZOJ2813 奇妙的Fibonacci
- [bzoj2813]奇妙的Fibonacci数列 (数论+线性筛)
- 【BZOJ2813】奇妙的Fibonacci
- 【bzoj2813】奇妙的Fibonacci
- [BZOJ2813]奇妙的Fibonacci(数论)
- leetcode459奇妙做法的数学解释
- BZOJ 1002 - 奇妙的题目 + 高精度
- BZOJ2813--奇妙的斐波那契【线性筛】【斐波那契数列】
- 【bzoj 十连测】[noip2016十连测第八场]Problem A: 神炎皇(数学相关:线性筛)
- 非负矩阵分解:数学的奇妙力量
- c语言实现单链表面试题之进阶篇
- Error running app: Instant Run requires 'Tools | Android | Enable ADB integration' to be enabled.
- 3 Steps(二分图初步)
- 合并分支到master上
- freemarker的基本使用和常用标签
- bzoj 2813: 奇妙的Fibonacci 数学+线性筛
- 第三周-项目2
- [bzoj-1047][HAOI2007]理想的正方形 题解
- hbase
- 线程的使用及相关代码
- HttpClient的用法
- JavaWeb知识体系梳理和学习资料汇总
- XML解析(Dom /DOM4J)
- C++中两个类简单互相包含会形成死循环