小a与黄金街道
题目描述:
题目连接:https://ac.nowcoder.com/acm/contest/317/D
小a和小b来到了一条布满了黄金的街道上。它们想要带几块黄金回去,然而这里的城管担心他们拿走的太多,于是要求小a和小b通过做一个游戏来决定最后得到的黄金的数量。
游戏规则是这样的:
假设道路长度为n米(左端点为0,右端点为n),同时给出一个数k(下面会提到k的用法)设小a初始时的黄金数量为A,小b初始时的黄金数量为B。小a从1出发走向n−1,小b从n−1出发走向1,两人的速度均为1m/s。假设某一时刻(必须为整数)小a的位置为x,小b的位置为y,若gcd(n,x)=1且gcd(n,y)=1,那么小a的黄金数量A会变为A∗kx(kg),小b的黄金数量B会变为B∗ky(kg),当小a到达n−1时游戏结束小a想知道在游戏结束时A+B的值,答案对10^9+7取模。
输入描述:
一行四个整数n,k,A,B
输出描述:
输出一个整数表示答案
备注:
保证3⩽n⩽ 10^8,1⩽A,B,k⩽ 10^13
分析:这道题考察了欧拉数和快速幂。
代码:
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const LL X=1e9+7;
LL quickpower(LL k,LL sum)//快速幂模板
{
int res=1;
while(sum)
{
if(sum&1)res=res*k%X;
k=k*k%X;
sum>>=1;
}
return res;
}
int main()
{
LL n,k,a,b;
cin>>n>>k>>a>>b;
LL sum=n,euler=n;
for(LL p=2;p*p<=n;p++)//欧拉数求法,euler=n*(1-1/a)*(1-1/b)···即,数n的欧拉数为:n乘所有的n的因子的倒数被一减去的差。欧拉数的含义:1-n所有与n互质的数的个数。一定为偶数,因为n/2与n不互质(n=2的情况特殊但也符合偶数)。
{
if(n%p==0)
{
euler=euler/p*(p-1);
while(n%p==0)n/=p;
}
}
if(n>1)euler=euler/n*(n-1);
sum=sum*euler/2;
cout<<(a+b)*quickpower(k,sum)%X<<endl;
return 0;
}
2月14日。