2016年06月_cervoliu

原创【JZOJ 3466】选课 select

Description你真的认为选课是那么容易的事吗？HYSBZ的ZY同志告诉你，原来选课也会让人产生一种想要回到火星的感觉。假设你的一周有n天，那么ZY编写的选课系统就会给你n堂课。但是该系统不允许在星期i和星期i+1的时候选第i堂课，也不允许你在星期n和星期一的时候选第n堂课。然后连你自己也搞不清哪种选课方案合法，哪种选课不合法了。你只想知道，你到底有多少种合法的选课方案。对于第i组数据，n

2016-06-27 19:51:33 869

原创【JZOJ 3463】军训 training

DescriptionHYSBZ 开学了！今年HYSBZ 有n 个男生来上学，学号为1…n，每个学生都必须参加军训。在这种比较堕落的学校里，每个男生都会有Gi 个女朋友，而且每个人都会有一个欠扁值Hi。学校为了保证军训时教官不会因为学生们都是人生赢家或者是太欠扁而发生打架事故，所以要把学生们分班，并做出了如下要求： 1.分班必须按照学号顺序来，即不能在一个班上出现学号不连续的情况。 2.每个学生

2016-06-25 14:23:42 811

原创【JZOJ 3851】reward

DescriptionBsny最近公司运作不佳，本年度利润才m元，但员工的奖金还是要发的，公司有n个员工，怎么发奖金这个完全由老板Bsny自己决定。Bsny想要么把这m元全发了，激励一下员工，但具体怎么分配方案有很多。比如m=1, n=2, 那么可以员工1发1元，员工2发0元；也可以员工1发0元，员工2发1元，有两种方案。但其实，Bsny还是有点吝啬的，他想这m元不一定全部作为奖金，可以部分留给

2016-06-24 20:26:53 1870

原创【组合数取模】【SDOI2013】【BZOJ 3129】【JZOJ 3214】方程

Description给定方程 X1+X2+…+Xn=mX_1+X_2+…+X_n=m 我们对第1…n11\ldots n_1个变量进行一些限制： X1≤A1X_1≤A_1 X2≤A2X_2≤A_2 …… Xn1≤An1X_{n1}≤A_{n1} 我们对第n1+1,n1+2…n1+n2n_1+1,n_1+2\ldots n_1+n_2个变量进行一些限制： Xn1+1≥An1+1

2016-06-23 19:29:01 1197

原创拉格朗日插值法

前言在数值分析中，拉格朗日插值法是一种多项式插值方法。拉格朗日插值法下文，本蒟蒻可能有说的不妥甚至不对的地方，欢迎大神来打脸。定义多项式函数f(x)f(x)f(x)，则点值运算就是给定xxx，求y=f(x)y=f(x)y=f(x) 插值是点值的逆运算，就是给定nnn次多项式的n+1n+1n+1个点值表达(xi,yi)(xi,yi)(x_i,y_i)，要求出f(x)f(x)f...

2016-06-20 21:31:43 3392

原创【NOI 2010】【BZOJ 2005】【JZOJ 2225】能量采集

Description 对于100%的数据：1≤n,m≤100,0001 ≤ n, m ≤ 100,000Analysis首先，让我们想一想，线段上会碰到多少个植物其实是什么？假设坐标为(dn,dm),d=gcd(n,m),x=nd,y=md(dn,dm),d=gcd(n,m),x=\dfrac{n}{d},y=\dfrac{m}{d}。那么对于所有的xi,yi(xi≤n,yi≤m)x

2016-06-17 17:49:13 511

原创【JZOJ 4388】染色

Description N≤105N\leq 10^5Analysis暴力竟然跑得极快极快，理论O(n2)O(n^2)的比O(nlog2n)O(nlog_2n)的快得多得多。。。其实因为这里的NN只有10510^5，而我的暴力又是找过重心以后的，所以实际情况比O(n2)O(n^2)要小得多。 GDOI2015有一道题，也是N≤105N\leq 10^5，可是标程是O(n2)O(n^2)的，而

2016-06-16 21:25:33 889

原创【JZOJ 3430】 DY引擎

BOSS送给小唐一辆车。小唐开着这辆车从PKU出发去ZJU上课了。众所周知，天朝公路的收费站超多的。经过观察地图，小唐发现从PKU出发到ZJU的所有路径只会有N（2<=N<=300）个不同的中转点，其中有M（max(0, N-100) <=M<=N）个点是天朝的收费站。N个中转点标号为1…N，其中1代表PKU，N代表ZJU。中转点之间总共有E（E<=50,000）条双向边连接。每个点还有一个

2016-06-16 12:25:25 677

原创差分约束系统

差分约束系统的学习小记前言这是一篇学习小记。（如果你很赶时间直接跳过前言）由于本蒟蒻之前没有接触过这类题目，所以最近被一道题虐了一下，然后学习了。其实这个名字很耳熟，只是在我小的时候，听说了，却没有去学习。好吧，其实现在才这个东西是个非常基础的东西，但是它的思路极为巧妙。例题首先来看一道题。有一段长度为n的序列，点分为特殊点和一般点。有q个提示，提示形如(x,y,c)的三元组，表示序

2016-06-15 21:58:33 579

原创【BZOJ 2818】 gcd（附φ的线性筛法预处理）

