----------数论---------
Werky_blog
这个作者很懒,什么都没留下…
展开
-
模板:快速幂
快速幂:ll quick_pow(ll x, ll n){ //快速幂 ll ans=1; while(n){ if(n&1){ ans=ans*x; } x=x*x; n>>=1; } ...原创 2017-11-05 19:20:10 · 201 阅读 · 0 评论 -
模板:数学2
中国剩余定理:原创 2017-10-22 00:43:19 · 177 阅读 · 0 评论 -
模板:数学1
素数打表: const int maxn=1001005; int p[maxn]; //第 i 个质数 bool pp[maxn]; //记录是不是质数 int get_prime(){ int pnum=0; for(int i=2;i<=maxn;i++){ if(!pp[i]) p[pnum++]=i;原创 2017-10-18 19:21:11 · 202 阅读 · 0 评论 -
模板:莫比乌斯反演
如果有: F(d)=∑i|df(i) 那么有: f(d)=∑i|dμ(i)F(di) 另一种形式是如果 F(d)=∑d|if(i) 那么 f(d)=∑d|iμ(id)F(i) const int maxn=1e6+10; ll prime[maxn],mob[maxn],vis[maxn],cnt; void M原创 2017-10-25 23:53:45 · 245 阅读 · 0 评论 -
模板:miller-rabin素数测试(哥德巴赫猜想)
含卡迈尔数:#include<bits/stdc++.h> #define ll unsigned long long int using namespace std; ll ModMul(ll a,ll b,ll n){//快速积取模 a*b%n ll ans=0; while(b){ if(b&1) ans=(ans+a)%...原创 2018-05-05 14:14:51 · 178 阅读 · 0 评论