首先要知道什么是逆元
逆元定义:
对于ab≡1(mod m),b是a在模m下a的逆元。
首先对于同余符号≡,如 a≡b%m 表示的意思是a%m=b%m,a、b模m后余数相同。
逆元可以看作一个数的倒数 c * b≡1(mod m)(b为c的逆元)
可得(a/c)%m=(ab)%m
推导过程:(a/c)%m=(a/c)*1%m=(a/c) * c * b%m=(a * b)%m
由费马小定理:a、m是互质的正整数,有a^(m-1)≡1(mod m)
可得:a的逆元为:a^(m-2)
推导过程:a^(m-1)≡1(mod m)=a * a^(m-2)≡1(mod m),由逆元定义可得 a逆元=a^(m-2)
其实就是转变为通过快速幂找到逆元。
根据逆元的定义,gcd(i,p)==1才有逆元,而i的值是1~p-1,所以对于非质数来说,一定存在i与p不互质,也就是不存在逆元,直接AKCniubi,对于质数来说,可发现规律:其逆元之和为前 p-1 个数的和,然后直接用等差数列公式;
#include<iostream>
#include<cstdio>
#include<iomanip>
#include<cstdlib>
#include <algorithm>
#include<string.h>
#include<queue>
#include<math.h>
#include<set>
#define ll long long
using namespace std;
bool isPrime(ll n) {
if (n <= 3) {
return n > 1;
}
for (int i = 2; i <= sqrt(n); i++) {
if(n % i == 0) {
return false;
}
}
return true;
}
int main()
{
ll p;
while( cin >> p )
{
if(!isPrime(p))cout << "AKCniubi" << endl ;
else cout << (p-1)*p/2 << endl ; //等差数列通项公式
}
return 0;
}