Description给定整数N，求1<=x,y<=N且Gcd(x,y)为素数的数对(x,y)有多少对. 1<=N<=10^7Analysis设gcd(px,py)=p,pgcd(px,py)=p,p为素数。满足上式的x,yx,y必然互质。因此 Ans=∑1≤p≤n,p为素数(1+2∗∑x=2⌊np⌋φ(x))Ans=\sum_{1\leq p\leq n,p为素数}(1+2*\su

2016-06-15 12:32:26 849

原创【JZOJ 4489】【GDOI 2016 Day2】第一题 SigemaGO

Description Analysis这题在GDOI上想到了，可是只有20，so sad 对于那个插近道，只需要做limlim次最短路即可。把图分层，插近道相当于从上一层的图走到这一层的图。当然，这一层本身的最短路也要跑的。 SP(b)FA的时间复杂度O(lim∗n∗一个神奇的常数)O(lim*n*一个神奇的常数),dijkstra的是O(lim∗n∗log2n)O(lim*n*log_

2016-06-13 20:04:00 632

原创【JZOJ 4496】【GDSOI 2016】第一题互补约数（两种解法）

Description求∑i=1n∑d|igcd(d,id)\sum_{i=1}^n\sum_{d|i}gcd(d,\dfrac id) 对于100%的数据，n≤1011n\leq 10^{11}。Analysis方法一这是一个不需要莫比乌斯反演的方法。 ∑i=1n∑d|igcd(d,id)\sum_{i=1}^n\sum_{d|i}gcd(d,\dfrac id) =∑i∗j≤ngcd(i,

2016-06-13 19:33:05 891

原创【JZOJ 4488】【GDOI 2016 Day1】第四题疯狂动物城

Description给出一个N个节点的数，和M次操作。每次操作的类型如下： 1,x,y,z，将x到y的路径上的ai加上z 2,x,y，询问x到y的路径上，ai*(1+2+..+n-i)的和 3,x，将所有的a变更回第x次修改之后的状态。强制在线。 N,M<=10^5.Analysis裸题——码农一道！可持久化线段树+树链剖分兹磁树上的区间查询&修改。代码量飞飞飞飞起~

2016-06-12 20:58:39 671

原创【JZOJ 4486】【GDOI 2016 Day1】第二题最长公共子串

Description给你两个字符串S,T 给你N个区间，S串中这些区间内的字符可以随意改变在区间中的位置。求S,T可以达到的的最长公共子串长度。Analysis可以用双指针移动，我打的是dp。设f[i][j]f[i][j]表示S串的第i个区间与T串的第j个字符的LCS。预处理一堆东西，转移的话自行脑补。Code#include<cstdio>#include<cstring>

2016-06-12 20:50:24 800

原创 JZOJ 4485【GDOI 2016 Day1】第一题中学生数学题

Description给定商品成本 p0 和系数 k ，常数 n0 ，设价格为 p ，购买商品人数满足n1=⌊n0−k∗p⌋n1 =⌊ n0 -k*p⌋ 其中 p ， k 为实数， n0 为整数。题目分为两部分： 1）可以设置一个价格，求最大利润。 2）可以设置两个价格，p1，p2 p1 ， p2 ,有 n1=⌊n0−k∗p1⌋n1 = ⌊n0-k*p1 ⌋ 的人按 p1 购买，n2=

2016-06-12 20:42:31 669

原创【JZOJ 4503】异或树动态点分治

Description 这只熊孩子会搞TT次破坏,点权为EiEi。对于100%的数据2<=n,T<=30000,Ei<=163842<=n,T<=30000,Ei<=16384AnalysisEi<=16384=214Ei<=16384=2^{14}，这启示了我们什么？ xorxor怎么解决？可以拆位做，最后每个二进制位合并。这题如果不搞破坏，点分治可做。可是搞破坏，难道要打动态点分

2016-06-12 20:13:49 1150 2

原创【JZOJ 3432】服务器斜率优化DP常见问题&详细解答

Description我们需要将一个文件复制到n个服务器上，这些服务器的编号为S1, S2, …, Sn。首先，我们可以选择一些服务器，直接把文件复制到它们中；将文件复制到服务器Si上，需要花费ci > 0的置放费用。对于没有直接被复制文件的服务器Si来说，它依次向后检查Si+1, Si+2, …直到找到一台服务器Sj：Sj中的文件是通过直接复制得到的，于是Si从Sj处间接复制得到该文件，这种复

2016-06-11 19:22:22 829 1

Description某城市的街道呈网格状，左下角坐标为A(0, 0)，右上角坐标为B(n, m)，其中n >= m。现在从A(0, 0)点出发，只能沿着街道向正右方或者正上方行走，且不能经过图示中直线左上方的点，即任何途径的点(x, y)都要满足x >= y，请问在这些前提下，到达B(n, m)有多少种走法。 100%的数据中，1 <= m <= n <= 5 000Analysis因为我是蒟蒻

2016-06-11 19:18:04 882

原创 [JZOJ junior 2042] SuperPow

题目描述众所周知，a^b表示a的b次幂。例如：2^3=2*2*2=8。一天，某只肥皂很无聊，于是在纸上写了形如a^b的式子玩。FL见到了，过来一起玩。突然，FL脑洞一开：我给你普及一个符号“^^”，叫做超级幂。a^^b表示a^(a^(a^(a^a^(…))))，共b个a。例如2^^3=2^(2^2)=2^4=16。那么你知道a^^b的个位数是多少么？肥皂表示一脸懵逼，于是一旁看热闹的YDL出

2016-06-04 16:34:37 629

原创 [JZOJ 3427] 归途与征程

Description对于100%的测试点，1<=N<=100,1<=M<=100000。Analysis我们可以把A串视为一些串中间隔着一些星号。显然，A串头尾都没有星号是有星号的特殊情况。因为无星号可以跳过头尾的串变成头尾都是星号。下图，上者A，下者B。我们可以对于A串中的小串与B串做一次KMP，搞一个bz[i][j]bz[i][j]表示第i个小串是否能与B的第j个位置往后相应长

2016-06-01 12:41:58 630

LYD729