20210813 a,b,c
考场
稍微想了想发现 T1 是 sb 题，枚举矩形的三个边界，右边界双指针扫就行了，T2 八成 DP，T3 感觉非常不可做
T1 犹豫了一下要不要算补集，感觉直接写也不难，就打消了这个念头
T2 只会 \(O(nma\log a)\) 的 DP，有 65pts，如果卡常+数据水也许能再混 15pts。感觉剩下的分要用数学，T3 部分分太香就没多想
T3 手玩样例发现不考虑重边的话就是一棵树，那么可以i拎出来唯一的路径 DP，\(O(qn)\) 有 47pts，放弃梦想了。
大概 7.50 开始码 T1 ，8.05 过了大样例，拍上后就没管
T2 又想了一会怎么优化，然后写了暴力 DP。发现 sub4 要跑个 8s 左右，试图卡常，然而并没有什么用，拍上看 T3 了
写完暴力还有 1.5h+，感觉发现树的部分是《[SDOI2011]染色》的弱化版，推广到仙人掌上，每条边有两种选择，带 \(2^2\) 常数分类讨论就行了，还有 1h+ 开始码。
码到剩 40min 时发现分类讨论的细节特别多，保树的分吧。
然后发现边权转点权后倍增的区间非常阴间，加上时间不多导致脑子一片混乱，最后 10min 还没过样例，弃了
res
rk5 90+65+46
T1 在 \(l=0\) 时可能一段 \(s_i\) 都是 \(0\)，右边界可能比左边界还小。。。
T3 数组开小 RE 了 1pt
rk1 高俊垚,张泽阳 100+100+21
rk33 刘荣信 30+0+47
总结
这场的部分分拿得比较全，主要就是 T2 低太多的（前几都 A 了），事实上最近几场打的都不太好，不仅稳定不在前 3，还翻车考过 rk10，rk1 更是好久没拿了。趁机反思一下。
自己感觉最大的问题就是想不出来题，除了 sb 题和没有思维含量的 DS 外几乎没有在考场上 A 题，DS 有时还会调不出来导致整场爆炸。
从上次颓废后状态有所松懈，还是要绷紧啊。DP 有几天没写了，最近的早晨也比较荒废，弄完数论开组合数学吧。
最近高强度考试没时间刷题，那就更要珍惜考场上思考和考完改题的时间，不仅想清楚每个题的做法，更要清楚的是怎么想到这么做，这个题用到的思考方式/性质/结论能不能迁移到其他题上，tricks 的博客不要流于形式。
昨天和今天还暴露出一个问题，就是总在比赛快结束时 rush 一个题，随着时间推移得分不断退化，最终爆 \(0\)。除了考试快结束时心态不稳外，码力不行也是很重要的一个原因（从今天下午调 T3 也能看出来）。
前几天 hkh 写猪国杀，我认为写大模拟没有意义，考场上更重要的还是在较短时间内写/调出中等长度的代码，不过或许可以找几道毒瘤 DS 做做，最近没有时间就等回了 sdfz，先记在这里。
努力，奋斗！——《喜剧之王》
sol
叒没时间写了
T1
const int N = 35, M = 5e4+5;
int n,m,mn,mx;
char s[N][M];
int cnt[M],sum[M];
LL ans;
signed main() {
// freopen("a.in","r",stdin);
// freopen("a.out","w",stdout);
read(n,m);
For(i,1,n) scanf("%s",s[i]+1);
read(mn,mx);
For(i,1,n) {
mem(cnt,0,m);
For(j,i,n) {
For(k,1,m) cnt[k] += s[j][k]-'0', sum[k] = sum[k-1]+cnt[k];
int l = 1, r = 1;
For(k,1,m) {
ckmax(l,k), ckmax(r,k);
while( l <= m && sum[l]-sum[k-1] < mn ) ++l;
while( r < m && sum[r+1]-sum[k-1] <= mx ) ++r;
if( l > r ) break;
ans += r-l+1;
}
}
}
write(ans);
return iocl();
}
T2
容斥
const int N = 25, M = 1e5+5, mod = 1e9+7;
int n,m,a[N][M];
int ans,cnt[N][M],f[M];
void ckadd(int &x,int y) { x += y; if( x>=mod)x-=mod; else if(x<0)x+=mod;}
void ckmul(int &x,int y) { x = (LL)x * y %mod; }
signed main() {
read(n,m);
For(i,1,n) For(j,1,m) read(a[i][j]), ++cnt[i][a[i][j]];
For(i,1,n) For(j,1,1e5)
for(int k = j+j; k <= 1e5; k += j) cnt[i][j] += cnt[i][k];
rFor(i,1e5,1) {
f[i] = 1;
For(j,1,n) ckmul(f[i],cnt[j][i]+1);
--f[i];
for(int j = i+i; j <= 1e5; j += i) ckadd(f[i],-f[j]);
ckadd(ans,(LL)f[i]*i%mod);
}
write(ans);
return iocl();
}
T3
点分治
typedef vector<int> VI;
const int N = 1e5+5;
int n,m,qn,mm=1;
vector<PII> to[N];
PII q[N];
int ans[N];
namespace divi {
int rt,siz[N],mx[N],fa[N],dep[N],top[N],f[N][3][3];
bool vis[N];
VI q[N],col[N];
vector<pair<int,VI>> e[N];
void findrt(int u,int fa) {
siz[u] = 1, mx[u] = 0;
for(auto i : e[u]) {
int v = i.fi; if( v == fa || vis[v] ) continue;
findrt(v,u);
siz[u] += siz[v], ckmax(mx[u],siz[v]);
}
ckmax(mx[u],siz[0]-siz[u]);
if( mx[u] < mx[rt] ) rt = u;
}
void build(int u,int fa) {
vis[u] = 1, divi::fa[u] = fa, dep[u] = dep[fa]+1;
for(auto i : e[u]) {
int v = i.fi; if( vis[v] ) continue;
siz[ rt=0 ] = siz[v], findrt(v,u), build(rt,u);
}
}
int lca(int u,int v) {
if( dep[u] < dep[v] ) swap(u,v);
while( dep[u] > dep[v] ) u = fa[u];
while( u != v ) u = fa[u], v = fa[v];
return u;
}
void dfs(int u,int fa,int top) {
divi::top[u] = top;
for(auto i : e[u]) {
int v = i.fi; if( vis[v] || v == fa ) continue;
col[v] = i.se;
For(j,0,col[top].size()-1) For(k,0,col[v].size()-1) {
f[v][j][k] = -1e9;
For(l,0,col[u].size()-1)
ckmax(f[v][j][k],f[u][j][l]+(col[v][k]!=col[u][l]));
}
dfs(v,u,top);
}
}
void work(int u) {
vis[u] = 1;
for(auto i : e[u]) {
int v = i.fi; if( vis[v] ) continue;
col[v] = i.se;
For(j,0,col[v].size()-1) For(k,0,col[v].size()-1)
f[v][j][k] = j==k ? 1 : -1e9;
dfs(v,u,v);
}
for(int i : q[u]) {
int x = ::q[i].fi, y = ::q[i].se; if( x == y ) continue;
if( y == u ) swap(x,y);
if( x == u )
For(j,0,col[top[y]].size()-1) For(k,0,col[y].size()-1)
ckmax(ans[i],f[y][j][k]);
else
For(j,0,col[top[x]].size()-1) For(k,0,col[x].size()-1)
For(l,0,col[top[y]].size()-1) For(r,0,col[y].size()-1)
ckmax(ans[i],
f[x][j][k]+f[y][l][r]-(col[top[x]][j]==col[top[y]][l]));
}
for(auto i : e[u]) {
int v = i.fi; if( vis[v] ) continue;
siz[ rt=0 ] = siz[v], findrt(v,u), work(rt);
}
}
void main() {
siz[ rt=0 ] = mx[0] = n, findrt(1,0), build(rt,0);
For(i,1,qn) q[lca(::q[i].fi,::q[i].se)].pb(i);
mem(vis,0,n);
siz[ rt=0 ] = n, findrt(1,0), work(rt);
}
}
void dfs(int u,int fa) {
sort(to[u].begin(),to[u].end()),
to[u].erase(unique(to[u].begin(),to[u].end()),to[u].end());
for(int i = 0, j,v; j = i, i < to[u].size(); i = j+1)
if( (v=to[u][i].fi) != fa ) {
while( j+1 < to[u].size() && to[u][j+1].fi == v ) ++j;
divi::e[u].pb(MP(v,VI())), divi::e[v].pb(MP(u,VI()));
for(int k = i; k <= j && k < i+3; ++k)
divi::e[u].back().se.pb(to[u][k].se),
divi::e[v].back().se.pb(to[u][k].se);
dfs(v,u);
}
}
signed main() {
read(n,m);
For(i,1,m) {
int x,y,z; read(x,y,z);
to[x].pb(MP(y,z)), to[y].pb(MP(x,z));
}
read(qn);
For(i,1,qn) read(q[i].fi,q[i].se);
dfs(1,0);
divi::main();
For(i,1,qn) write(ans[i]);
return iocl();
}
