UOJ 55 【WC2014】紫荆花之恋—

点分治。在点分树上每个点上用 splay 维护管辖的点的情况。做几次就重构点分树。TLE。只能过 20 分。

#include<cstdio>

#include<cstring>

#include<algorithm>

#define ll long long

using namespace std;

int rdn()

{

  int ret=;bool fx=;char ch=getchar();

  while(ch>''||ch<''){if(ch=='-')fx=;ch=getchar();}

  while(ch>=''&&ch<='')ret=ret*+ch-'',ch=getchar();

  return fx?ret:-ret;

}

int Mn(int a,int b){return a<b?a:b;}

int Mx(int a,int b){return a>b?a:b;}

const int N=1e5+,base=,K=base+,M=N*K,mod=1e9;//K not log!!!//K=base+log

int n,hd[N],xnt,to[N<<],nxt[N<<],w[N<<],r[N],sz[N],mn,Rt; ll ans;

int rt[][N],tot,c[M][],vl[M],siz[M],fa[M],pre[N][K],dis[N][K],dep[N];

bool vis[N];

void add(int x,int y,int z){to[++xnt]=y;nxt[xnt]=hd[x];hd[x]=xnt;w[xnt]=z;}

void pshp(int cr){siz[cr]=siz[c[cr][]]+siz[c[cr][]]+;}

void rotate(int x)

{

  int y=fa[x],z=fa[y],d=(x==c[y][]);

  if(z)c[z][y==c[z][]]=x;

  fa[x]=z;

  fa[c[x][!d]]=y; fa[y]=x;

  c[y][d]=c[x][!d]; c[x][!d]=y;

  pshp(y);pshp(x);

}

void splay(int x)

{

  while(fa[x])

    {

      int y=fa[x],z=fa[y];

      if(z)

    ((x==c[y][])^(y==c[z][]))?rotate(x):rotate(y);

      rotate(x);

    }

}

void inst(int &cr,int v,int pr)

{

  if(!cr){cr=++tot;c[cr][]=c[cr][]=;fa[cr]=pr;vl[cr]=v;siz[cr]=;return;}

  siz[cr]++;

  if(v<=vl[cr])inst(c[cr][],v,cr);

  else inst(c[cr][],v,cr);

}

void ins(int x,int v,bool fx)

{

  inst(rt[fx][x],v,);splay(tot);rt[fx][x]=tot;

}

int fnd(int cr,int v)

{

  if(!cr)return -;

  int ret=(vl[cr]>=v?fnd(c[cr][],v):fnd(c[cr][],v));

  if(ret==-&&vl[cr]>=v)ret=cr;

  return ret;

}

void dfsx(int cr)

{

  int ls=c[cr][],rs=c[cr][];

  printf("  cr=%d(%d %d) c0=%d(%d %d) c1=%d(%d %d)\n",

     cr,vl[cr],siz[cr],ls,vl[ls],siz[ls],rs,vl[rs],siz[rs]);

  if(ls)dfsx(ls); if(rs)dfsx(rs);

}

int qry(int x,int v,bool fx)

{

  int cr=fnd(rt[fx][x],v);

  if(cr==-)return ; splay(cr); rt[fx][x]=cr;

  return siz[c[cr][]]+;

}

void getrt(int cr,int fa,int s)

{

  sz[cr]=;int mx=;

  for(int i=hd[cr],v;i;i=nxt[i])

    if(!vis[v=to[i]]&&v!=fa)

      getrt(v,cr,s),sz[cr]+=sz[v],mx=Mx(mx,sz[v]);

  mx=Mx(mx,s-sz[cr]); if(mx<mn)mn=mx,Rt=cr;

}

void dfs(int cr,int fa,int d)

{

  pre[cr][++dep[cr]]=Rt;dis[cr][dep[cr]]=d;

  ins(Rt,r[cr]-d,);if(dep[cr]>)ins(Rt,r[cr]-dis[cr][dep[cr]-],);

  for(int i=hd[cr],v;i;i=nxt[i])

    if(!vis[v=to[i]]&&v!=fa)dfs(v,cr,d+w[i]);

}

void solve(int cr,int s)

{

  vis[cr]=;dfs(cr,,);

  for(int i=hd[cr],v;i;i=nxt[i])

    {

      if(vis[v=to[i]])continue;

      int ts=(sz[v]<sz[cr]?sz[v]:s-sz[cr]);

      mn=N;getrt(v,cr,ts);solve(Rt,ts);

    }

}

