素数

素数

一些定理性质

• 唯一分解定理

• 威尔逊定理：

  $p$是质数，则有$(p-1)!\equiv -1(\mathrm{m}\mathrm{o}\mathrm{d}p)$

  逆定理同时成立，如果有$(p-1)!\equiv -1(\mathrm{m}\mathrm{o}\mathrm{d}p)$，则$p$一定是质数。

  有$((p-1)!+1)\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$可以通过$\sin (\pi ((n-1)!+1)/n)=0$来得到素数所在点函数。

• 费马小定理：

  假设$p$是一个质数，$a$不是$p$的倍数，则有$a^{p-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$，证明如下

  我们分别取$i\ast a(\mathrm{m}\mathrm{o}\mathrm{d}p)$的值，$foriinrange(1,p)$，由于$i$与$p$互质，$a$与$p$互质，所以，这些取模后的值是一个$(1n-1)$的排列，那么另这$p-1$个数相乘，得到的就是$a^{p-1}\ast (p-1)!(\mathrm{m}\mathrm{o}\mathrm{d}p)=(p-1)!$，因为$(p-1)!(\mathrm{m}\mathrm{o}\mathrm{d}p)!=0$，所以两边同时约去得到$a^{p-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$。

Miller_Rabin素数测试

利用了$a^p\equiv a(\mathrm{m}\mathrm{o}\mathrm{d}p)$，这一性质来对素数进行判定，但是具体的证明我也没仔细看，就是测试一个数字之后就能有$1/4$的概率证明其不是质数，然后通过多测试几个数，就可以几乎精确的得到这个数是不是质数。

#include <bits/stdc++.h>
#define mp make_pair
#define pb push_back

using namespace std;

typedef pair<int, int> pii;
typedef long long ll;
typedef unsigned long long ull;

const double eps = 1e-7;
const double pi = acos(-1.0);
const int inf = 0x3f3f3f3f;

inline ll read() {
    ll f = 1, x = 0;
    char c = getchar();
    while(c < '0' || c > '9') {
        if(c == '-') f = -1;
        c = getchar();
    } 
    while(c >= '0' && c <= '9') {
        x = (x << 1) + (x << 3) + (c ^ 48);
        c = getchar();
    }
    return f * x;
}


ll qpow(ll a, ll n, ll mod) {
    ll ans = 1;
    while(n) {
        if(n & 1)   ans = (ans * a) % mod;
        n >>= 1;
        a = (a * a) % mod;
    }
    return ans;
}

bool miller_rabin(ll a) {
    if(a == 2)  return true;
    for(int i = 1; i <= 10; i++) {
        ll temp = rand() % (a - 2) + 2;
        if(qpow(temp, a, a) != temp)    return false;
    }
    return true;
}

int main() {
    // freopen("in.txt", "r", stdin);
    // freopen("out.txt", "w", stdout);
    // ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
    srand(time(0));
    ll a;
    while(scanf("%lld", &a) != EOF) {
        if(miller_rabin(a)) puts("prime");
        else    puts("not");
    }
    return 0;
}

欧拉公式与欧拉定理

欧拉公式$\varphi (n)=x$，$x$是小于$n$并且与$n$互质的数字的个数，有$\varphi (8)=4(1,3,5,7)$。

欧拉公式成立的几个定理：

• 如果$p$是质数，则有$\varphi (p)=p-1$，这是显然成立的。

• 如果$p$是质数，则$\varphi (p^i)=(p-1)\ast p^{i-1}$，证明：

  在$[1,p^i)$之间，与$p$互质的数有$p\ast t(t=1,2,3,\dots \dots p^{i-1}-1)$，所以$\varphi (p^i)=p^i-1-(p^{i-1}-1)=(p-1)\ast p^{i-1}$

• $\varphi (n)=n\ast (1-1/p_1)\ast (1-1/p_2)\dots \dots (1-1/p_{n-1})\ast (1-1/p_n)$，$p_i$是$n$的质因数，例如：$8=2^3，\varphi (8)=8\ast (1-1/2)=4$

• 如果$a,b$互质，则有$\varphi (a\ast b)=\varphi (a)\ast \varphi (b)$，证明：

  可用上述的第三个定理来说明，$a,b$中的质因数是互质的，所以，任何两个$p_i$是互质的，然后显然成立。

欧拉定理

$a,m$互质有$a^{\varphi (m)}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}m)$，当$m$时质数的时候，也就是$a^{m-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}m)$

欧拉筛(素数+欧拉函数)

最后给上素数筛和欧拉函数的线性求解法。

#include <bits/stdc++.h>
#define mp make_pair
#define pb push_back

using namespace std;

typedef pair<int, int> pii;
typedef long long ll;
typedef unsigned long long ull;

const double eps = 1e-7;
const double pi = acos(-1.0);
const int inf = 0x3f3f3f3f;

inline ll read() {
    ll f = 1, x = 0;
    char c = getchar();
    while(c < '0' || c > '9') {
        if(c == '-') f = -1;
        c = getchar();
    } 
    while(c >= '0' && c <= '9') {
        x = (x << 1) + (x << 3) + (c ^ 48);
        c = getchar();
    }
    return f * x;
}

const int N = 1e5 + 10;

int prime[N], phi[N], cnt = 0;
bool st[N];

void init() {
    phi[1] = 1;
    for(int i = 2; i < N; i++) {
        if(!st[i]) {
            // cout << i << endl;
            prime[cnt++] = i;
            phi[i] = i - 1;
        }
        for(int j = 0; j < cnt && i * prime[j] < N; j++) {
            st[i * prime[j]] = 1;
            if(i % prime[j] == 0) {
                phi[i * prime[j]] = phi[i] * prime[j];
                break;
            }
            phi[i * prime[j]] = phi[i] * (prime[j] - 1);
        }
    }
}


int main() {
    // freopen("in.txt", "r", stdin);
    // freopen("out.txt", "w", stdout);
    // ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
    init();
    return 0;
}