题解:

https://files.cnblogs.com/files/clrs97/old-solution.pdf

Code:

#include<cstdio>
#include<cstdlib>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=100010,M=1010,E=5500000,BUF=10000000,OUT=10000000;
unsigned int SA,SB,SC;
int Case,n,p,m,lim,X,i,j,k,op,x,y,a[N],g[N],nxt[N],st[N],en[N],dfn,s[N<<1],q[N];
int b[M],cb,f[M][M],ga[N<<1],gq[N<<1],vl[E],vr[E],w[E],NXT[E],ED;
ll ans[N<<1],bit[N<<1],bitx[N<<1],bity[N<<1],bitxy[N<<1];
char Buf[BUF],*buf=Buf,Out[OUT],*ou=Out;int Outn[30],Outcnt;
inline void read(int&a){for(a=0;*buf<48;buf++);while(*buf>47)a=a*10+*buf++-48;}
inline void read(unsigned int&a){for(a=0;*buf<48;buf++);while(*buf>47)a=a*10+*buf++-48;}
inline void write(ll x){
if(!x)*ou++=48;
else{
for(Outcnt=0;x;x/=10)Outn[++Outcnt]=x%10+48;
while(Outcnt)*ou++=Outn[Outcnt--];
}
}
inline void writeln(ll x){write(x);*ou++='\n';}
inline unsigned int rng61(){
SA^=SA<<16;
SA^=SA>>5;
SA^=SA<<1;
unsigned int t=SA;
SA=SB;
SB=SC;
SC^=t^SA;
return SC;
}
inline void addedge(int x,int y){nxt[y]=g[x];g[x]=y;}
inline void add(int&x,int l,int r,int z){vl[++ED]=l;vr[ED]=r;w[ED]=z;NXT[ED]=x;x=ED;}
void dfs(int x){
lim=st[x]=++dfn;
s[dfn]=1;
for(int i=g[x];i;i=nxt[i])dfs(i);
en[x]=++dfn;
s[dfn]=-1;
}
inline bool cmp(int x,int y){return a[x]<a[y];}
inline void initrect(){
int _=0,i,j;
sort(b+1,b+cb+1);
for(i=1;i<=cb;i++)if(b[i]!=b[i-1])b[++_]=b[i];
cb=_;
for(i=1;i<=cb;i++)for(j=1;j<=cb;j++)f[i][j]=0;
}
inline void maketag(int xl,int xr,int yl,int yr,int p){
xl=lower_bound(b,b+cb+1,xl-1)-b+1;
xr=lower_bound(b,b+cb+1,xr)-b+1;
yl=lower_bound(b,b+cb+1,yl-1)-b+1;
yr=lower_bound(b,b+cb+1,yr)-b+1;
f[xl][yl]+=p;
f[xl][yr]-=p;
f[xr][yl]-=p;
f[xr][yr]+=p;
}
inline void getrect(){
int i,j;
for(i=1;i<=cb;i++)for(j=1;j<=cb;j++){
f[i][j]+=f[i-1][j]+f[i][j-1]-f[i-1][j-1];
if(f[i][j]){
add(ga[b[i-1]+1],b[j-1]+1,b[j],1);
add(ga[b[i]+1],b[j-1]+1,b[j],-1);
}
}
}
inline void solve(int x,int mode){
if(mode==0){
b[++cb]=st[x]-1;
b[++cb]=en[x];
for(int i=g[x];i;i=nxt[i]){
b[++cb]=st[i]-1;
b[++cb]=en[i];
}
}else{
maketag(st[x],en[x],1,lim,1);
maketag(1,lim,st[x],en[x],1);
for(int i=g[x];i;i=nxt[i])maketag(st[i],en[i],st[i],en[i],-2);
}
}
inline void ins(int y,int p){
ll B=1LL*s[X-1]*p,C=1LL*s[y-1]*p,D=1LL*s[X-1]*s[y-1]*p;
for(;y<=lim;y+=y&-y){
bit[y]+=p;
bitx[y]+=B;
bity[y]+=C;
bitxy[y]+=D;
}
}
inline ll ask(int y){
ll A=0,B=0,C=0,D=0;int o=y;
for(;y;y-=y&-y){
A+=bit[y];
B+=bitx[y];
C+=bity[y];
D+=bitxy[y];
}
return A*s[X]*s[o]-B*s[o]-C*s[X]+D;
}
void work(){
read(n),read(p),read(SA),read(SB),read(SC);
for(i=1;i<=n;i++)g[i]=0;
for(i=1;i<=n+n;i++)ga[i]=gq[i]=bit[i]=bitx[i]=bity[i]=bitxy[i]=0;
dfn=ED=0;
for(i=2;i<=p;i++)addedge(i-1,i);
for(i=p+1;i<=n;i++)addedge(rng61()%(i-1)+1,i);
for(i=1;i<=n;i++)a[i]=rng61()%n+1;
dfs(1);
for(i=1;i<=lim;i++)s[i]+=s[i-1];
for(i=1;i<=n;i++){
q[i]=i;
if(en[i]>lim)en[i]=lim;
}
sort(q+1,q+n+1,cmp);
for(i=1;i<=n;i=j){
b[0]=0;
b[cb=1]=lim;
for(j=i;j<=n&&a[q[i]]==a[q[j]];j++)solve(q[j],0);
initrect();
for(k=i;k<j;k++)solve(q[k],1);
getrect();
}
read(m);
for(i=1;i<=m;i++){
read(op),read(x),read(y);
ans[i]=0;
int xl,xr,yl,yr;
if(op==1){
xl=xr=st[x];
yl=yr=st[y];
}else{
xl=1,xr=st[x];
yl=1,yr=st[y];
}
add(gq[xl-1],yl,yr,-i);
add(gq[xr],yl,yr,i);
}
for(X=1;X<=lim;X++){
for(j=ga[X];j;j=NXT[j]){
ins(vl[j],w[j]);
ins(vr[j]+1,-w[j]);
}
for(j=gq[X];j;j=NXT[j]){
k=w[j];
if(k>0)ans[k]+=ask(vr[j])-ask(vl[j]-1);
else ans[-k]-=ask(vr[j])-ask(vl[j]-1);
}
}
for(i=1;i<=m;i++)writeln(ans[i]);
}
int main(){
fread(Buf,1,BUF,stdin);
read(Case);
while(Case--)work();
fwrite(Out,1,ou-Out,stdout);
return 0;
}

  

