zoukankan      html  css  js  c++  java
  • POJ 1094差分约束系统拓扑排序

    Sorting It All Out
    Time Limit: 1000MS   Memory Limit: 10000K
    Total Submissions: 32933   Accepted: 11445

    Description

    An ascending sorted sequence of distinct values is one in which some form of a less-than operator is used to order the elements from smallest to largest. For example, the sorted sequence A, B, C, D implies that A < B, B < C and C < D. in this problem, we will give you a set of relations of the form A < B and ask you to determine whether a sorted order has been specified or not.

    Input

    Input consists of multiple problem instances. Each instance starts with a line containing two positive integers n and m. the first value indicated the number of objects to sort, where 2 <= n <= 26. The objects to be sorted will be the first n characters of the uppercase alphabet. The second value m indicates the number of relations of the form A < B which will be given in this problem instance. Next will be m lines, each containing one such relation consisting of three characters: an uppercase letter, the character "<" and a second uppercase letter. No letter will be outside the range of the first n letters of the alphabet. Values of n = m = 0 indicate end of input.

    Output

    For each problem instance, output consists of one line. This line should be one of the following three:

    Sorted sequence determined after xxx relations: yyy...y.
    Sorted sequence cannot be determined.
    Inconsistency found after xxx relations.

    where xxx is the number of relations processed at the time either a sorted sequence is determined or an inconsistency is found, whichever comes first, and yyy...y is the sorted, ascending sequence.

    Sample Input

    4 6
    A<B
    A<C
    B<C
    C<D
    B<D
    A<B
    3 2
    A<B
    B<A
    26 1
    A<Z
    0 0
    

    Sample Output

    Sorted sequence determined after 4 relations: ABCD.
    Inconsistency found after 2 relations.
    Sorted sequence cannot be determined.
    

    题意:

    题意:给定一组字母的大小关系判断它们是否能组成唯一的拓扑序列。

    算法:拓扑排序:先判断是否有环,再判断是否有序,最后才能判断是否能得出结果。必须遍历完整个图才能知道是否有环!!!


    第一次做拓扑排序类的题,需要知道几个概念:


    拓扑排序:

    一.定义

        对一个有向无环图(Directed Acyclic Graph简称DAG)G进行拓扑排序,是将G中所有顶点排成一个线性序列,使得图中任意一对顶点u和v,若<u,v> ∈E(G),则u在线性序列中出现在v之前。

       通常,这样的线性序列称为满足拓扑次序(Topological Order)的序列,简称拓扑序列。

       注意:

       1)只有有向无环图才存在拓扑序列;

       2)对于一个DAG,可能存在多个拓扑序列;

       如:

       该DAG的拓扑序列为A B C D或者A C B D

     而此有向图是不存在拓扑序列的,因为图中存在环路

    二.拓扑序列算法思想

     (1)从有向图中选取一个没有前驱(即入度为0)的顶点,并输出之;

     

     (2)从有向图中删去此顶点以及所有以它为尾的弧;

         重复上述两步,直至图空,或者图不空但找不到无前驱的顶点为止。
    三.代码实现
        采用邻接矩阵实现,map[i][j]=0,表示节点i和j没有关联;map[i][j]=1,表示存在边<i,j>,并且j的入度加1;

    C语言代码实现:(无队列输出)

    #include<iostream>
    #include<stdlib.h>
    #include<stdio.h>
    #define MAX 100
    usingnamespace std;
    
    void toposort(int map[MAX][MAX],int indegree[MAX],int n)
    {
        int i,j,k;
        for(i=0;i<n;i++) //遍历n次
        {
            for(j=0;j<n;j++) //找出入度为0的节点
            {
                if(indegree[j]==0)
                {
                    indegree[j]--;
                    cout<<j<<endl;
                    for(k=0;k<n;k++) //删除与该节点关联的边
                    {
                        if(map[j][k]==1)
                        {
                            indegree[k]--;
                        }
                    }
                    break;
                }
            }
        }
    }
    
    
    int main(void)
    {
        int n,m; //n:关联的边数,m:节点数
        while(scanf("%d %d",&n,&m)==2&&n!=0)
        {
            int i;
            int x,y;
            int map[MAX][MAX]; //邻接矩阵
            int indegree[MAX]; //入度
            memset(map,0,sizeof(map));
            memset(indegree,0,sizeof(indegree));
            for(i=0;i<n;i++)
            {
                scanf("%d %d",&x,&y);
                if(!map[x][y])
                {
                    map[x][y]=1;
                    indegree[y]++;
                }
            }
            toposort(map,indegree,m);
        }
        return0;
    }

    差分约束系统:


    就是给出一些形如x-y<=b不等式的约束,问你是否满足有解的问题,可以转换成图论里的最短路径问题,下面开始详细介绍下

    比如给出三个不等式,b-a<=k1,c-b<=k2,c-a<=k3,求出c-a的最大值,我们可以把a,b,c转换成三个点,k1,k2,k3是边上的权,如图

    由题我们可以得知,这个有向图中,由题b-a<=k1,c-b<=k2,得出c-a<=k1+k2,因此比较k1+k2和k3的大小,求出最小的就是c-a的最大值了

    根据以上的解法,我们可能会猜到求解过程实际就是求从a到c的最短路径,没错的....简单的说就是从a到c沿着某条路径后把所有权值和k求出就是c -a<=k的一个

    推广的不等式约束,既然这样,满足题目的肯定是最小的k,也就是从a到c最短距离...

    理解了这里之后,想做题还是比较有困难的,因为题目需要变形一下,不能单纯的算..

    首先以poj3159为例,这个比较简单,就是给出两个点的最大差,然后让你求1到n的最大差,直接建图后用bellman或者spfa就可以过了

    稍微难点的就是poj1364,因为他给出的不等式不是x-y<=k形式,有时候是大于号,这样需要我们去变形一下,并且给出的还是>,<没有等于,都要变形

    再有就是poj1201,他要求出的是最长距离,那就要把形式变换成x-y>=k的标准形式

    注意点:

    1. 如果要求最大值想办法把每个不等式变为标准x-y<=k的形式,然后建立一条从y到x权值为k的边,变得时候注意x-y<k =>x-y<=k-1

       如果要求最小值的话,变为x-y>=k的标准形式,然后建立一条从y到x的k边,求出最长路径即可

    2.如果权值为正,用dj,spfa,bellman都可以,如果为负不能用dj,并且需要判断是否有负环,有的话就不存在


    拓扑排序的核心思想:

    拓扑排序的核心就是每次找入度为0的点 进入输出队列 然后将与此点相连的节点入度减1 重复做当做n-1 次后还有点没进输出队列 那么这些点就是环上的 因为环上的各点入度都为1 没有0的 就不能更新。

    AC:

    #include <iostream>
    #include <cstring>
    #include <cstdio>
    using namespace std;
    
    int indegree[27];     //入度
    int map[27][27];     //图
    int queue[27];       // 每次找到入度为0的点后,存储要输出的顺序
    
    int TopoSort(int n)     //拓扑排序
    {
        int count=0;     //记录解空间中零入度顶点的个数
        int temp[27];     //对入度顶点备份
        int pos;      //记录一个零入度顶点的数目
        int i,j;
        int m;       //零入度顶点的个数
        int flag=1;      //flag=1:有序  flag=-1:不能确定
        for (i=1; i<=n; i++)
        {
            temp[i] = indegree[i];  //备份
        }
        for (i=1; i<=n; i++)  //遍历n遍,必须把图全部遍历完(根据当前输入的,尽管map全部初始化为0了,只要根据当前输入的所有条件判断)
        {
            m = 0;
            for (j=1; j<=n; j++)  //查找零入度顶点的个数
            {
                if (temp[j] == 0)
                {
                    m ++;
                    pos = j;   //记录一个零入度顶点的位置
                }
            }
            if (m == 0)     //零入度顶点的个数==0:有环
            {
                return 0;
            }
      if (m > 1)     //零入度顶点的个数>1,说明还有其他两点之间没确定关系,判断是无序(只有1到n个点确定关系,且没有环,根据拓扑排序才能确定有序)
      {
                flag=-1;    //当知道无序时,并不一定知道该图是否有环,因此要继续遍历,可别急着退出...
            }
            queue[count++] = pos;   //零入度顶点入队
            temp[pos] = -1;    //将零入度顶点的入度置为-1
            for (j=1; j<=n; j++)  //删除以pos为起点的边
            {
                if (map[pos][j] == 1)
                {
                    temp[j] --;   //相应定点的入度减1
                }
            }
        }
        return flag;
    }
    int main()
    {
        int n,m;
        int sign;      //当sign=1时,程序已经得出结果,不需再考虑后面的输入
        string str;
        while (scanf("%d%d",&n,&m) && n!=0 && m!=0)
        {
            memset(map,0,sizeof(map));
            memset(indegree,0,sizeof(indegree));
            sign=0;
            for (int i=1; i<=m; i++)
            {
                cin>>str;
                if (sign)
                {
                    continue;   //一旦得出结果,对后续的输入置之不理!
                }
                int u = str[0]-'A'+1;
                int v = str[2]-'A'+1;
                map[u][v] = 1;
                indegree[v] ++;
                int s = TopoSort(n);
                if (s == 0)    //有环
                {
                    printf("Inconsistency found after %d relations.
    ",i);
                    sign=1;
                }
                if (s == 1)    //有序
                {
                    printf("Sorted sequence determined after %d relations: ",i);
                    for (int j=0; j<n; j++)
                    {
                        printf("%c",queue[j]+'A'-1); //根据进队顺序输出字符,整型转字符型
                    }
                    printf(".
    ");//每句话的结尾的点号,记得输出!
                    sign=1;
                }
            }
            if(!sign)     //无法得出结果
            {
                printf("Sorted sequence cannot be determined.
    ");
            }
        }
    
        return 0;
    }
    


  • 相关阅读:
    产品经理的十宗罪,你犯了几宗?
    产品经理的10大顾虑
    【FastAPI 学习 七】GET和POST请求参数接收以及验证
    【FastAPI 学习 六】异常处理
    【FastAPI 学习 五】统一响应json数据格式
    前端展示(三)
    前端展示(二)
    前端设计(一)
    后端流程分析
    生成词云图
  • 原文地址:https://www.cnblogs.com/zhangmingzhao/p/7256387.html
Copyright © 2011-2022 走看看