欧拉降幂问题

题目连接

        D-幂运算_2023河南萌新联赛第（四）场：河南大学 (nowcoder.com)

        djz的数学题 (nowcoder.com)

 

欧拉降幂

当 b 过大时，快速幂的时间复杂度超过题目要求时，考虑欧拉降幂。 

 

 §（c）是 欧拉函数 里求小于他的跟他互质的数的个数

伪代码

long long oula(long long b)
{
	long long ans=b;
	for(int i=2;i*i<=b;i++)
	{
		if(b%i==0)
		{
			ans=ans/i*(i-1);
			while(b%i==0) b/=i;
		}
	}
	
	if(b>1) ans=ans/b*(b-1);
	return ans;
}

求出后套板子即可。

看不懂板子看下面代码该如何操作就行拉

题-F-数学题

 这题通过手磨，拆两次括号可以得到，

原式=1+x^1+x^2+x^3+x^4……x^(2*2^n-1).

也就是求首相为1，公比为x的前2^(n+1)项和。

这个题目的余数为 998244353 是个质数，可以求出他的 oula 值为 998244352

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int M=998244353;
ll ksm(ll a,ll x,ll mod)
{
	ll cnt=1;
	while(x)
	{
		if(x%2)	cnt=cnt*a%mod;
		a=a*a%mod;
		x/=2;
	}
	return cnt;
}

int main()
{
	ll n,x;
	cin>>n>>x;
	x%=M;
	if(x==1)	cout<<ksm(2,n+1,M);
	else
	{
		ll q=ksm(2,n+1,M-1);
		q=(ksm(x,q,M)-1+M)%M;
		ll p=ksm(x-1,M-2,M);
		cout<<(p*q%M+M)%M;
	}
	return 0;
}

题-D-幂运算

 题目的mod需要输入，我们先求出其oula值，在套公式

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll n,p;
long long oula(long long b)
{
	long long int ans=b;
	for(int i=2;i*i<=b;i++)
	{
		if(b%i==0)
		{
			ans=ans/i*(i-1);
			while(b%i==0) b/=i;
		}
	}
	
	if(b>1) ans=ans/b*(b-1);
	return ans;
}
ll ksm(ll a,ll x,ll mod)
{
	ll cnt=1;
	while(x)
	{
		if(x%2)	cnt=cnt*a%mod;
		a=a*a%mod;
		x/=2;
	}
	return cnt;
}

int main()
{
    cin>>n>>p;
    ll pp=oula(p);
    n=ksm(2,n,pp);
    cout<<ksm(2,n,p);
	return 0;
}