摘要：
今年的重庆省选？ 具体就是，对于每次修改，A[p,q]这个位置， 设d=gcd(p,q) ,则 gcd为d的每一个格子都会被修改，且他们之间有个不变的联系 A[p,q]/p/q==A[k,t]/k/t 所以只要记录对于gcd为d的所有格子，只要保存A[d][d]的值就可以了。 那么求前k行k列的值a 阅读全文
摘要：
省选前。发现某题czt的代码又一次奇短无比。。 然后我就学了一个新的stl库。 【根本原因是 我正处于并将长期处于C++萌新状态。。来熟练打一句话证明： for i:=1 to n do if i mod 2=1 then blabla....; priority_queue <int> s[100 阅读全文
摘要：
明明用了那么久了、 然而突然才发现tarjan真是个好东西。。 回顾了一下 存在博客里以后用吧。。 首先tarjan 主要是三个用途： 【有向图 强连通分量】 一个强连通分量G中的点可以互相到达，而一个G中的点可以到达的G外的点不能到达G中任意一个点。 {dfn[x]==low[x]} 1 void 阅读全文
摘要：
【逆元】： 若 a*b=1(%p) 则a是b在%p意义下的逆元。 则在%p意义下 一个数 除以a就等价于乘b 两种求逆元方法：1、若p为质数时，则有性质则a^(p-1)=1(%p) 。 即a*a^(p-2)=1(%p)。 所以a的逆元就是a^(p-2) 2、exgcd(a,b,x,y) 可求 ax+ 阅读全文
摘要：
难得学习一下C++的库。。【至今连map,vector都没用过的我。。 首先#include<bitset>或<bits/stdc++.h> 定义函数： { bitset <15> b; 15是bitset长度 之后我就写成n了，b是变量名 或： string s='01001'; bitset < 阅读全文
摘要：
【啊 首先 这是道权限题，然后本人显然是没有权限的 23咳3】 最近数据结构做的越来越少。。然后 就跟上次一样 ，一做就是三四种不同写法。 等价的题面： 最近GY大神在sc2的天梯中被神族虐得很惨，表示很不爽。ryz决定帮助GY大神练习散枪兵技术。GY生产了n*m个枪兵，并站成了一个大小为n*m的方 阅读全文
摘要：
题面自己上网查。 学了一下单纯形。当然 证明什么的 显然是没去学。不然估计就要残废了 上学期已经了解了 什么叫标准型。 听起来高大上 其实没什么 就是加入好多松弛变量+各种*(-1)，使得最后成为一般形式： 给定A[][]，求满足A[i][j]*Xj<=A[i][0];（0<i<=n,0<j<=m） 阅读全文
摘要：
n种珠宝。每种各1个。有价格ci元,美度vi。 要求分别输出1元到m元 可买的最大优美度。 整数 ：0<n<=10000000, 0<ci<=300,0<=vi<=10^9, 0<m<=50000; 之前 系统的看过有关背包的题目。。然而这个做法还没见过。 首先复杂度是 300*m*log(m)+n 阅读全文
摘要：
4589: Hard Nim Description Claris和NanoApe在玩石子游戏，他们有n堆石子，规则如下： 1. Claris和NanoApe两个人轮流拿石子，Claris先拿。 2. 每次只能从一堆中取若干个，可将一堆全取走，但不可不取，拿到最后1颗石子的人获胜。 不同的初始局面， 阅读全文