vjudge 题面传送门 首先我们知道斐波那契数列的 lcm 是不太容易计算的,但是它们的 gcd 非常容易计算--\(\gcd(f_x,f_y)=f_{\gcd(x,y)}\),该性质已在我的这篇博客中给出了详细证明,这里就不再赘述了. 考虑怎样将 LCM 转化为 gcd,注意到有个东西叫 Min-Max 容斥,即对于集合 \(S\),\(\max(S)=\sum\limits_{\varnothing\ne T\subseteq S}(-1)^{|T|+1}\min(T)\),该性质同样可以…
[51nod1355]斐波那契的最小公倍数(min-max容斥) 题面 51nod 题解 显然直接算还是没法算的,所以继续考虑\(min-max\)容斥计算. \[lcm(S)=\prod_{T\subset S}gcd(T)^{(-1)^{|T|+1}}\] 而斐波那契数列满足\(gcd(f(a),f(b))=f(gcd(a,b))\), 于是和最小公倍佩尔数一样的类似处理 \[lcm(S)=\prod_{i=1}^{\infty}f(i)^{\sum_{T\subset S}[gcd(T)=…
之前一直没敢做矩阵一类的题目 其实还好吧 推荐看一下 : http://www.cnblogs.com/SYCstudio/p/7211050.html 但是后面的斐波那契 推导不是很懂  前面讲的挺好的 后来看到了 http://blog.csdn.net/flyfish1986/article/details/48014523 相当于  是一个那个东西的k-1次方  而且由于 F(1) = 1 所以直接求k-1次方就可以了 #include<bits/stdc++.h> using nam…
斐波那契数列的定义如下:   F(0) = 0 F(1) = 1 F(n) = F(n - 1) + F(n - 2) (n >= 2)   (1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, 233, 377, ...) 给出n,求F(n),由于结果很大,输出F(n) % 1000000009的结果即可.   Input输入1个数n(1 <= n <= 10^18).Output输出F(n) % 1000000009的结果.Sample Input…
直接斐波那契... #include<stdio.h> #include<queue> #include<string.h> #include<iostream> #include<algorithm> using namespace std; typedef long long LL; const int INF=0x3f3f3f3f; const LL mod=1e9+7; LL a[1010]; int main() { a[1]=1; a…
斐波那契数列的定义如下:   F(0) = 0 F(1) = 1 F(n) = F(n - 1) + F(n - 2) (n >= 2)   (1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, 233, 377, ...) 给出n,求F(n),由于结果很大,输出F(n) % 1000000009的结果即可.   输入 输入1个数n(1 <= n <= 10^18). 输出 输出F(n) % 1000000009的结果. 输入样例 11 输出样例 89解…
Description 给定 \(n\) 个正整数 \(a_1,a_2,...,a_n\),求 \(\text{lcm}(f_{a_1},f_{a_2},...,f_{a_n})\).其中 \(f_i\) 是斐波那契数列第 \(i\) 项. \(n\leq 50000,a_i\leq 10^6\). Sol 首先关于集合 \(S\) 的\(\text{lcm}\)可以用类似\(\text{min-max}\)容斥的式子搞一下,变成跟\(\gcd\)有关: \[ \text{lcm}(T)=\pr…
\(\mathcal{Description}\)   Link.   令 \(f\) 为 \(\text{Fibonacci}\) 数列,给定 \(\{a_n\}\),求: \[\operatorname{lcm}\{f_{a_1},f_{a_2},\cdots,f_{a_n}\}\bmod(10^9+7) \]   \(n\le5\times10^4\),\(a_i\le10^6\). \(\mathcal{Solution}\)   你得知道: \[\gcd(f_i,f_j)=f_{\gc…
传送门 题意 分析 我们发现该数列遵循下列规律: 1 1,2 1,2,2 1,2,2,2,3 1,2,2,2,3,2,3,3 我们令A[i]表示f[i]开始长为f[i-1]的i的最短表示和 那么得到A[i]=A[i-1]+A[i-2]+f[i-2] 那么先预处理出每一段和A[i],i不会超过84 先连续加A[i],对于剩余一段,递归处理,具体见代码 感谢qwb 代码 #include <bits/stdc++.h> using namespace std; #define ll long lo…
#include <iostream> #include <algorithm> using namespace std; typedef long long LL; ; ; struct Matrix { LL v[maxn][maxn]; }; //矩阵间的乘法 Matrix matrix_mul(Matrix A, Matrix B){ Matrix ans; ; i < maxn; i++){ ; j < maxn; j++){ ans.v[i][j] = ;…