void cz(int s)

{

  tot=;

  for(int i=;i<=s;i++)vis[i]=;for(int i=;i<=s;i++)dep[i]=;

  for(int i=;i<=s;i++)rt[][i]=rt[][i]=;

  mn=N;getrt(,,s);solve(Rt,s);

}

int main()

{

  int bh=rdn();n=rdn();int pr=rdn()%mod,tw=rdn();r[]=rdn();puts("");

  ins(,r[],);dep[]=;dis[][]=;pre[][]=;

  for(int i=,T=;i<=n;i++,T++)

    {

      if(T==base)cz(i-),T=;

      pr=rdn()^(ans%mod);tw=rdn();r[i]=rdn();

      add(pr,i,tw);add(i,pr,tw);

      for(int j=;j<=dep[pr];j++)

    {

      pre[i][j]=pre[pr][j];dis[i][j]=dis[pr][j]+tw;

      ins(pre[i][j],r[i]-dis[i][j],);

      if(j>)ins(pre[i][j],r[i]-dis[i][j-],);

    }

      dep[i]=dep[pr]+;pre[i][dep[i]]=i;

      ins(i,r[i],);if(dep[i]>)ins(i,r[i]-dis[i][dep[i]-],);

      for(int j=;j<dep[i];j++)

    {

      int u=pre[i][j],v=pre[i][j+],d=dis[i][j]-r[i];

      ans+=qry(u,d,);

      ans-=qry(v,d,);

    }

      printf("%lld\n",ans);

    }

  return ;

}

写对拍的话，可以这样生成，就是在 mk.cpp 里调用写好的程序 a.cpp 得到 ans 。可以把 a.cpp 删掉 main 函数改成 a.hpp 。

#include<cstdio>

#include<cstring>

#include<algorithm>

#define ll long long

using namespace std;

int rdn()

{

  int ret=;bool fx=;char ch=getchar();

  while(ch>''||ch<''){if(ch=='-')fx=;ch=getchar();}

  while(ch>=''&&ch<='')ret=ret*+ch-'',ch=getchar();

  return fx?ret:-ret;

}

int Mn(int a,int b){return a<b?a:b;}

int Mx(int a,int b){return a>b?a:b;}

const int N=,base=,K=base+,M=N*K,mod=1e9;//K not log!!!//K=base+log

int n,hd[N],xnt,to[N<<],nxt[N<<],w[N<<],r[N],sz[N],mn,Rt; ll ans;

int rt[][N],tot,c[M][],vl[M],siz[M],fa[M],pre[N][K],dis[N][K],dep[N];

bool vis[N];

void add(int x,int y,int z){to[++xnt]=y;nxt[xnt]=hd[x];hd[x]=xnt;w[xnt]=z;}

void pshp(int cr){siz[cr]=siz[c[cr][]]+siz[c[cr][]]+;}

void rotate(int x)

{

  int y=fa[x],z=fa[y],d=(x==c[y][]);

  if(z)c[z][y==c[z][]]=x;

  fa[x]=z;

  fa[c[x][!d]]=y; fa[y]=x;

  c[y][d]=c[x][!d]; c[x][!d]=y;

  pshp(y);pshp(x);

}

void splay(int x)

{

  while(fa[x])

    {

      int y=fa[x],z=fa[y];

      if(z)

    ((x==c[y][])^(y==c[z][]))?rotate(x):rotate(y);

      rotate(x);

    }

}

void inst(int &cr,int v,int pr)

{

  if(!cr){cr=++tot;c[cr][]=c[cr][]=;fa[cr]=pr;vl[cr]=v;siz[cr]=;return;}

  siz[cr]++;

  if(v<=vl[cr])inst(c[cr][],v,cr);

  else inst(c[cr][],v,cr);

}

void ins(int x,int v,bool fx)

{

  inst(rt[fx][x],v,);splay(tot);rt[fx][x]=tot;

}

int fnd(int cr,int v)

{

  if(!cr)return -;

  int ret=(vl[cr]>=v?fnd(c[cr][],v):fnd(c[cr][],v));

  if(ret==-&&vl[cr]>=v)ret=cr;

  return ret;

}

void dfsx(int cr)

