• 4010: [HNOI2015]菜肴制作


    4010: [HNOI2015]菜肴制作

    Time Limit: 5 Sec  Memory Limit: 512 MB
    Submit: 446  Solved: 251
    [Submit][Status][Discuss]

    Description

    知名美食家小 A被邀请至ATM 大酒店,为其品评菜肴。 

    ATM 酒店为小 A 准备了 N 道菜肴,酒店按照为菜肴预估的质量从高到低给予
    1到N的顺序编号,预估质量最高的菜肴编号为1。由于菜肴之间口味搭配的问题,
    某些菜肴必须在另一些菜肴之前制作,具体的,一共有 M 条形如“i 号菜肴‘必须’
    先于 j 号菜肴制作”的限制,我们将这样的限制简写为<i,j>。现在,酒店希望能求
    出一个最优的菜肴的制作顺序,使得小 A能尽量先吃到质量高的菜肴:也就是说,
    (1)在满足所有限制的前提下,1 号菜肴“尽量”优先制作;(2)在满足所有限制,1
    号菜肴“尽量”优先制作的前提下,2号菜肴“尽量”优先制作;(3)在满足所有限
    制,1号和2号菜肴“尽量”优先的前提下,3号菜肴“尽量”优先制作;(4)在满
    足所有限制,1 号和 2 号和 3 号菜肴“尽量”优先的前提下,4 号菜肴“尽量”优
    先制作;(5)以此类推。 
    例1:共4 道菜肴,两条限制<3,1>、<4,1>,那么制作顺序是 3,4,1,2。例2:共
    5道菜肴,两条限制<5,2>、 <4,3>,那么制作顺序是 1,5,2,4,3。例1里,首先考虑 1,
    因为有限制<3,1>和<4,1>,所以只有制作完 3 和 4 后才能制作 1,而根据(3),3 号
    又应“尽量”比 4 号优先,所以当前可确定前三道菜的制作顺序是 3,4,1;接下来
    考虑2,确定最终的制作顺序是 3,4,1,2。例 2里,首先制作 1是不违背限制的;接
    下来考虑 2 时有<5,2>的限制,所以接下来先制作 5 再制作 2;接下来考虑 3 时有
    <4,3>的限制,所以接下来先制作 4再制作 3,从而最终的顺序是 1,5,2,4,3。 
    现在你需要求出这个最优的菜肴制作顺序。无解输出“Impossible!” (不含引号,
    首字母大写,其余字母小写) 

    Input

     第一行是一个正整数D,表示数据组数。 

    接下来是D组数据。 
    对于每组数据: 
    第一行两个用空格分开的正整数N和M,分别表示菜肴数目和制作顺序限
    制的条目数。 
    接下来M行,每行两个正整数x,y,表示“x号菜肴必须先于y号菜肴制作”
    的限制。(注意:M条限制中可能存在完全相同的限制) 

    Output

     输出文件仅包含 D 行,每行 N 个整数,表示最优的菜肴制作顺序,或

    者”Impossible!”表示无解(不含引号)。 

    Sample Input

    3
    5 4
    5 4
    5 3
    4 2
    3 2
    3 3
    1 2
    2 3
    3 1
    5 2
    5 2
    4 3

    Sample Output

    1 5 3 4 2
    Impossible!
    1 5 2 4 3

    HINT

     【样例解释】 


    第二组数据同时要求菜肴1先于菜肴2制作,菜肴2先于菜肴3制作,菜肴3先于

    菜肴1制作,而这是无论如何也不可能满足的,从而导致无解。 

    100%的数据满足N,M<=100000,D<=3。 

     

    Source

    题解:乍一看,做菜必须遵循一定的顺序,那么基本可以断定和TopSort有关,然后再看题目,貌似还要小号尽可能靠前(HansBug:特别提醒注意——是小号尽可能靠前,不是字典序最小,两者有很大区别),所以一开始我的反应就是求字典序最小的的拓扑序列,然后做出来一运行就粗线了前面所说的情况= =

    然而,这个问题也不难解决,既然正着没办法,还不如反着来——1.反向建图 2.反向字典序最大,然后啪啪啪求了出来,然后倒着输出,就没了

    PS:1.本蒟蒻由于太弱了,一时没想起来普通的堆怎么写(HansBug:只会写左偏树有木有TT),于是逗比的种了一棵线段树,功能是单点修改+查询最大位置的0(HansBug:好吧我比较逗请无视QAQ)

    2.特别注意,这题数据貌似在每一行非Impossible!的输出最后貌似要保留一个空格,我一开始没保留结果PE了

      1 /**************************************************************
      2     Problem: 4010
      3     User: HansBug
      4     Language: Pascal
      5     Result: Accepted
      6     Time:2904 ms
      7     Memory:11180 kb
      8 ****************************************************************/
      9  
     10 type
     11     point=^node;
     12     node=record
     13                g:longint;
     14                next:point;
     15     end;
     16 var
     17    a:array[0..200000] of point;
     18    e,b,c,f:array[0..100005] of longint;
     19    d:array[0..1000005] of longint;
     20    i,j,k,l,m,n:longint;
     21 function max(x,y:longint):longint;
     22          begin
     23               if x>y then max:=x else max:=y;
     24          end;
     25 function min(x,y:longint):longint;
     26          begin
     27               if x<y then min:=x else min:=y;
     28          end;
     29 procedure swap(var x,y:longint);
     30           var z:longint;
     31           begin
     32                z:=x;x:=y;y:=z;
     33           end;
     34 procedure add(x,y:longint);
     35           var p:point;
     36           begin
     37                new(p);p^.g:=y;p^.next:=a[x];a[x]:=p;
     38           end;
     39 procedure built(z,x,y:longint);
     40           begin
     41                if x=y then
     42                   begin
     43                        d[z]:=e[x];
     44                        b[x]:=z;
     45                   end
     46                else
     47                    begin
     48                         built(z*2,x,(x+y) div 2);
     49                         built(z*2+1,(x+y) div 2+1,y);
     50                         d[z]:=min(d[z*2],d[z*2+1]);
     51                    end;
     52           end;
     53 function getzero(z,x,y:longint):longint;
     54          begin
     55               if d[z]<>0 then exit(-1);
     56               if (x=y) and (d[z]=0) then exit(x);
     57               if d[z*2+1]=0 then exit(getzero(z*2+1,(x+y) div 2+1,y))
     58                  else exit(getzero(z*2,x,(x+y) div 2));
     59          end;
     60 procedure beep(x,y:longint);
     61           begin
     62                x:=b[x];
     63                if y=maxlongint then d[x]:=maxlongint else inc(d[x],y);
     64                while x>1 do
     65                      begin
     66                           x:=x div 2;
     67                           d[x]:=min(d[x*2],d[x*2+1]);
     68                      end;
     69           end;
     70 procedure pd;
     71           var i,j,k,l,n,m:longint;p:point;
     72           begin
     73                readln(n,m);
     74                fillchar(e,sizeof(e),0);
     75                fillchar(f,sizeof(f),0);
     76                fillchar(c,sizeof(c),0);
     77                for i:=0 to n do a[i]:=nil;
     78                for i:=1 to m do
     79                    begin
     80                         readln(j,k);
     81                         add(k,j);inc(e[j]);
     82                    end;
     83                built(1,1,n);
     84                for i:=1 to n do
     85                    begin
     86                         f[i]:=getzero(1,1,n);
     87                         if f[i]=-1 then
     88                            begin
     89                                 writeln('Impossible!');
     90                                 exit;
     91                            end;
     92                         c[f[i]]:=1;beep(f[i],maxlongint);
     93                         p:=a[f[i]];
     94                         while p<>nil do
     95                               begin
     96                                    if c[p^.g]=0 then beep(p^.g,-1);
     97                                    p:=p^.next;
     98                               end;
     99                    end;
    100                for i:=n downto 1 do write(f[i],' ');
    101                writeln;
    102                i:=0;
    103           end;
    104 begin
    105      readln(n);
    106      for i:=1 to n do pd;
    107      readln;
    108  
    109 end.
     
  • 相关阅读:
    Linux内核的异常修复原理
    sudo: insmod: command not found
    在Qemu+ARM上运行Minix3内核
    2021.34 面对干扰
    2021.33 实践
    selenium+python自动化106
    python测试开发django-111.模型管理器(models.Manager)
    python笔记64
    python笔记63
    python笔记62
  • 原文地址:https://www.cnblogs.com/HansBug/p/4491369.html
Copyright © 2020-2023  润新知