2015年07月_二分抄代码

原创 Hello world

来自HN.长沙市一中高一 (from 2014) 因为GB是无上的力量，所以就叫这个名字了小学和初中时候就接触过OI但一直没认真学，看着现在某礼某郡等同省竞赛强校的神犇们初中省队，初中金牌什么的，很后悔啊… 然而还是比较颓，我也没有同机房的神犇们聪明，马上就要滚粗辣参加了NOIP2014 当时学了一些算法但考得很差，与TBTB等同校神犇形成鲜明对比（可以去她的博客园看看）有些东西写

2015-07-01 22:20:12 1618 6

原创数学专题总结

数学总结在感受了picks的课后，发现OI中的数学是非常难的。虽然最近考试的题目并不难，但由于我天生推不出公式想不出数学方法智商捉急的基因导致各种考试想不到，主要原因还是做题做少了，方法知识不够健全，还是有必要总结一下的。求一个数对于mod的逆元：费马定理： P为素数，a为正整数，a与P互质，则a^(p-1) ≡ 1 (mod p) 所以一个数的逆元就是a^(p-2)应用：在取模的

2015-07-14 17:37:48 1085

原创 [DP]【最大全零矩阵】【2015.7.9TEST】E

E 0.9 seconds, 32 MB “ 于是乎，你至少证明了你智商比金天成高，也就说你证明了你不是低智儿童，不错不错。然而这次，我貌似也卡住了，你给我打下手吧。勇敢的少年啊快去创造奇迹！” ——-By Doctor Z 貌似 Z 博士正在解析 Zvangelion 初号机的一些问题，中间遇到了困难。 Zvangelion 初号机有一块 R*S 的电路模块被某种 UMA 感染了

2015-07-10 14:55:12 1075 1

原创 [NOIP2009][spfa]最优贸易

都7月份了，我连联赛题都想不出我真是可以滚粗了！这道题就是对1一个正向spfa记录下到i时能买到的最小价格xiao[i]，对n一个反向spfa，记录下到i时能卖出的最小价格da[i],然后枚举1到n号节点寻找最大的差值。#include<cstdio>#include<cstring>#include<vector>#define maxl 100001using namespace

2015-07-08 22:10:18 1168

原创 O（n）求素数，求欧拉函数，求莫比乌斯函数，求对mod的逆元，各种求