{

  int ls=c[cr][],rs=c[cr][];

  printf("  cr=%d(%d %d) c0=%d(%d %d) c1=%d(%d %d)\n",

     cr,vl[cr],siz[cr],ls,vl[ls],siz[ls],rs,vl[rs],siz[rs]);

  if(ls)dfsx(ls); if(rs)dfsx(rs);

}

int qry(int x,int v,bool fx)

{

  int cr=fnd(rt[fx][x],v);

  if(cr==-)return ; splay(cr); rt[fx][x]=cr;

  return siz[c[cr][]]+;

}

void getrt(int cr,int fa,int s)

{

  sz[cr]=;int mx=;

  for(int i=hd[cr],v;i;i=nxt[i])

    if(!vis[v=to[i]]&&v!=fa)

      getrt(v,cr,s),sz[cr]+=sz[v],mx=Mx(mx,sz[v]);

  mx=Mx(mx,s-sz[cr]); if(mx<mn)mn=mx,Rt=cr;

}

void dfs(int cr,int fa,int d)

{

  pre[cr][++dep[cr]]=Rt;dis[cr][dep[cr]]=d;

  ins(Rt,r[cr]-d,);if(dep[cr]>)ins(Rt,r[cr]-dis[cr][dep[cr]-],);

  for(int i=hd[cr],v;i;i=nxt[i])

    if(!vis[v=to[i]]&&v!=fa)dfs(v,cr,d+w[i]);

}

void solve(int cr,int s)

{

  vis[cr]=;dfs(cr,,);

  for(int i=hd[cr],v;i;i=nxt[i])

    {

      if(vis[v=to[i]])continue;

      int ts=(sz[v]<sz[cr]?sz[v]:s-sz[cr]);

      mn=N;getrt(v,cr,ts);solve(Rt,ts);

    }

}

void cz(int s)

{

  tot=;

  for(int i=;i<=s;i++)vis[i]=;for(int i=;i<=s;i++)dep[i]=;

  for(int i=;i<=s;i++)rt[][i]=rt[][i]=;

  mn=N;getrt(,,s);solve(Rt,s);

}

a.hpp

#include<cstdio>

#include<cstring>

#include<algorithm>

#include<ctime>

#define ll long long

#include "a.hpp"

using namespace std;

const int M1=,M2=;

int main()

{

  srand(time());srand(rand());

  n=N;//n=rand()%N+1;

  printf("%d\n%d\n",rand()%+,n);

  r[]=rand()%M2+;

  printf("0 0 %d\n",r[]);

  ins(,r[],);dep[]=;dis[][]=;pre[][]=;

  int pr,tw;

  for(int i=,T=;i<=n;i++,T++)

    {

      if(T==base)cz(i-),T=;

      pr=rand()%(i-)+;

      tw=rand()%M1+;r[i]=rand()%M2+;

      printf("%lld %d %d\n",pr^(ans%mod),tw,r[i]);

      add(pr,i,tw);add(i,pr,tw);

      for(int j=;j<=dep[pr];j++)

    {

      pre[i][j]=pre[pr][j];dis[i][j]=dis[pr][j]+tw;

      ins(pre[i][j],r[i]-dis[i][j],);

      if(j>)ins(pre[i][j],r[i]-dis[i][j-],);

    }

      dep[i]=dep[pr]+;pre[i][dep[i]]=i;

      ins(i,r[i],);if(dep[i]>)ins(i,r[i]-dis[i][dep[i]-],);

      for(int j=;j<dep[i];j++)

    {

      int u=pre[i][j],v=pre[i][j+],d=dis[i][j]-r[i];

      ans+=qry(u,d,);

      ans-=qry(v,d,);

    }

    }

  return ;

}

mk.cpp

尝试改一改。设一个 base ，重构点分树只重构不平衡的那部分。再把 splay 的一些操作写成非递归的，就能过 85 分。主要是 r <= 10 的数据过不了。

#include<cstdio>

#include<cstring>

#include<algorithm>

#include<vector>

#define ll long long

#define db double

#define pb push_back

#define ls c[cr][0]

#define rs c[cr][1]

using namespace std;

int rdn()

{

  int ret=;bool fx=;char ch=getchar();

  while(ch>''||ch<''){if(ch=='-')fx=;ch=getchar();}

  while(ch>=''&&ch<='')ret=ret*+ch-'',ch=getchar();

  return fx?ret:-ret;

}