SDOI2018:原题识别的更多相关文章

  1. #LOJ2564 SDOI2018 原题识别 主席树

    转载请注明原文地址:http://www.cnblogs.com/LadyLex/p/9057297.html 原题链接: 今天考试考了前天的SDOI考题 天啊我菜爆,只有T2拿了30分 然后考试后半 ...

  2. [SDOI2018]原题识别

    题解: ..感觉挺烦得 而且我都没有注意到树随机这件事情.. 就写个30分的莫队.. #include <bits/stdc++.h> using namespace std; #defi ...

  3. [CF676C]Vasya and String(尺取法,原题)

    题目链接:http://codeforces.com/contest/676/problem/C 原题题解链接:http://www.cnblogs.com/vincentX/p/5405468.ht ...

  4. (各个公司面试原题)在线做了一套CC++综合測试题,也来測一下你的水平吧(二)

    刚才把最后的10道题又看了下.也发上来吧. 以下给出试题.和我对题目的一些理解 前10道题地址 (各个公司面试原题)在线做了一套CC++综合測试题.也来測一下你的水平吧(一) 11.设已经有A,B,C ...

  5. NOIP2016原题终结测试(2017081801)

    NOIP2016还有几道原题没有写掉,今天就一并布置掉. 答案的问题,有部分会先放到NOIP题解中,是单独发布的. 最后会汇总放在答案中,各位不要急. 还有,后期会有原创题测试,这个不急,反正11月才 ...

  6. 2019.3.16 noiac的原题模拟赛

    RT,这太谔谔了,我不承认这是模拟赛 但是虽然是搬了三道题,题目本身也还能看,就这么着吧 (怎么机房里就我一道原题都没做过啊 T1 CF24D Broken Robot 比较简单地列出式子之后,我们发 ...

  7. 浅谈《剑指offer》原题:不使用条件、循环语句求1+2+……+n

    转载自:浅谈<剑指offer>原题:求1+2+--+n 如侵犯您的版权,请联系:windeal12@qq.com <剑指offer>上的一道原题,求1+2+--+n,要求不能使 ...

  8. P1048 采药(洛谷,动态规划递推,01背包原题)

    题目直接放链接 P1048 采药 这题只是01背包+背景故事而已 原题来的 PS:我写了一篇很详细的01背包说明,如果下面ac代码有看不懂的地方可以去看看 对01背包的分析与理解(图文) 下面上ac代 ...

  9. NOIP原题 斗地主(20190804)

    题目描述 牛牛最近迷上了一种叫斗地主的扑克游戏.斗地主是一种使用黑桃.红心.梅花.方片的A到K加上大小王的共54张牌来进行的扑克牌游戏.在斗地主中,牌的大小关 系根据牌的数码表示如下:3<4&l ...

随机推荐

  1. Mybatis的应用2 使用mybits+SpringBoot完成第一个查询的demo(随后加增加,更新,删除)

    首先在mapper下面新建一个mysql.xml mysql.xml <?xml version="1.0" encoding="UTF-8" ?> ...

  2. Mybatis的应用1 Mybatis和logback的应用配置

    首先新建一个module, 然后,在pom文件里面添加一些引用的项. pom.xml <?xml version="1.0" encoding="UTF-8&quo ...

  3. (二叉树 递归 DFS) leetcode 100. Same Tree

    Given two binary trees, write a function to check if they are the same or not. Two binary trees are ...

  4. [APIO2007] 风铃

    题目链接 可能是个树上 DP?指针真好玩 23333. 首先对于所有玩具如果有深度差超过 1 的就是无解(在这里贡献 WA * 3),所以 dfs 一遍记录深度是有必要的…… 然后如果有一个点的两颗子 ...

  5. 利用DOS命令窗口进行Mail通信(二)

    一:SMTP协议(对邮件进行发送) <SP>代表空格,<CRLF>代表回车和换行 SMTP命令格式 说明 ehlo<SP><domain><CRL ...

  6. Subversion配置

    1.下载Apache-Subversion-1.9.7:https://github.com/wangfajun/dev-tools 2.打开idea-->File-->Settings, ...

  7. Linux记录-grafana opentsdb安装

    wget https://dl.grafana.com/oss/release/grafana-6.1.4-1.x86_64.rpm sudo yum localinstall grafana-6.1 ...

  8. css样式关键字(initial,inherit,unset,revert,all)

    关键字 在CSS中,有4个关键字理论上可以应用于任何的CSS属性,它们是initial(初始).inherit(继承).unset(未设置).revert(还原).而all的取值只能是以上这4个关键字 ...

  9. CSS布局 两列布局之单列定宽,单列自适应布局思路

    前言 说起自适应布局方式,单列定宽单列自适应布局是最基本的布局形式.比如斗鱼的直播间,后台管理系统都是常用的 我们将从 float, inline-block, table, absolute, fl ...

  10. 常见的游戏AI技术对比(FSM,HFSM,BT,GOAP,HTN,Utilitay,机器学习)

    最近研究的内容,比较懒还没写文章,先占个坑,过几天补