筛素数void shai(){ no[1]=true;no[0]=true; for(int i=2;i<=r;i++) { if(!no[i]) p[++p[0]]=i; int j=1,t=i*p[1]; while(j<=p[0] && t<=r) { no

2015-07-08 10:29:11 2172 5

原创 [BZOJ 3884][欧拉定理]上帝与集合的正确用法

看看我们机房某畸形写的题解：http://blog.csdn.net/sinat_27410769/article/details/46754209此题为popoQQQ神犇所出，在此orz#include<cstdio>#include<cstring>#include<cmath>#define maxl 10000001long long ans;long long f[maxl];b

2015-07-04 11:15:33 842 2

原创【BZOJ1951】[中国剩余定理][SDOI2010]古代猪文

求g的p次方%mod，根据费马小定理，g^sigma(C(n,d))(d|n)%mod=g^(sigma(C(n,d))(d|n)%(mod-1))%mod,然而mod-1不是质数，只能用把它拆成4个质因数，然后对4个模方程分别求解，先用lucas定理和费马小定里求出对4个质数取模的sigma的值(num[i])，注意，枚举因数d的时候枚举到sqrt(n)就可以了，同时加上C(N,I)和C（n，

2015-07-04 11:11:00 1421

问第k大的无平方因子数，二分求[1,x]中有多少无平方因子数，根据容斥原理， Q=x-x内有一个平方因子的数+x内有两个平方因子的数.. =x-x内(4的倍数个数+9的倍数的个数+25的倍数的个数)+x内(36(2*3*2*3)的倍数+100(2*5*2*5)的倍数)…，刚好莫比乌斯函数对于奇数个质数相乘的数是-，偶数个质数相乘的数是+，所以枚举i，让i的平方<n，[1,x]中有i的平方的数就是n/

2015-07-01 21:47:05 498

原创 [递推]Bzoj3930 选数[CQOI2015]

在l到r的数字中，含有因数为i*k的个数为r/(i*k)-(l)/(i*k),方案数就为t^n种，但gcd=i*k，则选择相同的n个数和f[i*k]（2<=j<=m）都是不行的要减去，所以从大到小枚举i，最后输出f[1]就行啦(如果l<=k<=r ),就要给f[1]++#include<cstdio>#include<cstring>#define maxl 100010#define mod

2015-07-01 21:45:40 598

原创 [欧拉函数]Bzoj2705 Longge的问题[SDOI2012]

因为要求gcd之和，所以枚举约数k然后找有多少个gcd是k的，因为gcd(i,n)=k,所gcd(i/k,n/k)=1,所以只要与n/k互质的都行，就用欧拉函数找出与n/k互质的个数，因为i是n的约数，所以n/i也是i的约数，就只要枚举到sqrt（n）就行啦#include<cstdio>#include<cmath>long long n,ans;int m;long long phi(lon

2015-07-01 21:43:38 535

原创 Bzoj1485 有趣的数列[HNOI2009]

实就是从左往右扫每个数，把放在奇数项看作入栈，偶数看作出栈，然后就是一个求卡特兰数，但有些直接枚举质因数会超时,筛素数的时候对于每一个数只要记录它的最小质因子，然而出来在把每个数分解掉，比如60，最小质因子是2，变成30，又是2，边15，是3，变5，然后就没了，原来那种质因数分解发会枚举出许多不能再除的数，所以这种快一些#include<cstdio>#include<cstring>#defi

2015-07-01 21:41:34 586

原创 [容斥原理]Bzoj1042 硬币购物[HAOI2008]

先用递推算出硬币无限的所有方案数，然后,因为一个硬币超限的最小个数是d[i]+1,能使硬币超限的最小钱数是c[i]*(d[i]+1),对于一个硬币超限的所有情况就是f[sum-c[i]*(d[i]+1)]，然而4个硬币超限的情况并是有重叠关系，这就用到容斥原理的关系，奇数个就-,偶数个就加#include<cstdio>#include<cstring>#define maxn 100010i

2015-07-01 21:40:31 567

原创 [费马小定理]Bzoj3260 跳

坑定先走长的再走短的，当要转弯时，此时的方案书为（m为长边）c(m,m)，而c(m+i,m)=c(m+i-1,m)/i*(m+i)，而a^(mod-1)与1（% mod）是同余的，/i就是乘i^(mod-2)然后再mod就没区别了，这样就能的到向上（右）走的每个数，把他们加起来就行了#include<cstdio>#include<cstring>#define mod 1000000007l

2015-07-01 21:38:16 758

原创【组合】Bzoj3142 数列[HNOI2013]

设相邻两项的差值为a[i]=A[i+1]-A[i]，因为第一天股价多少是不知道的，就=n-a[1]-a[2]..a[k-1]，而对于每一个a[i]又有1到m种选择，于是就是所有的序列加起来，把n拆出来 n*m^k-1 - 后面序列的所有数，后面是有m^（k-1）*（k-1）的数，但m个数均出现相等的次数，则每个数均出现m^(k-2)*(k-1),m个数的和公式得综合为n*m^(k-1)-m(m+1)

2015-07-01 21:36:29 1001

原创 [组合]Bzoj1008 越狱[HNOI2008]

所有情况就是m^n,不可能都情况就是对于2-n个人都不能与前一个人一样每个人有m-1种选择，于是不可能的总数就是（m-1）^(n-1)，相减就行啦#include<cstdio>#include<cstring>#define mod 100003long long n,m;long long ans;void prework(){ scanf("%lld %lld",&m,&n)

2015-07-01 21:33:44 540

原创【莫比乌斯反演】【bzoj2301】problem b

f(i)x∈[1,n],y∈[1,m],所有gcd=i的点对的数量,我们发现并不好求，于是再设F(i)x∈[1,n],y∈[1,m],所有满足i整除gcd(x,y)的点对的数量，F显然就是(n/i)(m /i),f根据莫比乌斯反演+gcd化为1简化一下，枚举倍数就变成了1到n/k啦，然后就用下标分块优化就可以求啦 #include<cstdio>#include<cstring>#define

2015-07-01 21:31:54 522

原创 APIO 2015

DAY0 晚上7点上的火车，火车上玩了几把dota，然后看了一下电影，软卧的床真心舒服啊，只是房间里没有插座，我只能强行到厕所旁边给笔记本充电，12点就睡了。DAY1 我一睁开眼，看见手表5点半，我就想这我手表什么时候坏的，昨晚看时间还是对的啊，我去年买的表怎么就坏了。。。。但天太亮了，我也睡不着了，就出门去隔壁看看，一出门就碰到GB，我问几点了，结果他说5点半- -，我擦，北方这天亮

2015-07-01 21:25:09 1269 2

原创【莫比乌斯反演】[BZOJ 2820 YY的GCD]

从黄学长那盗来的图可以线筛预处理F,或者暴力枚举质数（这个复杂度我不能确定）,按照素数粗略个数n/logn以及调和级数求和logn来看暴力的复杂度接近On 处理完F以后就是喜闻乐见的下底函数分块#include#include#define maxl 10000001int n,m;int no[maxl],p[maxl],mu[maxl];long long sum[m

2015-07-01 21:00:05 686 2

原创 POI代码等式

Poi1999基因片段问题描述有这样一段遗传基因K，它是由一系列的自然数组成：K=a1,a2,a3,a4……am。在该段基因中，连续的两个自然数被称做它的“特征”。例如对基因段：8, 5, 1, 4, 2, 3(5,1)就是它的“特征”之一，而(4,3)则不是。L教授正在研究这样一段长度未知的奇特基因段，并且已经成功地分析出了该基因段的大多数“特征

2015-07-01 20:47:38 838 1

binarycopycode