whai的专栏

njoj 1251 zlly长了一张包子脸题意：zlly长了一张包子脸，他特别喜欢吃糖果，现在他手头有若干种糖果，每种糖果有个口味值，每种糖果有无数多个。然后娄童鞋也非常喜欢吃糖果，他的口味特别广泛，他喜欢各种各样的口味值，他要求zlly用现有的口味值拼出新的口味值。现在，娄童鞋想知道他不能吃到的口味值最大是多少？你能帮他的忙吗？举个例子，现在zlly手头有3,6,10三种糖果，他

1. 对于两个数A，B，如果A，B互质，那么Ax+By (x>=0 && y>=0) 最大的不能表示的数为AB-A-B，且不能表示数的个数为：(A-1)(B-1)/22. 对于3个数：定理一：设a,b,c为正整数，(a,b,c)=1，x,y,z为非负整数，ax+by+cz所不能表出的最大整数为M，那么当 c > ab/(a,b)^2 - a/(a,b) - b/(a,b) 时，M

bzoj 2186 [Sdoi2008]沙拉公主的困惑题意：大富翁国因为通货膨胀，以及假钞泛滥，政府决定推出一项新的政策：现有钞票编号范围为1到N的阶乘，但是，政府只发行编号与M!互质的钞票。房地产第一大户沙拉公主决定预测一下大富翁国现在所有真钞票的数量。现在，请你帮助沙拉公主解决这个问题，由于可能张数非常大，你只需计算出对R取模后的答案即可。R是一个质数。限制：数据组数T

