[HDOJ 4984] Goffi and Graph [最短路]

给定一个图，每条边的权值都是随时间变化的，为一个一次函数。定义f(i,t)为t时刻，从点1到点i的所有的路中，权值最小的边的最大值。F(i)为f(i)从0到T积分后的结果。问Sum_{i=2..n}{F(i)}。

我们可以求出每条边的一次函数的交点，然后在两个交点中间的位置，边的权重之间的大小关系是不变的。因为最多只有100条边，所以每次用dijstra求出f(i)即可。

据说可以求最大生成树。然后直接在树上dfs也可以求出f(i)。

#include <cstdio>
#include <cstring>
#include <set>
#include <queue>

inline double min(double a,double b) {
    return a<b?a:b;
}

struct Node {
    int fe,frome;
    double v;
};
struct Edge {
    int t,ne,a,b;
    double v;
    void calV(double ti) {
        v=a+b*ti;
    }
    double calAns(double bg,double ed) {
        return b*(ed*ed-bg*bg)/2+a*(ed-bg);
    }
};
struct QueNode {
    int i;
    double v;
    QueNode (int ii,double vv) {
        i=ii;v=vv;
    }
    friend bool operator < (const QueNode &a,const QueNode &b) {
        return a.v<b.v;
    }
};

std::set<double> c;
std::priority_queue<QueNode> d;
Node a[51];
Edge b[200];
int n,m,ti,bp;

void putEdge(int x,int y,int aa,int bb) {
    b[bp].t=y;
    b[bp].a=aa;
    b[bp].b=bb;
    b[bp].ne=a[x].fe;
    a[x].fe=bp++;
}
void addTime(const Edge &x,const Edge &y) {
    if (x.b==y.b) return;
    c.insert((double)(y.a-x.a)/(x.b-y.b));
}
double calAns(double bg,double ed) {
    double ti=(bg+ed)/2;
    int i,j;
    for (j=0;j<bp;j++) b[j].calV(ti);
    for (i=1;i<=n;i++) a[i].v=-1e-10;
    a[1].v=1e10;
    d.push(QueNode(1,a[1].v));
    while (!d.empty()) {
        int i=d.top().i;
        double v=d.top().v;
        d.pop();
        if (v==a[i].v) {
            //printf("--%d %lf\n",i,v);
            for (j=a[i].fe;j!=-1;j=b[j].ne) {
                int ii=b[j].t;
                double vv=min(v,b[j].v);
                if (a[ii].v<vv) {
                    a[ii].v=vv;
                    if (v>b[j].v) a[ii].frome=j;
                    else a[ii].frome=a[i].frome;
                    d.push(QueNode(ii,vv));
                    //printf("%d %lf %lf\n",ii,vv,b[j].v);
                }
            }
        }
    }
    //for (i=1;i<=n;i++) printf("%lf\n",a[i].v);
    double ans=0;
    for (i=2;i<=n;i++) {
        ans+=b[a[i].frome].calAns(bg,ed);
    }
    //printf("%lf %lf %lf\n",bg,ed,ans);
    return ans;
}

int main() {
    int t,i,j;
    scanf("%d",&t);
    while (t--) {
        scanf("%d%d%d",&n,&m,&ti);
        for (i=1;i<=n;i++) a[i].fe=-1;
        c.clear();
        bp=0;
        for (i=0;i<m;i++) {
            int x,y,a,b;
            scanf("%d%d%d%d",&x,&y,&a,&b);
            putEdge(x,y,a,b);
            putEdge(y,x,a,b);
        }
        for (i=0;i<bp;i+=2) {
            for (j=i+2;j<bp;j+=2) {
                addTime(b[i],b[j]);
            }
        }
        double ans=0;
        double lastTime=0;
        std::set<double>::iterator it=c.begin();
        while (*it<0) it++;
        for (;it!=c.end()&&*it<ti;it++) {
            ans+=calAns(lastTime,*it);
            lastTime=*it;
        }
        ans+=calAns(lastTime,ti);
        printf("%.3lf\n",ans/ti);
    }
    return 0;
}