2014年10月29日
算法竞赛中的数论经典定理
摘要： 素数定理：记为小于等于的素数个数，那么有定理：设，，那么有定理：设，，那么定理：设，那么的值为（1）为素数，那么答案就是（2）有多个素因子，那么答案就是（3）只有一个素因子，那么答案就是该素因子定理：设为Fib数，那么有定理：给定两个互素的正整数和，那么它们最大不能组合的数为，不能组合的数的个数为定... 阅读全文
posted @ 2014-10-29 17:11 Kiss our dream 阅读(864) 评论(1) 推荐(0) 编辑
Baby Step Gaint Step
摘要： 给定同余式，求它在内的所有解，其中总是素数。分析：解本同余式的步骤如下 （1）求模的一个原根 （2）利用Baby Step Giant Step求出一个，使得，因为为素数，所以有唯一解。 （3）设，这样就有，其中，那么得到。 （4）求出所有的，可以知道一共有个解，我们求出所有的，然后排个序即... 阅读全文
posted @ 2014-10-29 15:38 Kiss our dream 阅读(226) 评论(0) 推荐(0) 编辑
2014年10月28日
素数分组 哥德巴赫猜想
摘要： 题目描述最少把1~n 分成多少组，可以使得每组的数的和为素数输入有多组数据第一行是一个数T，表示数据组数每组数据共1 行，为正整数n输出有T 行，每行为该情况的最少组数，无法分组时，输出-1样例输入12样例输出1哥德巴赫猜想裸题首先如果sum(n)是偶数，即两个素数之和，writeln(2)如果su... 阅读全文
posted @ 2014-10-28 16:57 Kiss our dream 阅读(250) 评论(0) 推荐(0) 编辑
FIB数列
摘要： 斐波那契级数除以N会出现循环，此周期称为皮萨诺周期。下面给出证明必然会出现循环这是基于下面事实：1．R(n+2)=F(n+2)modP=(F(n+1)+F(n))modP=(F(n+1)modp+F(n)modp)modp2．斐波那契数列的最大公约数定理：gcd(F(m),F(n))=F(gcd(m... 阅读全文
posted @ 2014-10-28 16:39 Kiss our dream 阅读(1566) 评论(0) 推荐(0) 编辑
构造矩阵来解数列递推问题
摘要： 转自某大牛。Fibonacci数列：F(0)=1 , F(1)=1 , F(n)=F(n-1)+F(n-2)我们以前快速求Fibonacci数列第n项的方法是 构造常系数矩阵（一） Fibonacci数列f[n]=f[n-1]+f[n-2],f[1]=f[2]=1的第n项快速求法（不考虑高精度）解法... 阅读全文
posted @ 2014-10-28 16:32 Kiss our dream 阅读(288) 评论(0) 推荐(0) 编辑
CERC 2013 Magical GCD
摘要： 题目大意如下：给定一个序列，每个序列有值xi,现给定t个数列,对于每个长n的数列，求一段[l,r]使 [r-l+1]*gcd(l,r)最大，gcd(l,r)指的是该连续区间的最大公约数。不难想到n^3,n^2logx,n^2的暴力吧n^3DP，n^2logx暴力枚举，n^2DP可以这样考虑，每次我对... 阅读全文
posted @ 2014-10-28 16:07 Kiss our dream 阅读(296) 评论(1) 推荐(0) 编辑
2014年10月27日
BZOJ 1096
摘要： const maxm=1e100; maxn=1000001; var f,x,p,c,sum,cost:array[0..maxn] of int64; q:array[0..maxn] of longint; n,i,h,t:longint; fun... 阅读全文
posted @ 2014-10-27 17:08 Kiss our dream 阅读(186) 评论(0) 推荐(0) 编辑
BZOJ 1062
摘要： program candy bzoj1062;const maxc=1000100; maxm=4010; maxn=2010;var n,len,m,i,p,t,l,r,c,d,q:longint; s:array[0..1,0..maxn,0..maxm] of lon... 阅读全文
posted @ 2014-10-27 17:07 Kiss our dream 阅读(300) 评论(0) 推荐(0) 编辑
BZOJ 1022
摘要： program bzoj1022;var t,n,i,ans,k,j,k1,k2:longint; bo:boolean;begin read(t); for k:=1 to t do begin read(n); ans:=0; k2:=0; ... 阅读全文
posted @ 2014-10-27 17:04 Kiss our dream 阅读(135) 评论(0) 推荐(0) 编辑
BZOJ 1023
摘要： program bzoj1023;uses math;const inf=1000000000; maxn=50005; maxm=20000005; type edge=record togo,next:longint;end; var n,m,cnt,ind,ans... 阅读全文
posted @ 2014-10-27 17:04 Kiss our dream 阅读(169) 评论(0) 推荐(0) 编辑