Pentium.Labs

System全家桶：https://zhuanlan.zhihu.com/c_1238468913098731520

导航

统计

11 2014 档案

poj1811 数论
摘要：题意：判断一个数是否是质数+分解质因数sol：模板题分解质因数用xudyh模板，注意factor返回的是无序的，factorG返回是从小到大的顺序（包括了1）判断质数用kuangbin随机化模板 1 #include 2 #include 3 #include 4 #include ... 阅读全文

posted @ 2014-11-22 21:36 Pentium.Labs 阅读(234) 评论(0) 推荐(0) 编辑

POJ3735 矩阵
摘要：题意：有n只猫咪，开始时每只猫咪有花生0颗，现有一组操作，由下面三个中的k个操作组成： 1. g i 给i只猫咪一颗花生米 2. e i 让第i只猫咪吃掉它拥有的所有花生米 3. s i j 将猫咪i与猫咪j的拥有的花生米交换 现将上述一组操作做m次后，问每只猫咪有多少颗花生？sol: 可参考Mat... 阅读全文

posted @ 2014-11-18 20:55 Pentium.Labs 阅读(146) 评论(0) 推荐(0) 编辑

xudyh的gcd模板
摘要：hdu 5019 1 #include 2 #include 3 #include 4 #include 5 #include 6 #include 7 #include 8 #include 9 #include 10 #include 11 #include... 阅读全文

posted @ 2014-11-18 19:24 Pentium.Labs 阅读(692) 评论(0) 推荐(0) 编辑

[转]vi command summary
摘要：1 The following tables contain all the basic vi commands. 2 3 *Starting vi* 4 5 Command Description 6 vi file ... 阅读全文

posted @ 2014-11-18 13:11 Pentium.Labs 阅读(203) 评论(0) 推荐(0) 编辑

Uva11538 排列组合水题
摘要：画个图就很容易推出公式：设mn=min(m,n)，mx=max(m,n)对角线上：横向：m*C(n,2)纵向：n*C(m,2)因为所有的C函数都是只拿了两个，所以可以优化下。不过不优化也过了= = 1 #include 2 using namespace std; 3 #define LL lon... 阅读全文

posted @ 2014-11-17 21:42 Pentium.Labs 阅读(136) 评论(0) 推荐(0) 编辑

UVA11426 欧拉函数
摘要：大白书P125 1 #include 2 #include 3 using namespace std; 4 #define MMX 4000010 5 #define LL long long 6 int phi[MMX],f[MMX]; 7 LL S[MMX]; 8 9 void calc... 阅读全文

posted @ 2014-11-17 18:55 Pentium.Labs 阅读(172) 评论(0) 推荐(0) 编辑

关于gcd的几个问题
摘要：这两天刷了几个关于gcd的很类似的问题，总结一下：BZOJ2818 1<=x<=n，1<=y<=n，求满足gcd(x,y)=质数的个数BZOJ2190 1<=x<=n，1<=y<=n，求满足gcd(x,y)=1（x、y互质）的个数BZOJ2301 a<=x<=b，c<=x<=d，求... 阅读全文

posted @ 2014-11-17 16:34 Pentium.Labs 阅读(156) 评论(0) 推荐(0) 编辑

POI2007_zap 莫比乌斯反演
摘要：题意：http://hzwer.com/4205.html同hdu1695 1 #include 2 #include 3 #include 4 #include 5 using namespace std; 6 #define LL long long 7 #define MMX 5001... 阅读全文

posted @ 2014-11-16 17:25 Pentium.Labs 阅读(242) 评论(0) 推荐(0) 编辑

BZOJ2005 莫比乌斯反演
摘要：题意：http://www.lydsy.com/JudgeOnline/problem.php?id=2005实际上把这些被挡住的点的坐标和能量值列举出来可以发现有个公式：“对于坐标系第一象限任意的整点(即横纵坐标均为整数的点)p(n,m)，其与原点o(0,0)的连线上除过原点整点的个数为gcd(n... 阅读全文

