【欧拉函数】小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−1n−1出发走向11，两人的速度均为1m/s1m/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−时游戏结束
小a想知道在游戏结束时A+B的值
答案对109+7取模

---

【小结】

个人认为这个题的确不错，出题人的能力还是不错的。因为这个是考察了欧拉函数，而且是用了相应的等差数列的特点。

【题解】

题意不多说了，就算看不懂中文，其实看看图也理解了差不多了，就是两个人，一个从1开始出发，另一个是在n-1的地方出发。

然后两个人只要同时满足gcd(x,n) x是他所在的地方。

然后这个考察了欧拉函数的使用，为什么这么说呢，由于，gcd(x,n)明显就是求phi(n)

因为，如果某一个人的位置gcd(x,n)!=1，那么gcd(n-x,n)也同样是不等于1，因为其实就是，

如果x是n的一个因子，那么必定能线性表示。这个其实性质和更相减损法一样。

其实两个人走的地方都是一样，即需要乘法的地方都是一样，而且是满足等差中项的性质。

x+(n-x)=n，所以答案就是

A=A*K^(phi(n)*n/2)

B=B*K^(phi(n)*n/2)

ans=A+B

---

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=1e9+7;
ll qpow(ll a,ll b){
    ll ans=1;
    while(b){
        if(b&1) ans=ans*a%mod;
        a=a*a%mod;b>>=1;
    }
    return ans;
}
ll phi(ll x){
    ll tx=x,ans=1;
    for(ll i=2;i*i<=x;i++){
        int cnt=0;
        if(x%i==0){
            while(x%i==0){
                x/=i;
                cnt++;
            }
            ans=ans*(qpow(i,cnt-1)*(i-1)%mod)%mod;
        }
    }
    if(x!=1){
        ans=ans*(x-1)%mod;
    }
    return ans;
}
int main()
{
    ll n,A,B,k;
    scanf("%lld%lld%lld%lld",&n,&k,&A,&B);
    A=A*qpow(k,(n*(phi(n)*qpow(2,mod-2))%mod)%mod)%mod;
    B=B*qpow(k,(n*(phi(n)*qpow(2,mod-2))%mod)%mod)%mod;
    ll ans=(A+B)%mod;
    printf("%lld\n",ans);
    return 0;
}