欧拉降幂快速幂用。。。。（前天打比赛真的是长见识了。。。）

emmmm，昨天第一次做牛客的比赛，就出来两个（真的是没救了。。）做第三个的时候看见的一个用快速幂的题，发现("▔□▔)肿么跟我学的那个快速幂长得这么不一样，今天看别人写的代码，顺便总结一下快速幂的用法。。

1.先上一个普通的取模的快速幂。。。

ll pow(ll a, ll b, ll mod)
{
	ll r = 1;
	a = a % mod;
	while (b > 0)
	{
		if (b & 1 == 1)r = (r*a) % mod;
		b >>= 1;
		a = (a*a) % mod;
	}
	return r;
}

再来一个昨天看的大神的代码。。
2.

#pragma warning(disable:4996)
#include <iostream>
#include <cstring>
#include <stdio.h>
#include <string>
#include <algorithm>
using namespace std;
typedef long long ll;
const ll mod = 1e9 + 7;
ll oula(ll n)
{
	ll s = n;
	for (int i = 2; i*i <= n; i++)
	{
		if (n%i == 0)
		{
			s = s / i * (i - 1);
			while (n%i == 0)n /= i;
		}
	}
	if (n > 1)s = s / n * (n - 1);
	return s;
}
ll kuickpow(ll x, ll y)
{
	ll ans = 1;
	while (y)
	{
		if (y & 1)ans = (ans*x) % mod;
		y >>= 1;
		x = (x*x) % mod;
	}
	return ans;
}
int main()
{
	ll n, k, a, b, sum;
	scanf("%lld%lld%lld%lld", &n, &k, &a, &b);
	sum = n * oula(n) / 2 % mod;
	printf("%lld\n", (a + b)*kuickpow(k, sum) % mod);
	return 0;

}

这里运用了一个公式是 a^x=a(x%φ(m)+ φ(m))(mod m)
证明的话可以在这个链接里面看看
https://blog.csdn.net/synapse7/article/details/19610361。。。
感觉自己还差得远。。。。。