posted @ 2014-11-16 15:37 Pentium.Labs 阅读(563) 评论(0) 推荐(0) 编辑

BZOJ2301 莫比乌斯反演
摘要：题意：a m) swap(n,m); for(int i = 1, la = 0; i 2 #include 3 #include 4 #include 5 using namespace std; 6 #define LL long long 7 #define MMX 5010... 阅读全文

posted @ 2014-11-15 15:57 Pentium.Labs 阅读(203) 评论(0) 推荐(0) 编辑

hdu1695 莫比乌斯反演
摘要：莫比乌斯反演：可参考论文：《POI XIV Stage.1 《Queries》解题报告By Kwc-Oliver》求莫比乌斯函数mu[i]：(kuangbin模板)http://www.cnblogs.com/kuangbin/archive/2013/08/21/3273440.htmlvoid ... 阅读全文

posted @ 2014-11-15 10:04 Pentium.Labs 阅读(714) 评论(0) 推荐(0) 编辑

BZOJ2190 欧拉函数
摘要：题意： http://www.lydsy.com/JudgeOnline/problem.php?id=2190这个矩阵是对称的，因此只要算出上三角的一半就行了一开始不太好看，不妨让坐标从0开始:eg：对于N=6的情况，令起始点坐标为(0,0)。那么所有能看见的点就是：1,1 1,2 1,3 ... 阅读全文

posted @ 2014-11-14 18:15 Pentium.Labs 阅读(377) 评论(0) 推荐(0) 编辑

莫比乌斯反演
摘要：定理：其中的是莫比乌斯函数。若i可以被除1之外的完全平方数整除，u(i)=0否则，设i的质因数的个数为k，则u(i)=(-1)^k一个简单的应用：求欧拉函数Reference：http://www.isnowfy.com/mobius-inversion/ 阅读全文

posted @ 2014-11-13 17:50 Pentium.Labs 阅读(173) 评论(0) 推荐(0) 编辑

BZOJ2818 欧拉函数
摘要：题意：求1--n中满足gcd(x,y)的值为质数的数对(x,y)的数目 （ (x,y)和(y,x)算两个 ）sol：设p[i]是一个质数，那么以下两个命题是等价的：1.gcd(x,y)=12.gcd(x*p[i],y*p[i])=p[i]eg:gcd(36,25)=1，gcd(36*7,25*7)... 阅读全文

posted @ 2014-11-13 15:28 Pentium.Labs 阅读(739) 评论(0) 推荐(0) 编辑

BZOJ3172 后缀数组
摘要：题意：求出一篇文章中每个单词的出现次数 对样例的解释：原文是这样的：aaaaaa注意每个单词后都会换行所以a出现次数为6，aa为3 (aa中一次，aaa中两次)，aaa为1标准解法好像是AC自动机or后缀自动机，还有人用KMP暴力过的= =用后缀数组做的。原来没刷过这种类型，顺便复习一下~Refer... 阅读全文

posted @ 2014-11-05 19:04 Pentium.Labs 阅读(245) 评论(0) 推荐(0) 编辑

BZOJ2186 欧拉函数
摘要：欧拉函数：一般记作φ(n)，表示1-n中与n互质的数的数量。欧拉函数是积性函数，即φ(m*n)=φ(m)*φ(n) //这条定理基友面试时还遇到了= =欧拉函数的值φ(n)=n*(1-p[1])*(1-p[2])*...*(1-p[n]) //p[i]是小于等于n的所有素数若n是m的倍数，则... 阅读全文

posted @ 2014-11-03 23:51 Pentium.Labs 阅读(333) 评论(0) 推荐(0) 编辑

poj 1845 数论综合
摘要：题意：求A^B的所有因数的和 mod 9901sol：一开始毫无思路，因为很多定理都不知道-_-||1. 整数的唯一分解定理：任意正整数都有且只有一种方式写出其素因子的乘积表达式。A=(p1^k1)*(p2^k2)*(p3^k3)*....*(pn^kn) 其中pi均为素数2. 约数和公式：对于已经... 阅读全文

posted @ 2014-11-03 17:38 Pentium.Labs 阅读(177) 评论(0) 推荐(0) 编辑

hdu 1576 求逆元
摘要：题意：给出n=A mod 9973和B，求(A/B) mod 9973昨天用扩展欧几里得做过这题，其实用逆元也可以做。逆元的定义：例如a*b≡1 (mod m)，则b就是a关于m的逆元。求逆元方法也很简单，用扩展欧几里得解这个方程即可。逆元性质：若a是b的逆元，则(x/a)mod p=(x*b)mo... 阅读全文

posted @ 2014-11-03 13:16 Pentium.Labs 阅读(262) 评论(0) 推荐(0) 编辑

poj 2891 扩展欧几里得迭代解同余方程组
摘要：Reference: http://www.cnblogs.com/ka200812/archive/2011/09/02/2164404.html之前说过中国剩余定理传统解法的条件是m[i]两两互质，所以这题就不能用传统解法了= =其实还有种方法：先来看只有两个式子的方程组：c≡b1 (mod a... 阅读全文

posted @ 2014-11-02 20:44 Pentium.Labs 阅读(326) 评论(0) 推荐(0) 编辑

hdu 3068 最长回文子串 马拉车模板
摘要：前几天用后缀数组写过一次这题，毫无疑问很感人的TLE了-_-||今天偶然发现了马拉车模板，O(N)时间就搞定reference：http://acm.uestc.edu.cn/bbs/read.php?tid=3258 1 #include 2 #include 3 #include 4 us... 阅读全文

posted @ 2014-11-02 20:31 Pentium.Labs 阅读(220) 评论(0) 推荐(0) 编辑

poj 2142 扩展欧几里得解ax+by=c
摘要：原题实际上就是求方程a*x+b*y=d的一个特解，要求这个特解满足|x|+|y|最小套模式+一点YY就行了总结一下这类问题的解法：对于方程ax+by=c设tm=gcd(a,b)先用扩展欧几里得求出方程ax+by=tm的解x0、y0然后有a*x0+b*y0=tm令x1=x0*(c/tm)，y1=y0*... 阅读全文

posted @ 2014-11-02 16:03 Pentium.Labs 阅读(316) 评论(0) 推荐(0) 编辑

poj 1061 扩展欧几里得解同余方程（求最小非负整数解）
摘要：题目可以转化成求关于t的同余方程的最小非负数解：x+m*t≡y+n*t (mod L)该方程又可以转化成：k*L+(n-m)*t=x-y利用扩展欧几里得可以解决这个问题：eg：对于方程ax+by=c设tm=gcd(a,b)若c%tm!=0，则该方程无整数解。否则，列出方程：a*x0+b*y0=tm易... 阅读全文

posted @ 2014-11-01 18:33 Pentium.Labs 阅读(557) 评论(0) 推荐(0) 编辑

同余、中国剩余定理
摘要：同余：a≡b (mod m)，表示a % m==b % m同余式的运算法则：--------------------------------------------------------------------------------------------------------------中国... 阅读全文

posted @ 2014-11-01 16:12 Pentium.Labs 阅读(1108) 评论(0) 推荐(0) 编辑

poj 1006 中国剩余定理解同余方程
摘要：其实画个图就明白了，该问题就是求同余方程组的解：n+d≡p (mod 23)n+d≡e (mod 28)n+d≡i (mod 33) 1 #include "iostream" 2 using namespace std; 3 int a[5],m[5]; 4 int p,e,i,d,ans; 5 ... 阅读全文

posted @ 2014-11-01 11:52 Pentium.Labs 阅读(307) 评论(0) 推荐(0) 编辑



Pentium.Lab Since 1998

点击右上角即可分享
微信分享提示