Claris' Blog

摘要： 询问[l,r]实际就是询问$\varphi(该区间内所有数的乘积)$$\varphi(n)=n\times...\times(1-\frac{1}{p_i})$由于素数只有60个，所以只需用线段树维护区间内所有数字的积，以及压位表示的出现过的素数即可$O(\log n)$每次操作#includety... 阅读全文

摘要： 左右与右左是两个独立的问题设f[i]表示i时刻左右车道减少一条的答案g[i]表示i时刻右左车道增加一条的答案ans=min(f[i]+g[i+r])计算f[i]：首先暴力计算出f[m+1]，同时记录下每个时刻刚开走的车的数量now[i]从m到1计算f[i]，如果该时刻开走的车数不足n1+1，则无影响... 阅读全文

摘要： 左上角和右下角不四连通等价于左下角和右上角八连通枚举正方形的左上角，先二分出最大的边长，使得里面不含障碍物然后再二分出最小的边长，使得两部分连通，用前缀和判断这题WA了好久…一直对拍都没问题…于是去看原题，发现有SPJ…然后改了个枚举姿势就过了…#include#define N 1510int n... 阅读全文

摘要： 对于树边显然只需要减少权值，对于非树边显然只需要增加权值设i不为树边，j为树边X[i]:i增加量X[j]:j减少量C[i]:修改1单位i的代价对于每条非树边i(u,v)，在树上u到v路径上的所有边j都需要满足$W_i+X_i\geq W_j-X_j$即$X_i+X_j\geq W_j-W_i$最后我... 阅读全文

摘要： 写个cmp然后sort就好了。cmp的话，需要快速知道两个串的lcp，于是倍增+Hash即可。#include#includetypedef long long ll;const int N=100010,M=18,P=1000000009;int n,i,j,k,f[N][M],a[N];char... 阅读全文

摘要： 第一问：打表可得规律：当且仅当x&(xtypedef long long ll;struct mat{ ll a[2][2]; mat(){a[0][0]=a[0][1]=a[1][0]=a[1][1]=0;} mat operator*(mat b){ mat c; for(in... 阅读全文

摘要： 有上下界的网络流T向S连容量为正无穷的边，将有源汇转化为无源汇每条边容量减去下界，设in[i]表示流入i的下界之和减去流出i的下界之和新建超级源汇SS,TT对于in[i]>0的点，SS向i连容量为in[i]的边对于in[i]0)$，则有解再求出以S,T为源汇的最大流即为答案#includeconst... 阅读全文

摘要： 将S1与S2用#号拼接在一起形成S串将S3与S串跑KMP求出S3在S串中每次出现的位置l[i]对于S串每个后缀i，求出f[i]表示该串不包含S3串的最长前缀然后求出S串的后缀数组先从小到大扫描后缀数组，同时维护一个tmp表示S2中的串与现在的串的最长公共前缀，且没有出现S3，如果碰到一个S1的后缀，... 阅读全文

摘要： 怎么3776又换题目了…换题目了…题目了…目了…了…SCC缩点后只有入度或者出度为0的点必须要放警察局假设一共有t-1个入度或者出度为0的SCCq[1]-q[t-1]表示这些SCC中点的个数q[t]表示剩下的点的个数f[i][j]表示前i个部分放了j个警察局的方案数，f[i][j]=sum(f[i-... 阅读全文

摘要： 从大到小枚举答案与T串的lcp，然后贪心#include#includechar s[5010],t[5010],ans[5010];int n,m,i,j,a[128],b[128],flag;int main(){ scanf("%s%s",s,t); n=std::strlen(s),m=... 阅读全文

摘要： Manacher求出所有极长回文子串后，得到一堆线段，转化成线段覆盖问题预处理出g[i]表示左端点不超过i的右端点的最大值贪心地线段覆盖即可时间复杂度$O(n)$#include#include#include#define N 100010using namespace std;char a[N]... 阅读全文

摘要： 从左边不能到达右边当且仅当存在一条与上下底边相连的分割线将它们分开设下底边为S，上底边为T，每个元件作为点，有公共部分的两个点互相连边最后拆点求最小割#include#include#includeusing namespace std;const int N=620,inf=~0U>>2;int ... 阅读全文

摘要： $C_{N+M}^N=\frac{(N+M)!}{N!M!}$考虑求出$ans\bmod 10^9$的值$10^9=2^9\times5^9$以$2^9$为例，先预处理出$1$..$2^9$中不是2的倍数的数的前缀积s[]，显然$n!\bmod 2^9$有着长度为$2^9$的循环节将答案表示成$a\... 阅读全文

摘要： 这题显然是一个最小斯坦纳树的模型，直接上模板即可，就是正六边形比较恶心。#include#include#include#include#define rep(i,n) for(int i=0;iP;int R,n,m,k=4,f[N][1Q;int read(){ char ch; while... 阅读全文

摘要： 将序列分成$\sqrt{n}$块，预处理出每两块之间的逆序对数，以及ap[i]表示前i块内数字出现次数的树状数组预处理：$O(n\sqrt{n}\log n)$修改时，ap[i]可以在$O(\sqrt{n}\log n)$复杂度内完成修改，然后考虑修改的位置对答案的贡献，可以发现相当于某一行、某一列... 阅读全文

摘要： 由于没有换根操作，所以直接用Splay维护DFS括号序列即可。对于查询1到x路径和，就是已入栈的减去已出栈的。子树加相当于区间加。换父亲相当于一个区间的移动。本题非常卡常数，一开始直接扒3153的代码狂T不止，改写Splay后也是狂T不止，最后把Splay从数组储存改成结构体储存才AC。#inclu... 阅读全文

摘要： 从低位到高位数位DP，f[i][j][k]表示已经填了后i位，转化的数字为j，后i位与x后i位的大小关系为k的方案数。#includeconst int N=202,B=7,P=(1P)c.x[i+1]++,c.x[i]&=P; return c; } Num operator-(Num ... 阅读全文

摘要： 设第$i$条直线的解析式为$a_ix+b_iy+c_i=0$，$d_i=a_i^2+b_i^2$，则要求一个点$(x,y)$，使得$\sum\left(\frac{\left|a_ix+b_iy+c_i\right|}{\sqrt{d_i}}\right)^2$最小。 假设$x$固... 阅读全文

摘要： 树上莫队，有几个坑：a可能等于b数据范围是100000#include#include#include#define N 100010#define K 17using namespace std;inline void read(int&a){char c;while(!(((c=getchar(... 阅读全文

摘要： l[i],r[i]表示站在i点往左往右走能看到的最高峰，用栈维护凸壳求出h[i]表示i点能看到的最高峰的高度a[i],b[i]表示i点往左往右走时反悔的点，即第一个h[j]>h[i]的j，用单调栈求出然后建树DFS一遍求出答案#include#define N 200010typedef long ... 阅读全文