2015年07月_Lcomyn

原创【BZOJ 3679】数字之积

人生第一道数位DP，首先对于每位数的乘积，有一个很显然的转移方程 d[i][j]d[i][j]表示ii位数乘积为jj的方案数，则有 d[i][j]=∑1≤k≤9,k|jd[i−1][k/j]d[i][j]=\sum_{1 \le k \le 9,k|j} d[i-1][k/j] 然而我们发现j可能很大，但经过实验发现只有5000余个，于是我们可以吧第二维下标换成在数表中的排名，单个递推就可做了

2015-07-30 16:34:58 1841

原创【SDOI 2013】【BZOJ 3130】费用流

这题的第一问是一个裸最大流，不多说了，关键在于第二问。首先有一个结论，Bob一定把费用加在一条边上，于是我们就可以二分每条边的流量，验证是否能得到最大流即可code：#include#include#include#define inf 0x7fffffffusing namespace std;int n,m,p,e,ei;int point[101],next[2001];

2015-07-29 16:55:24 904

原创【BZOJ 1148】【CTSC 2008】挂缀【BZOJ 1 148】【CTSC 2008】挂缀

【BZOJ 1 148】【CTSC 2008】挂缀这题显然是个贪心，然而我们应该如何贪才能得到最优解= =。。。。假设我们按重量升序贪心，那我们可以得到反例（假设在挂缀底部）：设有i,j,j且W_j那么当C_i时，应该选i而不是j假设我们按拉力升序贪心，依旧可以得到反例（假设在挂缀顶部，S为链以下重量和）：设有i,j,j且C_j那么当C_i时，应该先选i当第二高的，而非链顶

2015-07-28 19:57:52 1052

原创【BZOJ 1146】【CTSC 2008】网络管理network

一句话题意，树链上带修改区间第k大感觉可以dfs+主席树O(nlog2n)O(n\log^2n)过掉，但我不会写= = 于是写的线段树套平衡树+链剖+二分（修改O(nlog3n)O(n\log^3n),查询O(nlog4n)O(n\log^4n)慢了好多啊QAQ）这里简单介绍一下区间第K大做法，对于每个线段树所”管辖“的范围，建一棵相应范围内的平衡树（我用的Treap）；修改时，修改

2015-07-27 17:04:32 1002

原创后缀数组学习笔记

后缀数组是一个处理字符串的有力工具，基本用途有模式串匹配和后缀最长公共前缀；例题【JSOI 2007】【BZOJ 1031】字符加密ciper 后缀数组的躶体，现将字符串复制成环，然后输出sa[i]<lensa[i] \lt len的s[sa[i]−1]s[sa[i]-1]即可，code：#include<iostream>#include<cstdio>#include<cstrin

2015-07-25 16:53:21 778

原创【NOI 2010】能量采集&&超级钢琴

其实这两个题一点关系都没有，同一天做的，写在一起= = 能量采集 ans=∑x=1n∑y=1m2∗gcd(x,y)+1ans=\sum_{x=1}^{n}\sum_{y=1}^{m} 2*gcd(x,y)+1 ∑x=1n∑y=1m[d|gcd(x,y)]=⌊nd⌋∗⌊md⌋\sum_{x=1}^{n}\sum_{y=1}^{m}[d|gcd(x,y)]=\lfloor \frac {n}{d}

2015-07-11 17:03:36 979

原创【BZOJ 2818】Gcd

这题一开始我竟然想要用与能量采集差不多的思路去做= =（no zuo no die,why you try?）有个显然的转化 ∑nx=1∑ny=1[gcd(x,y)==P]=∑⌊n/P⌋x=1∑⌊n/P⌋y=1[gcd(x,y)==1]=2∗∑⌊n/P⌋i=2ϕ(i)+1\sum_{x=1}^{n}\sum_{y=1}^{n} [gcd(x,y)==P]= \sum_{x=1}^{\lfl

2015-07-11 16:48:57 838

原创【SHOI 2008】【BZOJ 1018】堵塞的交通traffic

发现很多有关区间联通问题的题线段树都可以做，不论是1列还是2列。这题的易错点在于两个端点x,yx,y之间的可行路径不一定在区间[x,y][x,y]之间，有可能是从其他区间绕过来的，想明白这一点，区间的更新又与1列的情况基本类似，这题就很好做了。 code：#include<iostream>#include<cstdio>#include<cstring>#include<algorithm

2015-07-11 16:38:11 1665

原创【JSOI 2008】【BZOJ 1014】火星人prefix

这题其实很奇怪，一眼看出fhq treap+字符串hash。。。结果在BZOJ被卡常了，结果学长告诉我了一些优化，有些自己还不会写，只写了inline和吧unsigned long long改为unsigned int。。。结果时间少了一半，谁能告诉我为什么。。。#include<iostream>#include<cstdio>#include<cstring>#include<cs

2015-07-04 21:56:28 935

原创【BZOJ 4128】 matrix

题目其实很简单，BSGS可以一眼看出来，但这就需要矩阵求逆，然而我并不会。。。于是发现了一种BSGS的非求逆方法，借此介绍一下。 aix−y=b(modP)a^{ix-y}=b (mod P) 推出 aix=b∗ay(modP)a^{ix}=b*a^y(mod P) 因此我们只需要对每个b∗ayb*a^yhash一下就好了，然后穷举ii,验证是否存在即可。code:#include<cst

2015-07-01 20:45:30 1136

Lcomyn的专栏