int Mn(int a,int b){return a<b?a:b;}

int Mx(int a,int b){return a>b?a:b;}

const int N=1e5+,K=,M=N*K,mod=1e9; const db base=0.77;

int n,hd[N],xnt,to[N<<],nxt[N<<],w[N<<],r[N],sz[N],mn,Rt; ll ans;

int rt[][N],tot,c[M][],vl[M],siz[M],fa[M],pre[N][K],dis[N][K],dep[N];

int delpl[M],top; bool vis[N];

vector<int> vt[N];

void add(int x,int y,int z){to[++xnt]=y;nxt[xnt]=hd[x];hd[x]=xnt;w[xnt]=z;}

void pshp(int cr){siz[cr]=siz[c[cr][]]+siz[c[cr][]]+;}

void rotate(int x)

{

  int y=fa[x],z=fa[y],d=(x==c[y][]);

  if(z)c[z][y==c[z][]]=x;

  fa[x]=z;

  fa[c[x][!d]]=y; fa[y]=x;

  c[y][d]=c[x][!d]; c[x][!d]=y;

  pshp(y);pshp(x);

}

void splay(int x)

{

  for(int y,z;fa[x];rotate(x))

    {

      y=fa[x];z=fa[y];

      if(z) ((x==c[y][])^(y==c[z][]))?rotate(x):rotate(y);

    }

}

int nwnd()

{

  int ret=(top?delpl[top--]:++tot);

  c[ret][]=c[ret][]=;siz[ret]=;return ret;

}

void ins(int x,int v,bool fx,bool nd)

{

  if(!rt[fx][x]){rt[fx][x]=nwnd();vl[rt[fx][x]]=v;fa[rt[fx][x]]=;return;}

  int cr=rt[fx][x],pr=;

  while()

    {

      siz[cr]++;pr=cr;

      if(v<=vl[cr])

    {

      if(!ls){ls=nwnd();fa[ls]=cr;vl[ls]=v;break;}

      cr=ls;

    }

      else

    {

      if(!rs){rs=nwnd();fa[rs]=cr;vl[rs]=v;break;}

      cr=rs;

    }

    }

  if(nd)splay(cr),rt[fx][x]=cr;

}

void dfs_del(int cr)

{

  delpl[++top]=cr;

  if(ls)dfs_del(ls);if(rs)dfs_del(rs);

}

void del(int cr)

{

  dfs_del(rt[][cr]);if(dep[cr]>)dfs_del(rt[][cr]);

  rt[][cr]=rt[][cr]=;

}

int fnd(int cr,int v)

{

  int ret=-;

  while(cr)

    {

      if(vl[cr]>=v)ret=cr,cr=ls;

      else cr=rs;

    }

  return ret;

}

int qry(int x,int v,bool fx)

{

  int cr=fnd(rt[fx][x],v);

  if(cr==-)return ; splay(cr); rt[fx][x]=cr;

  return siz[rs]+;

}

void getrt(int cr,int fa,int s)

{

  sz[cr]=;int mx=;

  for(int i=hd[cr],v;i;i=nxt[i])

    if(!vis[v=to[i]]&&v!=fa)

      getrt(v,cr,s),sz[cr]+=sz[v],mx=Mx(mx,sz[v]);

  mx=Mx(mx,s-sz[cr]); if(mx<mn)mn=mx,Rt=cr;

}

void dfs(int cr,int fa,int d)

{

  pre[cr][++dep[cr]]=Rt;dis[cr][dep[cr]]=d;

  vt[Rt].pb(cr);

  ins(Rt,r[cr]-d,,);if(dep[cr]>)ins(Rt,r[cr]-dis[cr][dep[cr]-],,);

  for(int i=hd[cr],v;i;i=nxt[i])

    if(!vis[v=to[i]]&&v!=fa)dfs(v,cr,d+w[i]);

}

void solve(int cr,int s)

{

  vis[cr]=;vt[cr].clear();rt[][cr]=rt[][cr]=;

  dfs(cr,,);

  for(int i=hd[cr],v;i;i=nxt[i])

    {

      if(vis[v=to[i]])continue;

      int ts=(sz[v]<sz[cr]?sz[v]:s-sz[cr]);

      mn=N;getrt(v,cr,ts);solve(Rt,ts);

    }

}

