题目描述
Fibonacci 数列,f(n)=f(n−1)+f(n−2)
前n项为1 1 2 3 5 8 …
给出n,m,需要你计算出满足条件的对数(i,j)的个数,且i<=j。
条件是:
1<=gcd(f(i),f(j))<=n
i,j<=m
由于答案可能很大,所以对1e9+7取模
输入描述:
第一行为一个整数T,表示T组测试数据。
接下来共T行,每行为两个整数n,m。
1<=T<=1e5
1<=n<=1e18
1<=m<=1e6
输出描述:
输出满足条件的对数。
示例1
输入
2
5 5
4 4
输出
15
10
示例2
输入
3
1 6
3 2
3 1
输出
16
3
1
解决思路:
总结结论:
有关 F i b o n a c c i Fibonacci Fibonacci数的 g c d gcd gcd的题目,需要了解一个公式, g c d ( f ( n ) , f ( m ) ) = f ( g c d ( n , m ) ) gcd(f(n),f(m))=f(gcd(n,m)) gcd(f(n),f(m))=f(gcd(n,m))已知条件:
将题目中的要求写成公式形式。
即需要计算: ∑ g = 1 n ∑ i = 1 m ∑ j = i m [ g c d ( f ( i ) , f ( j ) ) = g ] \sum\limits_{g=1}^{n}\sum\limits_{i=1}^{m}\sum\limits_{j=i}^{m}[gcd(f(i),f(j))=g] g=1∑ni=1∑mj=i∑m[gcd(f(i),f(j))=g]
这个公式的含义是: g g g枚举 g c d gcd gcd的时候,看一下有多少对 ( i , j ) (i,j) (i,j)满足条件。
应用上面的公式,转换一下为: ∑ g = 1 n ∑ i = 1 m ∑ j = i m [ f ( g c d ( i , j ) ) = g ] \sum\limits_{g=1}^{n}\sum\limits_{i=1}^{m}\sum\limits_{j=i}^{m}[f(gcd(i,j))=g] g=1∑ni=1∑mj=i∑m[f(gcd(i,j))=g]
观察上面的公式,会发现只有当 g g g和 f ( g c d ( i , j ) ) f(gcd(i,j)) f(gcd(i,j))相等时才会有贡献。
且 F i b o n a c c i Fibonacci Fibonacci增长很快,所以 g g g的可能取值只有最小的 F i b o n a c c i Fibonacci Fibonacci数到小于等于n的最大的 F i b o n a c c i Fibonacci Fibonacci数。
因为正好取的是连续的 F i b o n a c c i Fibonacci Fibonacci数,所以 g c d ( i , j ) gcd(i,j) gcd(i,j)的取值也是连续的。
设小于等于 n n n的最大的 F i b o n a c c i Fibonacci Fibonacci数为 m x mx mx,且是第 c n t cnt cnt个 F i b o n a c c i Fibonacci Fibonacci数。
公式就可以转化为: ∑ g = 1 c n t ∑ i = 1 m ∑ j = i m [ g c d ( i , j ) = g ] \sum\limits_{g=1}^{cnt}\sum\limits_{i=1}^{m}\sum\limits_{j=i}^{m}[gcd(i,j)=g] g=1∑cnti=1∑mj=i∑m[gcd(i,j)=g]
需要计算的是 i ≤ j i \le j i≤j的对数,会发现等价于 j ≤ i j \le i j≤i的对数。
公式变为: ∑ g = 1 c n t ∑ i = 1 m ∑ j = 1 i [ g c d ( i , j ) = g ] \sum\limits_{g=1}^{cnt}\sum\limits_{i=1}^{m}\sum\limits_{j=1}^{i}[gcd(i,j)=g] g=1∑cnti=1∑mj=1∑i[gcd(i,j)=g]
g c d ( i , j ) = g gcd(i,j)=g gcd(i,j)=g意味着 i i i和 j j j都是 g g g的倍数,所以对 i , j i,j i,j除以 g g g,变换枚举上限。
公式为: ∑ g = 1 c n t ∑ i = 1 ⌊ m g ⌋ ∑ j = 1 i [ g c d ( i , j ) = 1 ] \sum\limits_{g=1}^{cnt}\sum\limits_{i=1}^{\left \lfloor \frac{m}{g}\right \rfloor}\sum\limits_{j=1}^{i}[gcd(i,j)=1] g=1∑cnti=1∑⌊gm⌋j=1∑i[gcd(i,j)=1]
发现里面就是欧拉函数的形式了。
公式为: ∑ g = 1 c n t ∑ i = 1 ⌊ m g ⌋ φ ( i ) \sum\limits_{g=1}^{cnt}\sum\limits_{i=1}^{\left \lfloor \frac{m}{g}\right \rfloor}\varphi (i) g=1∑cnti=1∑⌊gm⌋φ(i)
c n t cnt cnt很小,里面就是欧拉函数的前缀和,直接预处理即可。
思路里程:
看到 F i b o n a c c i Fibonacci Fibonacci的最大公约数,就试着套用上面的公式求解,这个公式可以降低枚举范围,本来枚举的是 F i b o n a c c i Fibonacci Fibonacci的范围,降低之后只需要枚举下标的范围了。
对数个数计算的几种方法。
第一种:和文章中一样,如果可以交换顺序的话,就用那种容易计算的就行。
第二种:可以容斥计算,小于等于的个数等于随便取的个数加上相等的对数个数除以2。
实现过程:
需要暴力计算出 c n t cnt cnt(小于等于 n n n的最大的 F i b o n a c c i Fibonacci Fibonacci数的下标)
欧拉筛预处理欧拉函数,再对欧拉函数求前缀和。
取模时处处取模,防止溢出。
代码
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
#ifdef DOIDO_DEBUG
#define debug(x) cout<< #x << " --> " << (x) << "\n";
#define debug_pair(x, y) cout<< #x << " " << #y \
<< " --> " << (x) << " " << (y) << "\n";
#else
#define debug(x)
#define debug_pair(x, y)
#endif
#define endl '\n'
#define N int(1e6 + 10)
#define MOD int(1e9 + 7)
int cnt;
int phi[N];
bool vis[N];
int prime[N];
ll pre_phi[N];
void init_prime(const int& n) {
phi[1] = 1;
for (int i = 2; i <= n; i++) {
if (!vis[i]) {
phi[i] = i - 1;
prime[++cnt] = i;
}
for (int j = 1; i <= n / prime[j]; j++) {
vis[i * prime[j]] = true;
if (i % prime[j] == 0) {
phi[i * prime[j]] = phi[i] * (prime[j]);
break;
}
phi[i * prime[j]] = phi[i] * (prime[j] - 1);
}
}
ll pre = 0;
for (int i = 1; i <= n; i++) {
// 处处取模,防止溢出
pre_phi[i] = (pre_phi[i - 1] + phi[i]) % MOD;
}
}
int main() {
ios::sync_with_stdio(false);
cin.tie(nullptr);
// 预处理欧拉函数前缀和
init_prime(1e6);
// 计算Fibonacci数,大于最大值就跳出。
static ll fi[N] = { 0, 1, 1 };
for (int i = 3; i <= 100; i++) {
fi[i] = fi[i - 1] + fi[i - 2];
if (fi[i] > 1e18) break;
}
int T;
cin >> T;
for (int j = 1; j <= T; j++) {
ll n, m;
cin >> n >> m;
int pos = 0;
for (int k = 1; k <= 100; k++) {
if (fi[k] > n) {
pos = k;
// 找到大于n的位置跳出即可。
break;
}
}
// 由于找的是大于n的最小位置,所以需要位置减一
pos--;
ll ans = 0;
for (int g = 1; g <= pos; g++) {
ans = (ans + pre_phi[m / g]) % MOD;
}
cout << ans << endl;
}
return 0;
}