Rainfoo的博客

const ll mod = 1e9 + 7;const int N = _____;//输入矩阵大小Nstruct Mat{ ll m[N+1][N+1];};Mat a,e;Mat Mul(Mat x,Mat y){ Mat ans; for(int i=1;i<=N;++i){ for(int j=1;j<=N;++j){ ans.m[i][j]=0; for(int k=1;k<=N;

ine[1]=1; for(int i=2;i<mod;i++){ int a=mod/i,b=mod%i; ine[i]=(ine[b]*(-a)%mod+mod)%mod; }

初步上市#include<iostream>using namespace std;int main(){ cout<<"Hello world"<<endl; return 0;}

ll eular(ll n){ ll ans = n; for(int i=2; i*i <= n; ++i) { if(n%i == 0) { ans = ans/i*(i-1); while(n%i == 0) n/=i; } ...

欧拉定理费马小定理推论

矩阵快速幂就是把快速幂的乘法变成矩阵乘法。应用：求斐波那契数取模（大数）斐波那契数列递推公式（这里取从第二项开始）：f(1)=1，f(2)=2，f(n)=f(n-1)+f(n-2)(n>=3)用矩阵表示为：进一步，可以得出直接推导公式：...

模板如下:#include<bits/stdc++.h>using namespace std;#define ll long longconst int maxn=110;const int mod=1e9+7;#define mod(x) ((x)%mod)int n;struct mat{ int m[maxn][maxn];}unit;mat oper...

运算规则：模运算与基本四则运算有些相似，但是除法例外。其规则如下：1.(a + b) % p = (a % p + b % p) % p （1）2.(a - b) % p = (a % p - b % p + p) % p （2）3.(a * b) % p = (a % p * b % p) % p （3）4. a ^ b % p = ((a % p)^b) % p （4）结合律：(...

#include<bits/stdc++.h>using namespace std;bool book[100000001];//用埃氏筛法生成素数表void prime(int b) { //初始化，默认全部都是质数 memset(book, true, sizeof(book)); book[1]=false;//1不是质数 int n=...

介绍：如果p是一个质数，而整数a不是p的倍数，则有a^（p-1）≡1（mod p）即a^(p-1)%p=1%p(p是素数，a不是p的倍数)。用法：这样的话 a*a^(p-2)=1则a^(p-2)是a的逆元，这样我们要求a的逆元的话，就用快速幂quick_power(a,p-2)即可快速就出来...

#include<bits/stdc++.h>using namespace std;const int maxn=1e6+5;int prime[maxn],phi[maxn],cnt=0;bool vis[maxn];void init(){ phi[1]=1; for(int i=2;i<maxn;i++){ if(!vis[i]){ prime[c...