void cz(int cr)

{

  int s=vt[cr].size();

  for(int i=,d,t=dep[cr]-;i<s;i++)

    {d=vt[cr][i];vis[d]=;del(d);dep[d]=t;}//t!!

  mn=N;getrt(cr,,s);solve(Rt,s);

}

int main()

{

  int bh=rdn();n=rdn();int pr=rdn()%mod,tw=rdn();r[]=rdn();puts("");

  ins(,r[],,);dep[]=;dis[][]=;pre[][]=;vt[].pb();vis[]=;

  bh=;

  for(int i=;i<=n;i++)

    {

      pr=rdn()^(ans%mod);tw=rdn();r[i]=rdn();

      add(pr,i,tw);add(i,pr,tw);

      for(int j=,u,v;j<=dep[pr];j++,v=u)

    {

      u=pre[i][j]=pre[pr][j];dis[i][j]=dis[pr][j]+tw;

      ins(u,r[i]-dis[i][j],,);

      if(j>)ins(u,r[i]-dis[i][j-],,);

      vt[u].pb(i);

      if(!bh&&j>&&vt[u].size()>vt[v].size()*base)bh=v;

    }

      dep[i]=dep[pr]+;pre[i][dep[i]]=i; vis[i]=;

      ins(i,r[i],,);if(dep[i]>)ins(i,r[i]-dis[i][dep[i]-],,);

      vt[i].pb(i);

      for(int j=;j<dep[i];j++)

    {

      int u=pre[i][j],v=pre[i][j+],d=dis[i][j]-r[i];

      ans+=qry(u,d,);

      ans-=qry(v,d,);

    }

      printf("%lld\n",ans);

      if(bh)cz(bh),bh=;

    }

  return ;

}

然后学习一下 SBT ，还是过不了，有 80 分。于是不想管了。

#include<cstdio>

#include<cstring>

#include<algorithm>

#include<vector>

#define ll long long

#define db double

#define pb push_back

#define ls c[cr][0]

#define rs c[cr][1]

using namespace std;

int rdn()

{

  int ret=;bool fx=;char ch=getchar();

  while(ch>''||ch<''){if(ch=='-')fx=;ch=getchar();}

  while(ch>=''&&ch<='')ret=ret*+ch-'',ch=getchar();

  return fx?ret:-ret;

}

int Mn(int a,int b){return a<b?a:b;}

int Mx(int a,int b){return a>b?a:b;}

const int N=1e5+,K=,M=N*K,mod=1e9; const db base=0.77;

int n,hd[N],xnt,to[N<<],nxt[N<<],w[N<<],r[N],sz[N],mn,Rt; ll ans;

int rt[][N],tot,c[M][],vl[M],siz[M],fa[M],pre[N][K],dis[N][K],dep[N];

int delpl[M],top; bool vis[N];

vector<int> vt[N];

void add(int x,int y,int z){to[++xnt]=y;nxt[xnt]=hd[x];hd[x]=xnt;w[xnt]=z;}

void pshp(int cr){siz[cr]=siz[c[cr][]]+siz[c[cr][]]+;}

void rotate(int &cr,int d)

{

  int y=c[cr][!d];

  c[cr][!d]=c[y][d];  c[y][d]=cr;

  pshp(cr); pshp(y);  cr=y;

}

void maintain(int cr,int d)

{

  if(siz[c[c[cr][d]][d]]>siz[cr])rotate(cr,!d);

  else if(siz[c[c[cr][d]][!d]]>siz[cr])rotate(c[cr][d],d),rotate(cr,!d);

  else return;///

  maintain(ls,);  maintain(rs,);

  maintain(cr,);  maintain(cr,);

}

int nwnd()

{

  int cr=(top?delpl[top--]:++tot);

  siz[cr]=;return cr;

}

void ins(int &cr,int v,int pr=)

{

  if(!cr){cr=nwnd();fa[cr]=pr;vl[cr]=v;return;}

  siz[cr]++;

  bool d=(v>vl[cr]); ins(c[cr][d],v,cr); maintain(cr,d);

}

void del(int &cr)

{

  delpl[++top]=cr;

  if(ls)del(ls);if(rs)del(rs); cr=;

}

int qry(int cr,int v)

{

  if(!cr)return ;

  if(vl[cr]>=v)return siz[rs]++qry(ls,v);

  return qry(rs,v);

}

