jzoj5430 【NOIP2017提高A组集训10.27】图

题目

Time Limits: 2000 ms Memory Limits: 524288 KB

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开始

Input

第一行三个数表示n,m,q
接下来m行，每行四个整数x,y,c,r描述一条边
接下来q行，每行四个整数u,v,l,r描述一组询问

Output

输出共q行，每行一个数表示对应询问的答案

Sample Input

输入1：
5 5 3
1 4 4 5
4 1 6 1
2 1 2 9
2 5 1 0
1 5 2 5
2 2 2 4
5 4 5 5
1 5 2 5
输入2：
4 8 6
2 4 5 8
2 4 4 8
2 3 6 4
1 4 5 0
2 4 10 10
1 3 5 2
3 2 2 9
3 4 1 1
3 2 1 5
3 1 2 2
1 1 1 7
2 3 2 4
3 3 1 7
1 2 2 5

Sample Output

输出1：
10
-1
9
输出2：
32
-1
41
14
36
27

Data Constraint

对于20%的数据，n<=10,m<=1000,q<=1000
对于40%的数据，n<=10,m<=10000,q<=30000
对于60%的数据，n<=20,m<=10000,q<=30000
对于80%的数据，n<=25,m<=10000,q<=200000
对于100%的数据，n<=30,m<=20000,q<=200000

题解

不妨把询问按照位置分类，然后用分治解决此题
假设当前的分治区间是(l,r)，那么我们现在要计算所有左端点在[l,mid]，右端点[mid,r]的所有区间
设lef[i][j][k]表示按照题目意思依次处理完i到mid的边，从j走到了k的最小代价
rig[i][j][k]表示依次处理完mid+1到i的边，从j走到了k的最小代价

直接转移就可以了

时间复杂度是m log m *n^2的
卡卡常数吧。。。
把min改成？语句就可以了

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<cstdio>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fo1(i,b,a) for(i=b;i>=a;i--)
#define ll long long
#define min(x,y) (x<y?x:y)
using namespace std;

const int maxn=33,maxm=2e4+5,maxq=2e5+5;

ll lef[maxm][maxn][maxn],rig[maxm][maxn][maxn];
ll as[maxq][5],bi[maxm][5],an[maxq];
ll dui[maxq],fi[maxm*20],qc[maxm],ne[maxq*20];
ll i,j,k,l,r,u,m,n,q,x,y,v,vv,now,na,mid;
char ch;

void add(int x,int y){
    if (fi[x]==0) fi[x]=++now; else ne[qc[x]]=++now;
    dui[now]=y; qc[x]=now;
}
int read(){
    int x=0;
    ch=getchar();
    while (ch<'0' || ch>'9') ch=getchar();
    while (ch>='0' && ch<='9'){
        x=x*10+ch-48;
        ch=getchar();
    }
    return x;
}
void dfs(int v,int l,int r){
/*  fo(i,l-1,r+1) fo(j,1,n) fo(k,1,n){
        rig[i][j][k]=na;
        lef[i][j][k]=na;
    }*/
    int mid=(l+r)/2;
    fo(j,1,n) fo(k,1,n){
        rig[mid][j][k]=na;
        lef[mid+1][j][k]=na;
    }
    fo(i,1,n){
        lef[mid+1][i][i]=0;
        rig[mid][i][i]=0;
    }
    fo1(i,mid,l){
        fo(j,1,n){
            fo(k,1,n){
                lef[i][j][k]=na;
                lef[i][j][k]=min(lef[i][j][k],lef[i+1][j][k]+bi[i][4]);
                if (j==bi[i][1])
                lef[i][j][k]=min(lef[i][j][k],lef[i+1][bi[i][2]][k]+bi[i][3]);
                if (j==bi[i][2])
                lef[i][j][k]=min(lef[i][j][k],lef[i+1][bi[i][1]][k]+bi[i][3]);
            }
        }
    }
    fo(i,mid+1,r){
        fo(j,1,n){
            fo(k,1,n){
                rig[i][j][k]=na;
                rig[i][j][k]=min(rig[i][j][k],rig[i-1][j][k]+bi[i][4]);
                if (k==bi[i][1]) 
                rig[i][j][k]=min(rig[i][j][k],rig[i-1][j][bi[i][2]]+bi[i][3]);
                if (k==bi[i][2])
                rig[i][j][k]=min(rig[i][j][k],rig[i-1][j][bi[i][1]]+bi[i][3]);
            }
        }
    }
    i=fi[v];
    while (i){
        x=as[dui[i]][3]; y=as[dui[i]][4]; u=as[dui[i]][1]; vv=as[dui[i]][2];
        an[dui[i]]=na;
        fo(j,1,n) an[dui[i]]=min(an[dui[i]],lef[x][u][j]+rig[y][j][vv]);
        if (an[dui[i]]==na) an[dui[i]]=-1;
        i=ne[i];
    }
    if (l!=r){
        dfs(v*2,l,mid);
        dfs(v*2+1,mid+1,r);
    }
}
int main(){
    freopen("graph.in","r",stdin);
    freopen("graph.out","w",stdout);
    n=read(); m=read(); q=read();
    fo(i,1,m){
        bi[i][1]=read(); bi[i][2]=read(); bi[i][3]=read(); bi[i][4]=read();
    }
    fo(i,1,q){
        as[i][1]=read(); as[i][2]=read(); as[i][3]=read(); as[i][4]=read();
        x=as[i][3]; y=as[i][4]; l=1; r=m; v=1;
        while (l<r){
            mid=(l+r)/2;
            if (x<=mid && y>=mid){
                break;
            }
            if (y<=mid){
                r=mid; v=v*2;
            } else{
                l=mid+1;
                v=v*2+1;
            }
        }
        add(v,i);
    }
    na=0x7fffffff; na=na*10000;
    dfs(1,1,m);
    fo(i,1,q) printf("%lld\n",an[i]);
    return 0;
}