小a与黄金街道欧拉数与快速幂

小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日。