void getrt(int cr,int fa,int s)

{

  sz[cr]=;int mx=;

  for(int i=hd[cr],v;i;i=nxt[i])

    if(!vis[v=to[i]]&&v!=fa)

      getrt(v,cr,s),sz[cr]+=sz[v],mx=Mx(mx,sz[v]);

  mx=Mx(mx,s-sz[cr]); if(mx<mn)mn=mx,Rt=cr;

}

void dfs(int cr,int fa,int d,int t)

{

  pre[cr][++dep[cr]]=Rt;dis[cr][dep[cr]]=d;

  vt[Rt].pb(cr);

  ins(rt[][Rt],r[cr]-d);if(t)ins(rt[][Rt],r[cr]-dis[cr][t]);

  for(int i=hd[cr],v;i;i=nxt[i])

    if(!vis[v=to[i]]&&v!=fa)dfs(v,cr,d+w[i],t);

}

void solve(int cr,int s,int t)//t:dep[Rt]-1

{

  vis[cr]=;vt[cr].clear();rt[][cr]=rt[][cr]=;

  dfs(cr,,,t);

  for(int i=hd[cr],v;i;i=nxt[i])

    {

      if(vis[v=to[i]])continue;

      int ts=(sz[v]<sz[cr]?sz[v]:s-sz[cr]);

      mn=N;getrt(v,cr,ts);solve(Rt,ts,t+);

    }

}

void cz(int cr)

{

  int s=vt[cr].size(),t=dep[cr]-;//t!!

  for(int i=,d;i<s;i++)

    {d=vt[cr][i];vis[d]=;del(rt[][d]);if(dep[d]>)del(rt[][d]);dep[d]=t;}

  mn=N;getrt(cr,,s);solve(Rt,s,t);

}

int main()

{

  int bh=rdn();n=rdn();int pr=rdn()%mod,tw=rdn();r[]=rdn();puts("");

  ins(rt[][],r[]);dep[]=;dis[][]=;pre[][]=;vt[].pb();vis[]=;

  bh=;

  for(int i=;i<=n;i++)

    {

      pr=rdn()^(ans%mod);tw=rdn();r[i]=rdn();

      add(pr,i,tw);add(i,pr,tw);

      vis[i]=; dep[i]=dep[pr]+; int lm=dep[i];//vis[]

      for(int j=,u,v;j<=lm;j++,v=u)

    {

      u=pre[i][j]=(j==lm?i:pre[pr][j]);

      dis[i][j]=(j==lm?:dis[pr][j]+tw);

      ins(rt[][u],r[i]-dis[i][j]);

      if(j>)ins(rt[][u],r[i]-dis[i][j-]);

      vt[u].pb(i);

      if(!bh&&j>&&vt[u].size()>vt[v].size()*base)bh=v;

    }

      for(int j=,u=pre[i][j],v,d;j<lm;j++,u=v)

    {

      v=pre[i][j+];d=dis[i][j]-r[i];

      ans+=qry(rt[][u],d);

      ans-=qry(rt[][v],d);

    }

      printf("%lld\n",ans);

      if(bh)cz(bh),bh=;

    }

  return ;

}

UOJ 55 【WC2014】紫荆花之恋——点分治+平衡树的更多相关文章

