华华给月月出题

传送门

题意:

求题中所给式子的答案。

思路:

题目唯一的难点在于想到 f ( x ) = x n f(x) = x^n f(x)=xn是一个完全积性函数,下面简单证明一下:
设 有 f ( x ) = x n , f ( y ) = y n 且 x ≤ n 并 且 y ≤ n , 则 : f ( x ∗ y ) = ( x ∗ y ) n = x n ∗ y n = f ( x ) ∗ f ( y ) 。 设有f(x) = x^n,f(y) = y^n且x\leq n 并且 y\leq n,则:f(x*y) = (x*y)^n = x^n*y^n = f(x)*f(y)。 f(x)=xn,f(y)=ynxnyn,f(xy)=(xy)n=xnyn=f(x)f(y)
证毕!
然后用欧拉筛求积性函数就好啦。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define endl '\n'
const int mod = 1e9+7;
ll qpow(ll a,ll b)
{
	ll res = 1;
	while(b)
	{
		if(b&1)res = res*a%mod;
		a = a*a%mod;
		b>>=1;
	}
	return res;
}
ll vis[13000010],pri[13000010];
int cnt;
void oula(int n)
{
	for(int i = 2; i <= n; i++)
	{
		if(!vis[i])pri[++cnt] = i,vis[i] = qpow(i,n);
		for(int j = 1; j <= cnt && pri[j]*i <= n; j++)
		{
			vis[pri[j]*i] = vis[pri[j]]*vis[i]%mod;
			if(pri[j]%i == 0)
			{
				break;
			}
		}
	}
}
int main()
{
	ll n;
	cin>>n;
	ll ans = 0;
	vis[1] = 1;
	oula(n);
	for(int i = 1; i <= n; i++)ans^=vis[i];
	cout<<ans;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
好的,这是一个比较简单的数学题,可以用矩阵快速幂求解。以下是 C++ 代码实现: ```c++ #include <iostream> #include <cstdio> #include <cstring> #include <algorithm> using namespace std; typedef long long LL; const int MAXN = 2; const int MOD = 1e9 + 7; struct Matrix { LL m[MAXN][MAXN]; Matrix() { memset(m, 0, sizeof(m)); } Matrix operator * (const Matrix& b) const { Matrix c; for (int i = 0; i < MAXN; ++i) { for (int j = 0; j < MAXN; ++j) { for (int k = 0; k < MAXN; ++k) { c.m[i][j] = (c.m[i][j] + m[i][k] * b.m[k][j]) % MOD; } } } return c; } } base, res; Matrix qpow(Matrix a, int b) { Matrix ans; for (int i = 0; i < MAXN; ++i) { ans.m[i][i] = 1; } while (b) { if (b & 1) { ans = ans * a; } a = a * a; b >>= 1; } return ans; } LL gcd(LL a, LL b) { return b == 0 ? a : gcd(b, a % b); } int main() { LL a, b, n; cin >> a >> b >> n; if (n == 1) { cout << a << endl; } else if (n == 2) { cout << b << endl; } else { base.m[0][0] = base.m[0][1] = base.m[1][0] = 1; res = qpow(base, n - 2); LL ans = gcd(a * res.m[0][0] % MOD + b * res.m[1][0] % MOD, b * res.m[1][0] % MOD + b * res.m[1][1] % MOD); cout << ans << endl; } return 0; } ``` 在这段代码中,我们定义了一个 `Matrix` 结构体,它表示一个 $2\times2$ 的矩阵。其中重载了 `*` 运算符,实现了矩阵乘法。 然后,我们定义了一个矩阵快速幂函数 `qpow`,用于求解矩阵的 $n$ 次方。 最后,在 `main` 函数中,我们通过快速幂求出矩阵 $base$ 的 $n-2$ 次方,然后根据题目要求求出 $\gcd(F_N, F_{N+1})$ 并输出即可。 需要注意的是,当 $n=1$ 或 $n=2$ 时,直接输出 $a$ 或 $b$ 即可。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值