THUPC2019/CTS2019/APIO2019/PKUSC2019游记

最新推荐文章于 2019-07-20 00:49:00 发布

diaopang1934

最新推荐文章于 2019-07-20 00:49:00 发布

阅读量327

点赞数

文章标签：数据结构与算法

原文链接：http://www.cnblogs.com/suika/p/10959533.html

版权

THUPC2019/CTS2019/APIO2019/PKUSC2019游记

5.10

中铺，火车好晃啊

5.11

打了THUPC2019的练习赛，华容道好评(四个小兵，~~杠鸭!~~)

5.12

打了THUPC2019

开场爆了好几发D，然后不想写，丢给ZCY，cerr居然wa了好多发

然后ZCY看出B是sb题，我们俩不想写，丢给ZCY

然后gmy和我分别写了M和J两个签到题，三人开始分别开C,I,L

ZCY丢给我一个\(O(m)\)的式子让我优化，我推了推之后发现k才20

你容斥是\(O(k^2)\)的又不是\(O(2^k)\)的你优化这个干嘛

然后我写容斥部分，他写组合数部分，交上去，WA ？？？

woc这个模数太小了可能没逆元啊，然后改改改改改，交上去，过了

此时gmy的I已经会了，但是我们俩又没人想写，依然丢给ZCY

写完一直WA，交暴力也WA

我们都在喷卡精的时候交了一发puts("nan");

然后过了，良心SPJ，最后rk36滚粗了

然后是讲评环节，广告俩小时，讲题？？？不会的还是不会

高调三人南，点了zcy的大三元，体验被飞的快感

晚上到宾馆才发现I题读入搞错了，我是傻逼

回去之后改了L和K。

5.13

在食堂里看到了rvalue,刚想打招呼就不见了QAQ,还好在考场旁又遇到了

然后打了CTS2019Day 1

玩题答3h+俩暴力

出来一看分数还行

orz \(\color{red}{\mathrm{Itst}}\) 切t2

5.14

打了CTS2019Day 2

写了30+10+0

一看成绩10+10+0

傻逼了，打铁滚粗了，凸包都能写错。

原因大概是我把左下角那个基准点也丢进去排序了吧

5.15

这都能领奖，还是个ag，不亏

吃顿神仙火锅平均一人90rmb，吓尿了

5.16

lzh来了，气氛突然变得？了起来

然后我被拉入lol坑

5.17

上午讲厉害东西完全听不懂，直接导致我下午去都没去。

5.18

打了APIO2019，自闭了

开场写完仨暴力，开始搞t1(device)，想了想好像挺傻逼的。

就是拆环做区间覆盖就行了。

我写了5k特判，然后因为有四个模数(A,B,B+1,A/gcd(A,B+1))，调不出来了。

最后20+43+60滚粗了。

出来一问大家都203，txc聚聚更是提前2小时ak离场

那我垫底稳了啊

5.19

上午又讲厉害东西，不过感觉比前天那个好玩？

晚上颁奖，我居然还有块Cu，不亏不亏

这下noi赛事的三个牌牌都齐了，就剩noi2019去打铁了（雾）

5.20

忘记干啥了

5.21 ~ 5.24

复习

5.25

全世界只有我在北京还没去tho，甚至连live party都去不上.jpg

下午打了PKUSC2019Day1

看t1，好像是个逆序对傻逼题

看t2，好像是个厉害题

看t3，好像还是个傻逼题

先写t3，对每个点搞出来61个数然后判加和是否能组成就行了吧

搞了个bitset优化背包，拿了60分感觉再卡卡常就能过了(flg)

然后写t1，写个主席树上去就过了，此时不过1h，感觉很稳(flg)

开始刚t3，写了这么个乱搞：

当可用数字个数大于22时直接视为可行 (22是二分出来的)

跑得飞快啊，直接2.几s

然后又yy出一个能使常数减半的做法，不过需要手写bitset控制位数= =

于是我就打开自带的mingw开始看里面的bitset

看了一个小时左右，突然意识到自己在考试，好像浪费一个小时不是很值(很傻逼)

看看t2吧，想了好长时间才发现可以枚举一个点

枚举一个点，把小于的视为0，大于的视为2，0022之间不影响。

然后状态数就非常少了，写了一发卡了卡常就过了。

回去卡t3，加了好多玄学优化(到最后也没手写bitset)

突然发现最后一个点艹过去了(1992ms)，然后第一个点wa了？？？不过还好取最高还有93分..

把第一个点特判掉之后就过不去最后一个点了(20xxms)，感觉评测机上下浮动很正常，那我多交几发就过了(flg)。

然后32发交到最后，流局了

提前出考场假装ak还行

问群里大佬，好多都提前2h ak啊

这个t3只要在trie树上做就好了啊，有用的点一共n个，按深度存log个bitset就行了，wtcl这都想不到。

明天继续加油QAQ

5.26

看题，woc三个树题，要爆零了

看t1，写完n^2的dp之后看了看链的做法，然后就会了正解

一开始写了个线段树合并然后发现并不能标记下传时新建节点，此时1h过去了QAQ

想着换成线段树启发式合并或者splay启发式合并能好写一点，先去看t2t3

看t2，只会8分

看t3，看了好久才看明白弦图是啥玩意，暴力不会打，告辞

只能去写t1啦，这个t1好像可以每次只维护子树内有限制的颜色集合，这样标记下传时不用新建节点比较妙，写一半发现我的做法在求点积的时候需要先乘再除，然后乘0就GG了，判掉乘0之后已经又过了2h。

