zoukankan      html  css  js  c++  java
  • 洛谷 P2770 航空路线问题【最大费用最大流】

    记得cnt=1!!因为是无向图所以可以把回来的路看成另一条向东的路。字符串用map处理即可。拆点限制流量,除了1和n是(i,i+n,2)表示可以经过两次,其他点都拆成(i,i+n,1),费用设为1,原图中的边(i,j)连接(i+n,j,1),注意特判掉iq&&jnst直接相连的情况,流量要设为2,然后跑最大费用最大流,如果流量小于2就是无解,否则分别从s和tdfs输出方案。
    记得输出方案前先输出一次起点。

    #include<iostream>
    #include<cstdio>
    #include<cstring>
    #include<queue>
    #include<map>
    #include<string>
    using namespace std;
    const int N=1000005,inf=1e9;
    int n,m,s,t,h[N],cnt=1,dis[N],fr[N],a[N][2],ans,sum;
    bool v[N];
    string c1,c2,ha[N];
    map<string,int>mp;
    struct qwe
    {
    	int ne,no,to,va,c;
    }e[N<<2];
    void add(int u,int v,int w,int c)
    {
    	cnt++;
     	e[cnt].ne=h[u];
     	e[cnt].no=u;
     	e[cnt].to=v;
     	e[cnt].va=w;
     	e[cnt].c=c;
     	h[u]=cnt;
    }
    void ins(int u,int v,int w,int c)
    {//cout<<u<<" "<<v<<" "<<w<<endl;
     	add(u,v,w,c);
     	add(v,u,0,-c);
    }
    bool spfa()
    {
    	for(int i=s;i<=t;i++)
    		dis[i]=-inf;
    	queue<int>q;
    	q.push(s);
    	dis[s]=0;
    	v[s]=1;
    	while(!q.empty())
    	{
    		int u=q.front();
    		q.pop();
    		v[u]=0;
    		for(int i=h[u];i;i=e[i].ne)
    			if(e[i].va&&dis[e[i].to]<dis[u]+e[i].c)
    			{
    				fr[e[i].to]=i;
    				dis[e[i].to]=dis[u]+e[i].c;
    				if(!v[e[i].to])
    				{
    					v[e[i].to]=1;
    					q.push(e[i].to);
    				}
    			}
    	}//cout<<dis[t]<<endl;
    	return dis[t]!=-inf;
    }
    void mcf()
    {
     	int x=inf;
     	for(int i=fr[t];i;i=fr[e[i].no])
     	x=min(x,e[i].va);
    	sum+=x;
    	for(int i=fr[t];i;i=fr[e[i].no])
    	{
    		e[i].va-=x;
    		e[i^1].va+=x;
    		ans+=x*e[i].c;//cout<<e[i].c<<endl;
    	}
    }
    int main()
    {
    	ios::sync_with_stdio(false);
    	cin>>n>>m;
    	s=1,t=n+n;
    	for(int i=1;i<=n;i++)
    	{
    		cin>>c1;
    		mp[c1]=i;
    		ha[i]=c1;
    	}
    	for(int i=1;i<=m;i++)
    	{
    		cin>>c1>>c2;//cout<<c1<<endl<<c2<<endl;
    		int x=mp[c1],y=mp[c2];//cout<<x<<" "<<y<<endl;
    		if(x>y)
    			swap(x,y);
    		if(x==1&&y==n)
    			ins(x+n,y,2,0);
    		else
    			ins(x+n,y,1,0);
    	}
    	ins(1,1+n,2,1);
    	ins(n,n+n,2,1);
    	for(int i=2;i<n;i++)
    		ins(i,i+n,1,1);
    	while(spfa())
    		mcf();//cout<<sum<<" "<<ans<<endl;
    	if(sum<2)
    	{
    		puts("No Solution!");
    		return 0;
    	}
    	cout<<ans-2<<endl;
    	cout<<ha[1]<<endl;
    	for(int i=h[s+n],j,k;i;i=e[i].ne)
            if(!e[i].va&&!(i&1))
    		{
                k=e[i].to;
                while(k)
    			{
                    cout<<ha[k]<<endl;
                    v[k]=1;
                    for(j=h[k+n],k=0;j;j=e[j].ne)
                        if(!e[j].va&&!(j&1))
    					{
                            k=e[j].to;
    						break;
                        }
                }
                break;
            }
        for(int i=h[t-n],j,k;i;i=e[i].ne)
            if(!e[i^1].va&&(i&1)&&!v[e[i].to-n])
    		{
                k=e[i].to-n;
                while(k)
    			{
                    cout<<ha[k]<<endl;
    				v[k]=1;
                    for(j=h[k],k=0;j;j=e[j].ne)
                        if(!e[j^1].va&&(j&1))
    					{
                            k=e[j].to-n;
    						break;
                        }
                }
                break;
            }
    	return 0;
    }
    
  • 相关阅读:
    C#动态显示时间
    死锁问题
    TCP_NODELAY算法使用事项
    二叉搜索树的后序遍历
    从上到下打印二叉树
    栈的压入、弹出序列
    包含min函数的栈
    顺时针打印矩阵
    树的子结构
    合并两个排序链表
  • 原文地址:https://www.cnblogs.com/lokiii/p/8428547.html
Copyright © 2011-2022 走看看