欧拉函数定义：欧拉函数phi(n)，表示小于或等于n的数中与n互质的数的数目。欧拉函数的性质：1. phi(1)=12. 若n是素数p的k次幂：phi(n)=p^k-p^(k-1)=(p-1)p^(k-1)3. 若m,n互质，phi(mn)=phi(m)*phi(n)欧拉函数的递推式：令p为n的最小质因数若p^2|n，则phi(n)=phi(n/

codeforces 526 c Om Nom and Candies题意：给出一个容量为c背包，两种物品的价值hr,hb和容量cr,cb，求装满背包能装的最大价值。限制：1 思路：剩余系设第一种物品有x个，第二种物品有y个，则有：x*hr + y*hb 则有：(c-y*hb) % hr = t不难看出对于同一个剩余系t，最优值在两端。所以求

hdu 5201 The Monkey King 母函数 泰勒展开题意：有n个苹果，m个人，要求分给第一个人最多，其他人随意，求有多少种分法。最后结果模1000000007。限制：1 思路：母函数，泰勒展开枚举第一个人分到的苹果，设为u，剩下的苹果为n-u个，分成m-1份，则有：生成函数为：G(x)=(1+x+x^2+...+x^(u-1))^(m-

bzoj 2299 [HAOI2011]向量 裴蜀定理题意：给你一对数a,b，你可以任意使用(a,b), (a,-b), (-a,b), (-a,-b), (b,a), (b,-a), (-b,a), (-b,-a)这些向量，问你能不能拼出另一个向量(x,y)。限制：-2*1e9 思路：题目的操作可以化为：1. x +- 2a;2. y +- 2a;3.

裴蜀定理(1)a,b互质的充要条件是：存在整数x,y使ax+by=1(2)设a1,a2,a3......an为n个整数，d是它们的最大公约数，那么存在整数x1......xn使得x1*a1+x2*a2+...xn*an=d(3)如果a1...an互质（不一定是两两互质），那么存在整数x1......xn使得x1*a1+x2*a2+...xn*an=1(4)对任何整数a,b和它们的最

Falsyta in Tina TownTime Limit: 2000/1000 MS (Java/Others)    Memory Limit: 131072/131072 K (Java/Others)Total Submission(s): 548    Accepted Submission(s): 83Problem DescriptionFalsyt

离散对数,giant-step baby-step,拓展giant-step baby-step (poj 2417, hdu 2815)普通giant-step baby-step：poj 2417题意：B^L==N(MOD P), 给出P,B,N, 求最小的非负L。限制：2 思路：离散对数，用giant-step baby-step解

关于原根的一些知识点：定义1：设，，使得成立的最小的，称为对模的阶，记为。定理1：如果模有原根，那么它一共有个原根。定理2：若，，，则。定理3：如果为素数，那么素数一定存在原根，并且模的原根的个数为。定理4：设是正整数，是整数，若模的阶等于，则称为模的一个原根。   假设一个数对于模

E. Prefix Product Sequencetime limit per test1 secondmemory limit per test256 megabytesinputstandard inputoutputstandard outputConsider a sequence [a1, a2, ...

D. Moodular Arithmetictime limit per test1 secondmemory limit per test256 megabytesinputstandard inputoutputstandard outputAs behooves any intelligent schoolboy

C. Mike and Frogtime limit per test1 secondmemory limit per test256 megabytesinputstandard inputoutputstandard outputMike has a frog and a flower. His frog is n

1013: [JSOI2008]球形空间产生器sphereTime Limit: 1 Sec  Memory Limit: 162 MBSubmit: 3336  Solved: 1737[Submit][Status][Discuss]Description有一个球形空间产生器能够在n维空间中产生一个坚硬的球体。现在，你被困在了这个n维球体中，你只知道球面上n+1个点的坐

Widget FactoryTime Limit: 7000MS Memory Limit: 65536KTotal Submissions: 5082 Accepted: 1761DescriptionThe widget factory produces several different kinds of widge

Painter's ProblemTime Limit: 1000MS Memory Limit: 10000KTotal Submissions: 5079 Accepted: 2460DescriptionThere is a square wall which is made of n*n small square

EXTENDED LIGHTS OUTTime Limit: 1000MS Memory Limit: 10000KTotal Submissions: 7820 Accepted: 5090DescriptionIn an extended version of the game Lights Out, is a puz

hdu 1724 题意： 求椭圆面积, x^2/a^2+y^2/b^2=1, 给出a,b,l,r, 求椭圆被直线x=l, x=r所夹面积。思路： 1. 自适应simpson，精度不是很高，所以eps要设小一点。 2. romberg，精度很好，但效率不是很高。自适应simpson： /*hdu 1724 题意： 求椭圆面积, x^2/a^2+y^2/b^2=1, 给出a,b,

hdu 4498 自适应simpson 题意： 给出k1,k2,…,kn, a1,a2,…,an 和 b1,b2,…,bn 求函数： F(x)=min{100,min{ki*(x-ai)^2+bi | 0 < i <= n}} 在坐标上画出的曲线的长度。限制： 1 <= n <= 50; 0 <= ai,bi < 100; 0 < ki < 100思路： 先求出所有交点，然后排序，把函

五边形数定理 (hdu 4651 , hdu 4658)五边形数定理说得很详细的博文点击打开链接总结一下五边形数几个要点1.第n个五边形数的公式：五边形数 0  1  2  5  7  12  15对应下标 0  1  -1 2  -2  3  -32.五边形数定理：即：3.分拆数：4.

Miller_Rabin算法 (poj 3641)二次探测法：定理：如果p是素数，x是小于p的正整数，且(x^2)%p=1，那么要么x=1，要么x=p-1。过程：要测试N是否为素数，首先将N-1分解为d*2^s。在每次测试开始时，先随机选一个介于[1,N-1]的整数a，如果对于所有的r(0 题目：poj 3641题意：给出p,a, 判断p是不是对于a的伪素数。限

N次剩余 (hdu 3930)任务：给定N, a, p, 求出(x^N)%p=a 在模p意义下的所有解x。说明：令g为p的原根，因为p为素数，所以phi(p)=p-1。由原根的性质得：如果g为p的原根，则：g^i mod p != g^j mod p (p为素数), 其中i != j且i, j介於1至(p-1)之间所以，可以设g^y=x, g^t=a，则有：g^

指数循环节题集有个公式x >= Phi(C), A^x = A ^ (x%Phi(C) + Phi(C)) (mod C)不会证1.hdu 3221 Brute-force Algorithm题目链接点击打开链接题意：Function Find(int n,function func)  If n=1  For i=1 to a do func()

母函数 整数拆分 ( hdu 1028 , poj 1221 )介绍母函数的博文点击打开链接hdu1028 这道题是母函数的基本题目题目链接点击打开链接题意：整数拆分，例如4的情况4 = 4;4 = 3 + 1;4 = 2 + 2;4 = 2 + 1 + 1;4 = 1 + 1 + 1 + 1;限制：1 poj1221 这道题

poj 2689题意：给出(l,r), 求出区间[l,r]之间差最小，和差最大的相邻的素数对。限制：1 思路：筛法，对于2^32 只需要 2^16范围内的素数，就可以把后面的素数筛出。/*poj 2689 题意： 给出(l,r), 求出区间[l,r]之间差最小，和差最大的相邻的素数对。 限制： 1 <= l < r <= 21474836

poj 3696题意：给出一个数字L,求出最短的888...8能被L整除，输出最短的长度。限制：1 思路：设x为最小长度888...8=(10^x-1)/9*8由题意得：(10^x-1)/9*8 % L=0-> (10^x-1)*8 % (9L) = 0-> (10^x-1) % (9L/gcd(L,8)) = 0-> 10^x % (9L/gcd(

平方剩余 (poj 1808)题意：判断平方剩余，即判断(x^2)%p=a是否有解。限制：|a| 思路：用欧拉准则计算勒让德符号(用来判断平方剩余)代码，内含求平方剩余的解函数/*poj 1808 题意： 判断平方剩余，即判断(x^2)%p=a是否有解。 限制： |a| <= 1e9 && a % p !=0; 2 < p < 1

poj 1811题意：给出一个整数n，判断n是不是素数，如果不是素数，输出最小的质因子。限制；2 思路：miller_rabin算法判素数pollard_rho算法求质因子复杂度O(log(n))/*poj 1811 题意： 给出一个整数n，判断n是不是素数，如果不是素数，输出最小的质因子。 限制； 2 <= N < 2^54

poj 1707 伯努利数题目链接：点击打开链接题意：令自然数幂和为S(k,n)=1^k+2^k+...+n^kS(k,n)可以表述成：S(k,n)=1/M*(a[k+1]*n^(k+1)+a[k]*n^k+...+a[1]*n+a[0])给出k，求最小的正整数M，使得a[0]...a[k+1]都为整数。限制：1 题解：伯努利数S(k,n)=

自然数幂和与伯努利数 小结最近在看自然数幂和这类问题，觉得ACdreamer这两篇博客写得不错，可以看下。http://blog.csdn.net/acdreamers/article/details/38929067http://blog.csdn.net/acdreamers/article/details/39941387先对伯努利数的一些性质总结一下：1. 伯努

素数筛一个有用的性质：1~n的素数只需要用1~sqrt(n)的素数来筛，这个性质在局部筛的情况下很有用。两种素数筛的代码：/*素数筛 一个有用的性质：1~n的素数只需要用1~sqrt(n)的素数来筛，这个性质在局部筛的情况下很有用。 */#include #include #include #include using namespace std;

spoj 4491 莫比乌斯反演题意：给出a,b，求gcd(x,y)=prime的方案数，其中：1 限制：1 思路：先把问题拆成一个一个来考虑，然后问题就变成gcd(x,y)=k的方案数。设f(k)为gcd(x,y)=k的方案数，设F(k)为gcd(x,y)为k的倍数的方案数，显然F(k)=floor(a/k)*floor(b/k)。由莫比乌斯反演得：

中国剩余定理小结 (互质,非互质) (poj 1006,hdu 3579)先证明下中国剩余定理条件：x%m_1=a_1x%m_2=a_2...x%m_n=a_nm_1,m_2,...,m_n两两互质证明：设M=m_1*m_2*m_3*...*m_nM_i=M/m_i因为gcd(M_i,m_i)=1，所以M_ix+m_iy=1(t_i*M_

zoj 1864题目链接：点击打开链接题意：求自然数幂和。限制：0 思路：k不大，而且答案不取模，直接搞(n+1)^(k+1) - n^(k+1) = C(k+1,1)*n^k + C(k+1,2)*n^(k-1) + ... + C(k+1,k)*n + 1;(n+1)^(k+1) - 1 = ( (n+1)^(k+1) - n^(k+1) ) +

poj 2115 (解单变元模线性方程)http://poj.org/problem?id=2115题意：给出a,b,c,k,求x,使得(a+c*x)%(2^k)=b限制：0 思路：拓展欧几里得单变元模线性方程令 A=c;C=((b-a)%(2^k)+2^k)%(2^k);B=2^k则这道题就化为Ax%n=B对于Ax%B=C-> Ax

51Nod 1228 伯努利数题目链接：点击打开链接题意：S(k,n)=1^k+2^k+...+n^k求自然数幂和对1e9+7取模。限制：1思路：伯努利数S(k,n)=S(k,n)=1/(k+1) * ( C(k+1,k)*B[k]*(n+1)^1 + C(k+1,k-1)*B[k-1]*(n+1)^2 + ... + C(k+1,0)*B[0]*(

poj 3090 (欧拉函数,找规律)题目：给出一个n*n的点阵，求从(0,0)出发斜率不相等的直线有多少条。限制：1 <= n <= 1000思路：先定义sum[i]sum[i] = 0, if(i == 1)sum[i] = sum[i-1] + phi[i], if(i >= 2)ans = sum[n] * 2 + 3

hdu 5152 A Strange Problem题目链接点击打开链接'题意：给你一个长度为N的序列，序列为A1,A2,A3,...AN，然后有M个操作，每个操作为以下三种操作的其中一个：1. 输出操作。给你l,r，输出∑i=lrAi的值。2. 修改操作。给你x，把Ax修改为2Ax3. 加法操作。给你l,r,x，给Ai(l≤i≤r)加上x由于输出操作的结果可能