BZOJ 2459:[BeiJing 2011]謎の善人

3177 ワード

セグメントツリーの各ノードメンテナンスd[4][4]は、4つの頂点間の最短ルートを表し、マージ時にFloyedでマージし、時分割3セグメントをクエリーしてマージする.
 
#include<cstdio>

#define N 100010

struct P{int d[4][4];}T[N<<2],tmp;

int n,m,op,x,y,u,v,i,j,k,a[N],b[N],c[N],f[8][8],g[8][8],inf=~0U>>2,flag,ans,d1,d2;

inline void up(int&x,int y){if(x>y)x=y;}

inline void floyed(int n){for(k=0;k<n;k++)for(i=0;i<n;i++)for(j=0;j<n;j++)up(f[i][j],f[i][k]+f[k][j]);}

inline void cal(P&x,int p){

  for(i=0;i<4;i++)for(j=0;j<4;j++)if(i==j)f[i][j]=0;else f[i][j]=inf;

  f[0][1]=f[1][0]=b[p];

  f[2][3]=f[3][2]=b[p+1];

  f[0][2]=f[2][0]=a[p];

  f[1][3]=f[3][1]=c[p];

  floyed(4);

  for(i=0;i<4;i++)for(j=0;j<4;j++)x.d[i][j]=f[i][j];

}

inline void up(P&x,P l,P r){

  for(i=0;i<6;i++)for(j=0;j<6;j++)if(i==j)f[i][j]=0;else f[i][j]=inf;

  for(i=0;i<4;i++)for(j=0;j<4;j++)up(f[i][j],l.d[i][j]),up(f[i+2][j+2],r.d[i][j]);

  floyed(6);

  for(i=0;i<4;i++)for(j=0;j<4;j++)x.d[i][j]=f[i>1?i+2:i][j>1?j+2:j];

}

void build(int x,int a,int b){

  if(a+1==b){cal(T[x],a);return;}

  int mid=(a+b)>>1;

  build(x<<1,a,mid),build(x<<1|1,mid,b),up(T[x],T[x<<1],T[x<<1|1]);

}

void change(int x,int a,int b,int c){

  if(a+1==b){cal(T[x],a);return;}

  int mid=(a+b)>>1;

  if(c<mid)change(x<<1,a,mid,c);else change(x<<1|1,mid,b,c);

  up(T[x],T[x<<1],T[x<<1|1]);

}

void ask(int x,int a,int b,int c,int d){

  if(c<=a&&b<=d){

    if(!flag)flag=1,tmp=T[x];else up(tmp,tmp,T[x]);

    return;

  }

  int mid=(a+b)>>1;

  if(c<mid)ask(x<<1,a,mid,c,d);

  if(d>mid)ask(x<<1|1,mid,b,c,d);

}

int main(){

  scanf("%d",&n);

  for(i=1;i<n;i++)scanf("%d",&a[i]);

  for(i=1;i<=n;i++)scanf("%d",&b[i]);

  for(i=1;i<n;i++)scanf("%d",&c[i]);

  build(1,1,n);

  scanf("%d",&m);

  while(m--){

    scanf("%d%d%d",&op,&x,&y);

    if(!op){

      if(x==0)scanf("%d",&a[y]);

      if(x==1)scanf("%d",&b[y]);

      if(x==2)scanf("%d",&c[y]);

      if(y<n)change(1,1,n,y);

      if(y>1)change(1,1,n,y-1);

    }else{

      if(x==y){puts("0");continue;}

      if(x>y)i=x,x=y,y=i;

      u=(x+1)/2,v=(y+1)/2,d1=x&1^1,d2=y&1^1;

      if(u==v){

        ans=b[u];

        if(u>1)flag=0,ask(1,1,n,1,u),up(ans,tmp.d[2][3]);

        if(u<n)flag=0,ask(1,1,n,u,n),up(ans,tmp.d[0][1]);

      }else{

        for(i=0;i<8;i++)for(j=0;j<8;j++)if(i==j)g[i][j]=0;else g[i][j]=inf;

        flag=0,ask(1,1,n,u,v);

        for(i=0;i<4;i++)for(j=0;j<4;j++)up(g[i+2][j+2],tmp.d[i][j]);

        if(u>1)for(flag=0,ask(1,1,n,1,u),i=0;i<4;i++)for(j=0;j<4;j++)up(g[i][j],tmp.d[i][j]);

        if(v<n)for(flag=0,ask(1,1,n,v,n),i=0;i<4;i++)for(j=0;j<4;j++)up(g[i+4][j+4],tmp.d[i][j]);

        for(i=0;i<8;i++)for(j=0;j<8;j++)f[i][j]=g[i][j];

        floyed(8);

        ans=f[2+d1][4+d2];

      }

      printf("%d
",ans); } } return 0; }