【JZOJ 5430】 图

Description

有一个n个点的无向图，给出m条边，每条边的信息形如x,y,c,r
给出q组询问形如u,v,l,r
接下来解释询问以及边的意义
询问表示，一开始你在点u上，然后按顺序处理编号从l到r的边
对于一条边x,y,c,r，你可以进行两种操作：
如果你当前在x点或者y点上，那么你可以走这条边（从x到y或从y到x）并付出c的代价（当然你也可以不走，看操作2）
如果你不走这条边或者不可以走这条边（即你当前不在x或y上），那么你需要付出r的代价询问如果要从u点开始，按顺序处理完编号从l到r的边之后到达v点的最小代价，如果不能到达v，那么输出-1。
边和点的编号从1开始
对于100%的数据，n<=30,m<=20000,q<=200000

Analysis

首先，由于有l到r这种恶心的限制，我们感受出离线处理会方便些
设个状态(x,y,l,r)表示经过l到r条边，x到y的最小代价，状态复杂度O(n^2m^2)，肯定是不能直接存了，而且要能快速计算。
考虑状态之间的合并，(x,y,l,mid)+(y,z,mid+1,r)=(x,z,l,r)，需要O(n^3)
状态前后两端加入，(x,y,l,l)+(y,z,l+1,r)(前端) (x,y,l,r-1)+(y,z,r,r)(后端)=(x,z,l,r)，需要O(n^2)
无法支持两端删除
如果使用莫队，需要打不带删除的版本，而且据说块大小应该设为$m\sqrt {\frac{1}{q}}$
时间带根号，应该过不了
如果使用分块，（坑）。。。。
如果使用分治，每次处理l,r跨mid的询问，可以弄出l~mid的信息，与mid+1~r的信息再合并（需要枚举到mid时走到哪个点）
算一下复杂度应该是$O(qn+(n^2m+q)logm)$

Code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,b,a) for(int i=b;i>=a;i--)
#define max(x,y) ((x)>(y)?(x):(y))
#define min(x,y) ((x)<(y)?(x):(y))
#define mset(a,x) memset(a,x,sizeof(a))
using namespace std;
typedef long long ll;
char ch;
void read(int &n){n=0;int p=1;for(ch=getchar();ch<'0' || ch>'9';ch=getchar())if(ch=='-') p=-1;for(;'0'<=ch && ch<='9';ch=getchar()) n=n*10+ch-'0';n*=p;}
void writeln(int x){if(x==0) {puts("0");return;}char st[20];int t=0;for(;x;x/=10) st[++t]=(x%10)+'0';while(t) putchar(st[t--]);puts("");}
const int N=33,M=2e4+5,Q=2e5+5,INF=2139062143;
int n,m,q,ans[Q],le[M][N][N],ri[M][N][N];
struct lyd
{
    int x,y,l,r,id;
}b[M],c[Q],cl[Q],cr[Q];
bool pd(int x,int y,lyd e)
{
    return x==e.x && y==e.y || x==e.y && y==e.x;
}
void dfs(int l,int r,int ql,int qr)
{
    if(l>r || ql>qr) return;
    if(l==r)
    {
        fo(i,ql,qr) ans[c[i].id]=pd(c[i].x,c[i].y,b[l])?b[l].l:(c[i].x==c[i].y?b[l].r:INF);
        return;
    }
    int mid=l+r>>1;
    mset(le[mid+1],127);
    fo(i,1,n) le[mid+1][i][i]=0;
    mset(ri[mid],127);
    fo(i,1,n) ri[mid][i][i]=0;
    fd(k,mid,l)
    {
        mset(le[k],127);
        fo(i,1,n)
            if(i==b[k].x || i==b[k].y)
                fo(j,1,n) le[k][i][j]=min(b[k].r+le[k+1][i][j],b[k].l+le[k+1][b[k].x+b[k].y-i][j]);
            else
                fo(j,1,n) le[k][i][j]=b[k].r+le[k+1][i][j];
    }
    fo(k,mid+1,r)
    {
        mset(ri[k],127);
        fo(i,1,n)
            fo(j,1,n)
                if(j==b[k].x || j==b[k].y)
                    ri[k][i][j]=min(b[k].r+ri[k-1][i][j],b[k].l+ri[k-1][i][b[k].x+b[k].y-j]);
                else ri[k][i][j]=b[k].r+ri[k-1][i][j];
    }
    int nl=0,nr=0;
    fo(k,ql,qr)
        if(c[k].r<mid) cl[++nl]=c[k];
        else
        if(c[k].l>mid) cr[++nr]=c[k];
        else
        if(c[k].l<=mid && mid<=c[k].r)
        {
            int res=INF;
            fo(i,1,n)
                res=min(res,(ll)le[c[k].l][c[k].x][i]+ri[c[k].r][i][c[k].y]);
            ans[c[k].id]=res;
        }
    fo(i,ql,ql+nl-1) c[i]=cl[i-ql+1];
    fo(i,qr-nr+1,qr) c[i]=cr[i-(qr-nr)];
    dfs(l,mid-1,ql,ql+nl-1);
    dfs(mid+1,r,qr-nr+1,qr);
}
int main()
{
    freopen("graph.in","r",stdin);
    freopen("graph.out","w",stdout);
    read(n),read(m),read(q);
    fo(i,1,m) read(b[i].x),read(b[i].y),read(b[i].l),read(b[i].r);
    fo(i,1,q) read(c[i].x),read(c[i].y),read(c[i].l),read(c[i].r),c[i].id=i;
    dfs(1,m,1,q);
    fo(i,1,q)
        if(ans[i]==INF) puts("-1");else writeln(ans[i]);
    return 0;
}