アルゴリズムの基礎 - 数学的知識 - オイラー関数、高速べき乗、拡張ユークリッド、中国剰余定理
オイラー関数
共素とは、2 つの数値の最大公約数が 1 のみであることを意味し、これがコードに反映されていますa和b互质,则b mod a = 1 mod a
(現時点ではよく理解できませんが、次のように理解できます: a と b の最大公約数は次のようになります) 1、つまり、1 を除数として使用し、b を約数として使用する場合、被除数 a については、余りが同じです。つまり、1/a の余りは b/a と同じです。 )b mod a = 1 mod a
オイラー関数の役割は、
互いに素である 1-n と n の数を見つけることです。
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <string.h>
#include <string>
#include <math.h>
#include <vector>
#include <queue>
#include <map>
#include <unordered_map>
using namespace std;
void get_eura(int x)
{
int res = x;
for (int i = 2; i <= x / i; ++ i)
{
if (x % i == 0)
{
//res = res * (1 - 1/i);或者res = res * (i - 1) / i;都不行,前者是浮点数1 后者会溢出
res = res / i * (i - 1);
while (x % i == 0)
{
x /= i;
}
}
}
if (x > 1) res = res / x * (x - 1);
cout << res << endl;
}
void solve()
{
int n;
cin >> n;
while (n -- )
{
int x;
cin >> x;
get_eura(x);
}
}
int32_t main()
{
ios::sync_with_stdio(0);
cin.tie(0);
int T = 1;
//cin >> T;
while (T --) solve();
return 0;
}
AcWing 874. オイラー関数を求めるふるい法
リニアシーブ+オイラー関数(微押し式付き)
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <string.h>
#include <string>
#include <math.h>
#include <vector>
#include <queue>
#include <map>
#include <unordered_map>
using namespace std;
const int N = 1e6 + 10;
int primes[N], st[N], eulers[N];
int cnt;
void get_eulers(int x)
{
eulers[1] = 1;
for (int i = 2; i <= x; ++ i)//只是在线性筛的过程中顺便求了一下每个数的欧拉函数
{
if (!st[i])//1-n的质数
{
primes[cnt++] = i;
eulers[i] = i - 1;
}
for (int j = 0; primes[j] <= x / i; ++ j)//1-n的合数//任何合数都含有质因数,4 = 1 * 2 * 1 * 2;
{
st[primes[j] * i] = 1;
if (i % primes[j] == 0)
{
eulers[i * primes[j]] = eulers[i] * primes[j];
break;//其实也相当于一个else
}
//eulers[i * primes[j]] = eulers[i] * primes[j] / primes[j] * (primes[j] - 1);
eulers[i * primes[j]] = eulers[i] * (primes[j] - 1);
}
}
}
void solve()
{
int n;
cin >> n;
get_eulers(n);
long long res = 0;
for (int i = 1; i <= n; ++ i) res += eulers[i];
cout << res;
}
int32_t main()
{
ios::sync_with_stdio(0);
cin.tie(0);
int T = 1;
//cin >> T;
while (T --) solve();
return 0;
}
高速パワー
1 2 4 8 ログの時間計算量が指数関数的に増加します
AcWing 875. ファストパワー
long long qmi(int a, int b, int p)
{
long long res = 1;
while (b)
{
if (b & 1)
{
res = res * a % p;
}
a = a * (long long)a % p;
b >>= 1;
}
return res;
}
AcWing 876. 高速電力反転要素
オイラー関数 => フェルマーの定理 => フェルマーの定理計算結果の高速なべき乗実現
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <string.h>
#include <string>
#include <math.h>
#include <vector>
#include <queue>
#include <map>
#include <unordered_map>
using namespace std;
long long qmi(int a, int b, int p)
{
long long res = 1;
while (b)
{
if (b & 1) res = res * a % p;
a = (long long)a * a % p;
b >>= 1;
}
return res;
}
void solve()
{
int n;
cin >> n;
while (n --)
{
int a, p;
cin >> a >> p;
if (a % p == 0) cout << "impossible" << endl;
else cout << qmi(a, p - 2, p) << endl;//a需要与m互质,否则a不存在乘法逆元
}
}
int32_t main()
{
ios::sync_with_stdio(0);
cin.tie(0);
int T = 1;
//cin >> T;
while (T --) solve();
return 0;
}
拡張ユークリッド (Pei Shu の定理)
AcWing 877. 拡張ユークリッド アルゴリズム
再帰の性質を理解する:
Peishu の定理と線形合同方程式の証明:
#include <cstdio>
#include <iostream>
using namespace std;
int exgcd(int a, int b, int &x, int &y)
{
if (b == 0)
{
x = 1, y = 0;
return a;
}
//d就是最大公约数,本题其实用不到
int d = exgcd(b, a % b, y, x);//本题的精髓
/*
只是为了方便更改x和y的值,如果用
d = exgcd(b, a % b, x, y);
最后就解得 新x = y 新y = x - a / b * y
那么代码就得这么写
int t = y;
y = x - a / b * y;
x = t;
显然比只要写一句 新y -= a / b * x; 麻烦
*/
y -= a / b * x;
return d;
}
void solve()
{
int n;
cin >> n;
while (n -- )
{
int a, b, x, y;
cin >> a >> b;
exgcd(a, b, x, y);
cout << x << " " << y << endl;
}
}
int32_t main()
{
ios::sync_with_stdio(0);
cin.tie(0);
int T = 1;
//cin >> T;
while (T --) solve();
return 0;
}
AcWing 878. 線形合同方程式
線形合同方程式は拡張ユークリッドの定理を使用して解きます。
この問題の導出過程は上記の通りです。
なぜ % m が必要なのでしょうか?
#include <cstdio>
#include <iostream>
using namespace std;
int exgcd(int a, int b, int &x, int &y)
{
if (b == 0)
{
x = 1, y = 0;
return a;
}
else//其实不用else,上面满足直接return了,上面不满足也会走到下面
{
int d = exgcd(b, a % b, y, x);
y -= a / b * x;
return d;
}
}
void solve()
{
int n;
cin >> n;
while (n -- )
{
int a, b, m, x, y;
cin >> a >> b >> m;
int d = exgcd(a, -m, x, y);
if (b % d != 0) cout << "impossible" << endl;
else cout << (long long)b / d * x % m << endl;
}
}
int32_t main()
{
ios::sync_with_stdio(0);
cin.tie(0);
int T = 1;
//cin >> T;
while (T --) solve();
return 0;
}