费马小定理欧拉定理逆元

零衣贰

已于 2023-01-12 23:06:55 修改

阅读量192

点赞数 1

分类专栏：学习笔记文章标签：算法 c++ 开发语言

于 2022-11-17 20:02:19 首次发布

本文链接：https://blog.csdn.net/weixin_73113801/article/details/127909575

版权

学习笔记专栏收录该内容

45 篇文章 1 订阅

订阅专栏

费马小定理：
若 $p$ 为素数且 $a, p$ 互素，则 $a^{p-1} \equiv 1 \pmod p$
欧拉定理：
欧拉函数 $\varphi (m)$ 表示 $1$ 到 $m$ 中与 $m$ 互素的个数
若 $a, m$ 互素,则 $a^{\varphi(m)} \equiv 1\pmod m$

在这里插入图片描述

例题：汉堡与三明治

遇到计算 $\frac{1}{2} \pmod {10007}$ 等价于 $1*{(2的逆元)}$
在这里插入图片描述
计算时乘以 $5004$ 就相当于乘以 $\frac{1}{2}$

#include<bits/stdc++.h>
using namespace std;
const int mod=10007;
const int MAXN=2005;
int f[MAXN][MAXN];
bool mem[MAXN][MAXN];
int q(int n,int a,int b)
{
	if(mem[a][b]) return f[a][b];
	mem[a][b]=1;
	if(n==1)
	{
		if(a==0) return 0;
		if(b==0) return 1;
		return f[a][b]=5004;
	}
	else
	{
		if(a==0) return 0;
		if(b==0) return 1;
		return f[a][b]= ( 5004*q(n-1,a-1,b)%mod + 5004*q(n-1,a,b-1)%mod ) %mod;
	}
}
int n,a,b;
int main()
{
	cin>>n>>a>>b;
	cout<<q(n,a,b)%mod;
	return 0;
}

例题：无尽的循环

首先，根据题意可推出以下方程：
$a+cx\equiv b\space \pmod {2^k}$ , 若 $x$ 有解则循环次数有限，无解则无限
等价于
$cx\equiv b-a\pmod{2^k}$

在这里插入图片描述

最后使用快速幂即可

#include<bits/stdc++.h>
#define int long long
using namespace std;
int a,b,c,k,mod;
int get_pow(int x,int p)
{
	int res=1,base=x;
	while(p>0)
	{
		if(p & 1) res=(res*base)%mod;
		base=base*base%mod;
		p>>=1;
	}
	return res;
}
void print(int x)
{
	while(x<0)
		x+=mod;
	cout<<x%mod;
}
signed main()
{
	cin>>a>>b>>c>>k;
	mod=1LL<<k;	
	int p = (1<<(k-1)) - 1;	int t = b-a;
	if(c%2!=0)
		print( (t*get_pow(c,p))%mod );
	else
	{
		while(c%2==0)
		{
			if(t%2==0) c/=2, t/=2, mod/=2;
			else {puts("Inf"); return 0;}
		}
		int p = (mod>>1) - 1;
		print( (t*get_pow(c,p))%mod );
	}
	return 0;
}

例题：哈希碰撞

朴素的解法是枚举每一位的字母，时间复杂度为 $O(26^n)$ ，此题最高为 $26^{10}$
在这里插入图片描述

考虑折半搜索，复杂度降为 $26^5+26^5<<26^{10}$

#include<bits/stdc++.h>
#define int long long
using namespace std;
int m,n,t,mod,inv;
int get_pow(int a,int b)
{
	int base=a, res=1;
	while(b>0)
	{
		if(b&1) res*=base, res%=mod;
		base*=base, base%=mod;
		b>>=1;
	}
	return res;
}

int cnt[1<<25];

void f1(int x,int s)
{
	if(x==n/2) {cnt[s]++; return ;}
	for(int i=1;i<=26;i++)
		f1(x+1,((33*s)^i)%mod);
}

int f2(int x,int target)
{
	if(x==n/2)	return cnt[target];
	int res=0;
	for(int i=1;i<=26;i++)
	{
		int new_target=( (target^i)*inv ) % mod;
		res+=f2(x-1,new_target);
	}
	return res;	
}
signed main()
{
	cin>>m>>t>>n; mod=(1<<m);
	inv=get_pow(33,(1<<(m-1))-1);
	f1(0,0);
	cout<<f2(n,t);
	return 0;
}