[Wc2009]shortest

传送门

终于把这题过了，了了我两年前写堵塞的交通一晚上无果的心结

因为是6要注意蛇皮走位啊！！这种-> S

 //Achen
 #include<bits/stdc++.h>
 #define For(i,a,b) for(int i=(a);i<=(b);i++)
 #define Rep(i,a,b) for(int i=(a);i>=(b);i--)
 #define Formylove return 0
 const int N=1e5+;
 typedef long long LL;
 typedef double db;
 using namespace std;
 int n,q,a[][N],tpa[],mx;

 template<typename T> void read(T &x) {
     char ch=getchar(); T f=; x=;
     while(ch!='-'&&(ch<''||ch>'')) ch=getchar();
     if(ch=='-') f=-,ch=getchar();
     for(;ch>=''&&ch<='';ch=getchar()) x=x*+ch-''; x*=f;
 }

 #define inf 1e15
 LL dl[N<<][][],dr[N<<][][],dm[N<<][][];
 void copyit(int x,int y) {
     For(i,,) For(j,,) dl[y][i][j]=dl[x][i][j],dr[y][i][j]=dr[x][i][j],dm[y][i][j]=dm[x][i][j];
 }
 void upd(int y,int ll,int rr) {
     int x=mx+,lc=mx+,rc=mx+;
     copyit(ll,lc); copyit(rr,rc);
     For(i,,) For(j,,) For(k,,) {
         dm[lc][i][j]=min(dm[lc][i][j],dm[lc][i][k]+dl[rc][k][j]);
         dm[rc][i][j]=min(dm[rc][i][j],dm[rc][k][j]+dr[lc][k][i]);
     }
     For(i,,) For(j,,) {
         dm[x][i][j]=inf;
         For(k,,) dm[x][i][j]=min(dm[x][i][j],dm[lc][i][k]+dm[rc][k][j]);
     }
     For(i,,) For(j,,) {
         dl[x][i][j]=dl[lc][i][j];
         dr[x][i][j]=dr[rc][i][j];
         For(k,,) dl[x][i][j]=min(dl[x][i][j],dm[lc][i][k]+dm[lc][j][k]);
         For(k,,) dr[x][i][j]=min(dr[x][i][j],dm[rc][k][i]+dm[rc][k][j]);
     }
     copyit(x,y);
 }

 #define lc (x<<1)
 #define rc ((x<<1)|1)
 #define mid ((l+r)>>1)
 void updone(int x,int l) {
     For(i,,) tpa[i]=(i?tpa[i-]:)+a[i][l];
     For(i,,) For(j,,) {
         LL ds=i>j?(tpa[i]-(j?tpa[j-]:)):(tpa[j]-(i?tpa[i-]:));
         dl[x][i][j]=dr[x][i][j]=dm[x][i][j]=ds;
     }
 }

 void build(int x,int l,int r) {
     if(l==r) { updone(x,l); return ; }
     build(lc,l,mid); build(rc,mid+,r);
     upd(x,lc,rc);
 }

 void change(int x,int l,int r,int pos) {
     if(l==r) { updone(x,l); return; }
     if(pos<=mid) change(lc,l,mid,pos);
     else change(rc,mid+,r,pos);
     upd(x,lc,rc);
 }

 int fl;
 void qry(int x,int l,int r,int ql,int qr,int t) {
     if(ql>qr) return ;
     if(l>=ql&&r<=qr) {
         if(!fl) copyit(x,t),fl=;
         else upd(t,t,x); return ;
     }
     if(ql<=mid) qry(lc,l,mid,ql,qr,t);
     if(qr>mid) qry(rc,mid+,r,ql,qr,t);
 }

 int main() {
     //freopen("2104.in","r",stdin);
     //freopen("2104.out","w",stdout);
     read(n); mx=(n<<);
     For(i,,) For(j,,n) read(a[i][j]);
     build(,,n);
     read(q);
     For(i,,q) {
         int op,x,y,v,sx,sy,tx,ty;
         read(op);
         if(op==) {
             read(x); read(y); read(v);
             a[x-][y]=v; change(,,n,y);
         }
         else {
             read(sx); read(sy); sx--;
             read(tx); read(ty); tx--;
             if(sy>ty) swap(sx,tx),swap(sy,ty);
             fl=; qry(,,n,sy,ty,mx+);
             fl=; qry(,,n,,sy-,mx+);
             fl=; qry(,,n,ty+,n,mx+);
             LL ans=dm[mx+][sx][tx];
             if(sy!=) {
                 For(j,,) For(k,,) dr[mx+][sx][j]=min(dr[mx+][sx][j],dl[mx+][sx][k]+dr[mx+][k][j]-dl[mx+][sx][sx]);
                 For(k,,) ans=min(ans,dl[mx+][sx][sx]+dr[mx+][sx][k]+dm[mx+][k][tx]);
             }
             if(ty!=n) {
                 For(j,,) For(k,,) dl[mx+][j][tx]=min(dl[mx+][j][tx],dr[mx+][k][tx]+dl[mx+][j][k]-dr[mx+][tx][tx]);
                 For(k,,) ans=min(ans,dm[mx+][sx][k]+dl[mx+][k][tx]+dr[mx+][tx][tx]);
             }
             if(sy!=&&ty!=n) {
                 For(k,,) For(l,,)
                     ans=min(ans,dr[mx+][sx][k]+dm[mx+][k][l]+dl[mx+][l][tx]+dl[mx+][sx][sx]+dr[mx+][tx][tx]);
             }
             printf("%lld\n",ans);
         }
     }
     Formylove;
 }