交一发，AcAcWaTle ，这啥玩意啊

发现m的范围好像不是2e5是1e9啊，那我空间还得再开大一点..诶怎么就mle了。

在想要不要换成splay启发式合并来优化空间和时间(此时已经把一开始的线段树合并扔掉了QAQ)

想了想还是算了，打暴力吧，这题过的人应该不多(flg)

然后玩t2二叉树的分..，写了个暴力发现这个\(m=\lceil n/4\rceil\)，求出\(m\)之后直接搜就搜出来答案了。

那后面的点应该也是这样做吧，就开始猜这个\(m\)是什么，然后没猜出来凉了

最后0.5h刚t3，感觉应该能找找规律，然后发现线弦图一定不存在长度大于等于4的环，写写写暴力，没调出来，滚了。

出来问成绩，好多大佬说今天比昨天简单？？？

47+49+0=96，被高一的ysy和gmy吊打50分，怎么人均切t1t2啊QAQ

然后又听说有好多人600... 我今天成功被区分了

听大佬讲题，t1做法都一样只是我调不出来别人能调出来，t2那个m其实就是链覆盖(叶子个数/2), t3注意到结论之后变成树上dp，每次乘一个两项的多项式，用分治fft即可。

5.27

颁奖颁奖

拿了个二等跑路

给高二的奖好少啊

初三和高二的一等奖线大概差了200分

在火车上写d2t1的线段树合并，甚至比启发式合并好写？？完全不用处理乘0的情况

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cstdlib>
#include <cmath>
#include <vector>
#include <iostream>
using namespace std;
typedef long long ll;
#define mod 998244353
#define N 400050
#define M 20000050
#define db(x) cerr<<#x<<" = "<<x<<endl
#define fov(i,x) for(i=0;i<(int)x.size();i++)
ll qp(ll x,ll y=mod-2) {ll re=1;for(;y;y>>=1,x=x*x%mod)if(y&1)re=re*x%mod; return re;}
int head[N],to[N],nxt[N],cnt,n,m,K;
int ls[M],rs[M],tot,siz[M];
int sum[M],tag1[M],tag2[M];
vector<int>V[N];
struct node {
    int rt;
    ll O,s;
    void gs() {
        s=(sum[rt]+(m-siz[rt])*O)%mod;
    }
}tr[N];
node X,Y;
inline void add(int u,int v) {to[++cnt]=v; nxt[cnt]=head[u]; head[u]=cnt;}
inline void giv1(ll v,int p) {
    if(!p)return ;
    sum[p]=sum[p]*v%mod; tag1[p]=tag1[p]*v%mod; tag2[p]=tag2[p]*v%mod;
}
inline void giv2(ll v,int p) {
    if(!p)return ;
    sum[p]=(sum[p]+siz[p]*v)%mod; tag2[p]=(tag2[p]+v)%mod;
}
inline void pushdown(int p) {
    if(tag1[p]!=1||tag2[p]) {
        giv1(tag1[p],ls[p]);
        giv2(tag2[p],ls[p]);
        giv1(tag1[p],rs[p]);
        giv2(tag2[p],rs[p]);
        tag1[p]=1; tag2[p]=0;
    }
}
inline void pushup(int p) {
    sum[p]=(sum[ls[p]]+sum[rs[p]])%mod;
    siz[p]=siz[ls[p]]+siz[rs[p]];
}
void update(int l,int r,int x,int &p) {
    if(!p) { p=++tot; tag1[p]=1; siz[p]=1; }
    if(l==r) { sum[p]=0; return ;}
    int mid=(l+r)>>1; pushdown(p);
    if(x<=mid) update(l,mid,x,ls[p]);
    else update(mid+1,r,x,rs[p]);
    pushup(p);
}
int merge(int l,int r,int x,int y) {
    if(!x&&!y) return 0;
    if(!y) {
        giv1(Y.O,x);
        return x;
    }
    if(!x) {
        giv1(X.O,y);
        return y;
    }
    if(l==r) {
        sum[x]=(sum[x]*sum[y])%mod;
        return x;
    }
    int mid=(l+r)>>1; 
    pushdown(x); pushdown(y);
    ls[x]=merge(l,mid,ls[x],ls[y]);
    rs[x]=merge(mid+1,r,rs[x],rs[y]);
    pushup(x); 
    return x;
}
void dfs(int x,int y) {
    int i;
    tr[x].O=1;
    fov(i,V[x]) update(1,m,V[x][i],tr[x].rt);
    for(i=head[x];i;i=nxt[i]) if(to[i]!=y) {
        dfs(to[i],x);
        X=tr[x], Y=tr[to[i]];
        tr[x].rt=merge(1,m,tr[x].rt,tr[to[i]].rt);
        tr[x].O=tr[x].O*tr[to[i]].O%mod;
    }
    tr[x].gs();
    ll s=tr[x].s;
    if(x!=1) {
        tr[x].O=(s-tr[x].O)%mod;
        giv1(mod-1,tr[x].rt);
        giv2(s,tr[x].rt);
    }
}
int main() {
    scanf("%d%d%d",&n,&m,&K);
    int i,x,y;
    for(i=1;i<n;i++) scanf("%d%d",&x,&y),add(x,y),add(y,x);
    for(i=1;i<=K;i++) scanf("%d%d",&x,&y),V[x].push_back(y);
    dfs(1,0);
    printf("%lld\n",(tr[1].s+mod)%mod);
}