UOJ#55 [WC2014]紫荆花之恋
题目描述强强和萌萌是一对好朋友.有一天他们在外面闲逛,突然看到前方有一棵紫荆树.这已经是紫荆花飞舞的季节了,无数的花瓣以肉眼可见的速度从紫荆树上长了出来. 仔细看看的话,这个大树实际上是一个带权树. ...
bzoj 3435: [Wc2014]紫荆花之恋替罪羊树维护点分治 && AC400
3435: [Wc2014]紫荆花之恋 Time Limit: 240 Sec Memory Limit: 512 MBSubmit: 159 Solved: 40[Submit][Status] ...
【BZOJ3435】[Wc2014]紫荆花之恋替罪点分树+SBT
[BZOJ3435][Wc2014]紫荆花之恋 Description 强强和萌萌是一对好朋友.有一天他们在外面闲逛,突然看到前方有一棵紫荆树.这已经是紫荆花飞舞的季节了,无数的花瓣以肉眼可见的速度从 ...
BZOJ 3435: [Wc2014]紫荆花之恋
二次联通门 : BZOJ 3435: [Wc2014]紫荆花之恋二次联通门 : luogu P3920 [WC2014]紫荆花之恋 /* luogu P3920 [WC2014]紫荆花之恋怀疑人生 ...
luogu P3920 [WC2014]紫荆花之恋
LINK:紫荆花之恋每次动态加入一个节点统计有多少个节点和当前节点的距离小于他们的权值和. 显然我们不能n^2暴力. 考虑一个简化版的问题树已经给出每次求某个节点和其他节点的贡献. 不难想到 ...
BZOJ3435[Wc2014]紫荆花之恋——动态点分治(替罪羊式点分树套替罪羊树)
题目描述强强和萌萌是一对好朋友.有一天他们在外面闲逛,突然看到前方有一棵紫荆树.这已经是紫荆花飞舞的季节了,无数的花瓣以肉眼可见的速度从紫荆树上长了出来.仔细看看的话,这个大树实际上是一个带权树.每 ...
[WC2014]紫荆花之恋(动态点分治+替罪羊思想)
题目描述强强和萌萌是一对好朋友.有一天他们在外面闲逛,突然看到前方有一棵紫荆树.这已经是紫荆花飞舞的季节了,无数的花瓣以肉眼可见的速度从紫荆树上长了出来.仔细看看的话,这个大树实际上是一个带权树.每 ...
luoguP3920 [WC2014]紫荆花之恋动态点分治 + 替罪羊树
意外的好写..... 考虑点分 \(dis(i, j) \leq r_i + r_j\) 对于过分治中心一点\(u\),有 \(dis(i, u) - r_i = dis(j, u) + r_j\) ...
BZOJ 3435 / Luogu 3920 [WC2014]紫荆花之恋 (替罪羊树动态点分治套 Treap)
题意略分析引用PoPoQQQ的话吾辈有生之年终于把这道题切了...QAQ (蒟蒻狂笑) Orz PoPoQQQ,我又抄PoPoQQQ的题解了 - 突然发现有旋Treap没那么难写学习了一波C ...

随机推荐

source insight 4.0的基本使用方法（转）
源:source insight 4.0的基本使用方法 source insight 4设置
Python Web学习笔记之为什么设计GIL
GIL(global interpreter lock),全局解释器锁,是很多编程语言实现中都具有的特性,由于它的存在,解释器无法实现真正的并发.它也是 Python 中经常讨论的话题之一. Pyth ...
linux及安全第七周总结——20135227黄晓妍
实验部分首先clone最新的menu 我们可以看到,test.c里多了一个exec的功能,它的代码和fork基本一致,多了一项加载hello rootfs也有一些变化执行一下exec 让我们启动一 ...
net.sf.json和 com.fasterxml.jackson中对象转json的区别
近期做项目的时候,发现使用net.sf.json包中的JSONObject或JSONArray将对象转为json数据结构存在一个坑.当对String类型的属性赋值为null情况下,转为json结构为& ...
POJ 1185 炮兵阵地（状压DP）题解
思路:和上一篇思路一样,但是这里要求最大能排几个,这里要开三维,记录上次和上上次的状态,再一一判定,状态转移方程为 dp[i][j][k] = max(dp[i][j][k],dp[i - 1][k] ...
P4 PI库安装
repo:Github PI P4Runtime - a control plane framework and tools for the P4 programming language 这个是P4 ...
Gym 101147J Whistle's New Car（dfs）
https://vjudge.net/problem/Gym-101147J 题意: 有n个城市,每个城市有一个权值,表示在这个城市的加油站可以加多少油. 现在要计算每个城市i,有多少个城市j可以到达 ...
JQuery 实现下拉列表选中
html代码如下: <select id="category" name="category"> <option value="&q ...
MAC OSX Docker安装教程2018
1.首先打开Docker商店 Docker商店地址:https://store.docker.com/editions/community/docker-ce-desktop-mac 2.打开后界面如 ...
python 进制转换
print hex(),hex(-) #转换成十六进制 print oct(),oct(-) #转换成八进制 print bin(),bin(-) #转换成二进制 print int("字面 ...

UOJ 55 【WC2014】紫荆花之恋——点分治+平衡树

UOJ 55 【WC2014】紫荆花之恋——点分治+平衡树的更多相关文章

随机推荐

热门专题