bzoj 1880: [Sdoi2009]Elaxia的路线【spfa+拓扑排序】

有趣啊

先spfa分别求出以s1，t1，s2，t2为起点的最短路，然后把在s1-->t1或者s2-->t2最短路上的边重新建有向图，跑拓扑最长路即可

#include<iostream>

#include<cstdio>

#include<cstring>

#include<queue>

#include<vector>

using namespace std;

const int N=1505,inf=1e9;

int n,m,x1,x2,y1,y2,d[N],dis[N],h[N],cnt,s1[N],s2[N],t1[N],t2[N];

bool v[N];

struct qwe

{

    int ne,no,to,va;

}e[1000005];

struct bian

{

	int u,v,w;

	bian(int U=0,int V=0,int W=0)

	{

		u=U,v=V,w=W;

	}

};

vector<bian>ve;

void add(int u,int v,int w)

{

	cnt++;

    e[cnt].ne=h[u];

    e[cnt].no=u;

    e[cnt].to=v;

    e[cnt].va=w;

    h[u]=cnt;

}

void spfa(int s,int dis[])

{

	queue<int>q;

    for(int i=1;i<=n;i++)

		dis[i]=inf,v[i]=0;

    dis[s]=0;

	v[s]=1;

    q.push(s);

    while(!q.empty())

    {

        int u=q.front();

        q.pop();

        v[u]=0;

        for(int i=h[u];i;i=e[i].ne)

            if(dis[e[i].to]>dis[u]+e[i].va)

            {

                dis[e[i].to]=dis[u]+e[i].va;

                if(!v[e[i].to])

				{

					v[e[i].to]=1;

                    q.push(e[i].to);

				}

            }

    }

}

int main()

{

    scanf("%d%d%d%d%d%d",&n,&m,&x1,&y1,&x2,&y2);

    for(int i=1,x,y,v;i<=m;i++)

    {

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

        add(x,y,v),add(y,x,v);

    }

    spfa(x1,s1);

    spfa(y1,t1);

    spfa(x2,s2);

    spfa(y2,t2);

    for(int i=1;i<=cnt;i+=2)

        if(min(s1[e[i].no],s1[e[i].to])+min(t1[e[i].to],t1[e[i].no])+e[i].va==s1[y1]&&min(s2[e[i].no],s2[e[i].to])+min(t2[e[i].no],t2[e[i].to])+e[i].va==s2[y2])

        {

            if(s1[e[i].no]<s1[e[i].to])

				ve.push_back(bian(e[i].no,e[i].to,e[i].va));

            else

				ve.push_back(bian(e[i].to,e[i].no,e[i].va));

        }

	memset(h,0,sizeof(h));

	cnt=0;

	for(int i=0;i<ve.size();i++)

		add(ve[i].u,ve[i].v,ve[i].w),d[ve[i].v]++;

	queue<int>q;

    for(int i=1;i<=n;i++)

    {

        dis[i]=0;

        if(!d[i])

			q.push(i);

    }

    int ans=0;

    while(!q.empty())

    {

        int u=q.front();

        q.pop();

        for(int i=h[u];i;i=e[i].ne)

        {

            dis[e[i].to]=max(dis[e[i].to],dis[u]+e[i].va);

            ans=max(ans,dis[e[i].to]);

            d[e[i].to]--;

            if(!d[e[i].to])

				q.push(e[i].to);

        }

    }

    printf("%d\n",ans);

    return 0;

}

巴特西

bzoj 1880: [Sdoi2009]Elaxia的路线【spfa+拓扑排序】

最新文章

热门文章