单木的专栏

这题的思路还是比较神奇的。 虽然ai的取值范围比较大但与m取gcd之后取值范围就缩小到m的因子范围内。 然后在把0到m-1的数，按与m的最大公约数分类 如果存在ai（与m取最大公约数之后的ai），ai|ba_i|b，则与m的公约数是b的数都能取到。#include <bits/stdc++.h>#define pi acos(-1)#define pb push_back#define

给定一个n,∑ni=1n%i=n2−∑ni=1i∗⌊ni⌋\sum_{i=1}^n n\%i=n^2-\sum_{i=1}^n i*\left\lfloor\frac ni\right\rfloor由于⌊ni⌋\left\lfloor\frac ni\right\rfloor在一定范围内是保持不变的，所以我们接下来可以把时间复杂度降到n‾√\sqrt n，也可以直接跳 题目：51nod 1125

两种方法都不懂，留个纪念吧/* ********************************Author : danmuCreated Time : 2016年09月19日 星期一 17时29分04秒File Name : a.cpp******************************** */#include <algorithm>#inc

思路一开就想到多校第一场(hdu 5726 GCD）到后来一直没进展。 看了题解自己一开始的思路是对的，但后面的处理和本题不太一样。 固定右端点（都行），计算每个点到该端点的gcd（不同的gcd较少，相同的区间可以合并，时间和空间动能降下来）。然后在右端点固定的情况下，用树桩数组统计以每个点为起点的不同区间的不同gcd个数记到该起点中，对于每个询问的区间，求树桩数组对应区间的和，几位答案。代码h

题目cf 615D Multipliers

定理如果1<=b<a1<=b<a,a没有2或5的质因子，并且a与b互质，那么b/ab/a 的循环节位数恰好等于min(10e≡1(moda))min(10^e \equiv 1(mod a)),e是正整数。如果1<=b<a1<=b<a,a没有2或5的质因子，并且a与b互质，那么b/ab/a 的循环节位数必整除ϕ(a)\phi(a)。代码/* **************************

模板#include<stdio.h>#include<algorithm>#include<iostream>#include<string.h>#include<math.h>using namespace std;const int MAXN=50;int a[MAXN][MAXN];//增广矩阵int x[MAXN];//解集bool free_x[MAXN];//标记是否是不

求模素数P的原根的方法：对P-1素因子分解，即P-1=(P1^a1)(P2^a2)…..(Pk^ak)。，若恒有这里写图片描述成立，那么g就是P的原根(对于合数而言,只需要把p-1换成这里写图片描述即可)

相关题目poj2429#include <iostream> #include <stdlib.h> #include <string.h> #include <algorithm> #include <stdio.h> const int Times = 10; const int N = 5500; using namespace std; typedef long

int CRT(int a[],int m[],int n) { int M = 1; int ans = 0; for(int i=1; i<=n; i++) M *= m[i]; for(int i=1; i<=n; i++) { int x, y; int Mi = M

寻找大范围内的素数 二次筛法 poj2689

相关题目uvalive 20342.

1.

定理：设a>1a>1，m,n>0m,n>0，那么有gcd(am−1,an−1)=agcd(m,n)−1gcd(a^m-1,a^n-1)=a^{gcd(m,n)}-1。 题目：http://acm.hdu.edu.cn/showproblem.php?pid=2685定理：设a>ba>b，gcd(a,b)=1gcd(a,b)=1，那么gcd(am−bm,an−bn)=agcd(m,n)−bgcd

记π(n)\pi(n)为小于等于的素数个数，那么有limn→∞π(n)n/ln(n)=1\lim_{n\rightarrow\infty}\frac{\pi(n)}{n/\ln(n)}=1 。相关题目nefu 117

struct Mat{ int ma[35][35]; Mat unit(int n){ Mat res; memset(res.ma,0,sizeof res.ma); for(int i=0;i<n;++i) res.ma[i][i]=1; return res; }

代码#include&lt;cstdio&gt;#include&lt;cstdlib&gt;#include&lt;cstring&gt;struct Mat{ int ma[35][35]; Mat unit(int n){ Mat res; memset(res.ma,0,sizeof res.ma); for(int i=0;i&lt;n;++i)

思路∑0≤k1,k2,k3,⋯km≤n∏1≤j<mC(kj+1,kj)\sum_{0\leq k_1,k_2,k_3,\cdots k_m\leq n}\prod _{1\leq j<m}C(k_{j+1},k_j) =∑km=0nC(n,km)∑km−1=0kmC(km,km−1)⋯∑k1=0k2C(k2,k1)=\sum_{k_m=0}^nC(n,k_m)\sum_{k_{m-1}=0}^{k

#include <bits/stdc++.h>using namespace std;const int N = 80000+5;const int MOD = 1e9+7;int m;vector<int> vt[N], com[N];vector<pair<int,int> > prefix[N], suffix[N];bool mark[N];int g[N], pre_cnt[N

思路求连续范围内的最大公约数的和，为了降低时间复杂度，都是枚举最大公约数求和。/* ********************************Author : danmuCreated Time : 2016年08月01日 星期一 06时28分29秒File Name : hdu5780.cpp*************************

#include<cstdio>#include<cstdlib>#include<cmath>#include<algorithm>#define LL long longusing namespace std;bool devide(LL x){ for(LL i=2;i*i<=x;++i){ if(x%i==0){ x/=i;

代码#include<cstdio>#include<cstdlib>#include<cmath>const double constant_a = (1 + sqrt(5)) / 2;const double constant_b = (1 - sqrt(5)) / 2;const double constant_c = sqrt(5) / 5;int main(){ int

题目链接http://acm.hdu.edu.cn/showproblem.php?pid=2814题意给出a,b,n,c,求F(ab)F(ab)n−1%cF(a^{b})^{F(a^b)^{n-1}}\%c,F(x)表示斐波那契函数。思路先利用ab%p=aφ(p)+b%φ(p)a^{b}\%p = a^{\varphi (p) + b\%\varphi (p)}降幂，再用找斐波那契模一个数有循环节

题目链接http://acm.hdu.edu.cn/showproblem.php?pid=5768题意在给定范围内能够整除7并且不满足模pi等于ri的数的个数。思路代码#include<cstdio>#include<cstdlib>#define LL long longLL p[20],r[20];int n;LL sum;LL powMod(LL a,int b,int mod)

代码#include<cstdio>#include<cstdlib>#include<algorithm>#define LL long longusing namespace std;bool isprime[100010];int prime[100010],cnt;void init(){ memset(isprime,1,sizeof isprime); is

给定一个数n，log(n)+1log(n)+1就是n的阶乘的长度。 1. 对于小数据lenn!=logn!10+1=log110+log210⋯logn10+1len_{n!}=log_{10}^{n!}+1=log_{10}^1+log_{10}^2\cdots log_{10}^n+1 2. 对于大数据，用斯特林公式近似n!≈2πn−−−√(ne)nn!\approx \sqrt{2\pi

一个序列的逆序数的最大值是n×(n−1)/2n\times(n-1)/2。

性质： 当d与m互质时，同余方程满足消元性质。

LL solve(){ for(i=1;i<n;++i){ LL d,x,y; exgcd(m,ttm[i],d,x,y); if((tr[i]-r)%d) return -1; x*=(tr[i]-r)/d; x=x%(ttm[i]/d); r+=x*m; m=m/d*ttm[i

void eular(){ for(int i=1;i<=MAX;++i) phi[i]=i; for(int i=2;i<=MAX;i+=2) phi[i]/=2; for(int i=3;i<=MAX;i+=2) if(phi[i]==i) for(int j=i;j<=MAX;j+=i) phi[j

小学奥数，知乎上看到的。 原式为：(8111−9111)%172(8^{111}-9^{111})\%17^2 =>(8111+(17−8)111)%172(8^{111}+(17-8)^{111})\%17^2 二项式展开 =>(8111+∑111i=0×17111−i×(−8)i)%172(8^{111}+\sum^{111}_{i=0}\times17^{111-i}\times(-8

Simpson公式是一个数值积分公式，在计算一些多项式函数（三次或三次一下）的定积分时会得出精确值。但容易验证它对于)(xf=4x通常是不准确的。 ∫baf(x)dx=b−a6(f(a)+4f(a+b2)+f(b)) \int^b_af(x)dx=\frac{b-a}6(f(a)+4f(\frac{a+b}2)+f(b))

参见《具体数学第2版》244页。原谅我这个可耻的链接————————会整理一下在发出来相关性质：∑ni=0f(i)=f(n)∗f(n−1)\sum_{i=0}^nf(i)=f(n)*f(n-1);

不同为1,相同为零。 总结一些常用的性质吧，随时更新。 1. sis_i xor sjs_j = ai+1a_{i+1} xor…xor aja_j，s是a的前缀和。

Lucas定理：我们令n=sp+q , m=tp+r .（q ，r ≤p）那么：int Lucas (ll n , ll m , int p) { return m == 0 ? 1 : 1ll*comb (n%p , m%p , p) * Lucas (n/p , m/p , p) % p ;}//comb()函数中，因为q ， r < p , 所以

题目链接：http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3940题目大意：对于0到m区间的数模Ai,模A2...模An，later有q个询问，问前面的问题等于yi的数有多少个（1&lt;=i&lt;=q）最后在按照题意求一下最终结果。解题思路：做的时候最关键的部分没有想明白，没有找到好的处理方法，还是经验少吧，题解...

1.线性求所言逆元LL inv[MAX];inv[1]=1;for(int i=2;i<=MAX;++i) inv[i]=(mod-mod/i)*inv[mod%i]%mod;

6N+1法int prime[MAX],cnt;bool isPrime(int k){ if(k==2) return true; if(!(k&1)) return false; for(int i=3;i*i<=k;i+=2) if(k%i==0) return false; retur

定理：F(n)和f(n)是定义在非负整数集合上的两个函数，并且满足条件F(n)=∑d|nf(d)F(n)=\sum_{d|n}f(d)，那么我们得到结论f(n)=∑d|nμ(d)F(n/d)f(n)=\sum_{d|n}\mu(d)F(n/d)。 在上面的公式中有一个函数μ(d)\mu(d)，它的定义如下： （1）若d=1，那么μ(d)=1\mu(d)=1. （2）若d=p1p

a*x+b*y=d,d=gcd(a,b).模板代码：void exgcd(int a,int b,int& d,int& x,int& y){ if(b==0){ d=a; x=1; y=0; } else{ exgcd(b,a%b,d,y,x); y-=a/b*x; }}-----取自小白算法推导及证明：详见算导（第二版）528；具体数学