链接:https://ac.nowcoder.com/acm/contest/317/D
来源:牛客网
小a与黄金街道
时间限制:C/C++ 1秒,其他语言2秒
空间限制:C/C++ 32768K,其他语言65536K
64bit IO Format: %lld题目描述
小a和小b来到了一条布满了黄金的街道上。它们想要带几块黄金回去,然而这里的城管担心他们拿走的太多,于是要求小a和小b通过做一个游戏来决定最后得到的黄金的数量。
游戏规则是这样的:
假设道路长度为nn米(左端点为00,右端点为nn),同时给出一个数kk(下面会提到kk的用法)
设小a初始时的黄金数量为AA,小b初始时的黄金数量为BB
小a从11出发走向n−1n−1,小b从n−1n−1出发走向11,两人的速度均为1m/s1m/s
假设某一时刻(必须为整数)小a的位置为xx,小b的位置为yy,若gcd(n,x)=1gcd(n,x)=1且gcd(n,y)=1gcd(n,y)=1,那么小a的黄金数量AA会变为A∗kx(kg)A∗kx(kg),小b的黄金数量BB会变为B∗ky(kg)B∗ky(kg)
当小a到达n−1n−1时游戏结束
小a想知道在游戏结束时A+BA+B的值
答案对109+7109+7取模输入描述:
一行四个整数n,k,A,Bn,k,A,B输出描述:
输出一个整数表示答案示例1
输入
4 2 1 1输出
32说明
初始时A=1,B=1A=1,B=1
第一个时刻如图所示,小a在11,小b在33,满足条件,此时A=1∗21=2,B=1∗23=8A=1∗21=2,B=1∗23=8
第二个时刻小a在22,小b在22,不满足条件
第三个时刻小a在33,小b在11,满足条件,此时A=2∗23=16,B=8∗21=16A=2∗23=16,B=8∗21=16
此时游戏结束A=2∗23=16,B=8∗21=16A=2∗23=16,B=8∗21=16
A+B=32A+B=32
示例2
输入
5 1 1 1输出
2备注:
保证3⩽n⩽108,1⩽A,B,k⩽1013
//若gcd (?,?) = 1,那么gcd (?,?−?)一定等于1
又因为前?个数的欧拉函数之和为
又因为指数很大,所以需要求2关于1e9+7的逆元
然后需要有一个扩展欧拉函数定理
所以:
This is the code
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<iostream>
#include<iomanip>
#include<list>
#include<map>
#include<queue>
#include<sstream>
#include<stack>
#include<string>
#include<set>
#include<vector>
using namespace std;
#define PI acos(-1.0)
#define EPS 1e-8
#define MOD 1e9+7
#define LL long long
#define ULL unsigned long long //1844674407370955161
#define INT_INF 0x7f7f7f7f //2139062143
#define LL_INF 0x7f7f7f7f7f7f7f7f //9187201950435737471
const int dr[]= {0, 0, -1, 1, -1, -1, 1, 1};
const int dc[]= {-1, 1, 0, 0, -1, 1, -1, 1};
// ios::sync_with_stdio(false);
// 那么cin, 就不能跟C的 scanf,sscanf, getchar, fgets之类的一起使用了。
const int maxn=1e6;
const LL mod=1e9+7;
vector<LL > prime;
vector<LL > num1;
vector<LL > num2;
bool not_prime[maxn+10];//初始为false
LL Euler(LL x)
{
LL res=x;
for(LL i=2; i<=(LL)sqrt(x); i++)
{
if(x%i==0)
{
res=res/i*(i-1);
while(x%i==0)
x/=i;
}
}
if(x>1)
res=res/x*(x-1);
return res;
}
LL power_mod(LL a,LL n)//快速幂求模运算
{
LL ans=1;
a%=mod;
while(n)
{
if(n&1)
ans=ans*a%mod;
n>>=1;
a=a*a%mod;
}
return ans%mod;
}
LL mu(LL a,LL b)
{
LL ret=0;
while(b)
{
if(b&1)
ret=(ret+a)%mod;
b>>=1;
a=(a+a)%mod;
}
return ret%mod;
}
int main()
{
LL n,k,a,b;
cin>>n>>k>>a>>b;
/*数据量在1e13时
LL Mod=Euler(mod);//1e9+6
LL inv2=500000004;
LL tem=((n%mod*Euler(n)%mod*inv2)%mod)%Mod+Mod;
*/
//long long满足精度时
LL tem=n*Euler(n)/2;
LL K=power_mod(k,tem);
LL ans=(mu(a,K)+mu(b,K))%mod;
cout<<ans<<endl;
return 0;
}