BZOJ1562: [NOI2009]变换序列

1562: [NOI2009]变换序列

Time Limit: 5 Sec  Memory Limit: 64 MB
Submit: 1148  Solved: 599
[Submit][Status]

Description

Input

Output

Sample Input

5
1 1 2 2 1

Sample Output

1 2 4 0 3

HINT

30%的数据中N≤50；
60%的数据中N≤500；
100%的数据中N≤10000。

Source

题解：

不愧是NOI的题目，果然是好题。

具体题解建BYVOID的博客：https://www.byvoid.com/blog/noi-2009-transform/ 

讲的非常清楚，个人认为第二种解法是容易想到的，这其中有很多细节需要注意：

1.在寻找最小字典序的时候，加入一个当前查询点，使得在寻找增广路的时候如果x<=cur 那么返回false

   这是我在思考的时候没有想到的，就是如何能让后面的点在更新解的时候不妨碍到前面已经更新过的点

   这是一个非常好的思路，十分巧妙

2.我在刚开始想的时候想到每个点一定有两个匹配点吗？答案是不一定的

   于是我在程序中进行了进一步的验证，但在后面的解题中我并没有再考虑a[i]=-1 or b[i]=-1的情况，同样得到了正确的答案

   应该是cur的限制卡掉了这些-1

代码：

 

const maxn=10000+100;maxm=20000+100;
type node=record
     go,next:longint;
     end;
var e:array[0..maxm] of node;
    head,a,b,px,py,c:array[0..maxn] of longint;
    v:array[0..maxn] of boolean;
    cur,i,n,tot,ans,tmp:longint;
    procedure swap(var x,y:longint);
     var t:longint;
         begin
         t:=x;x:=y;y:=t;
         end;

procedure insert(x,y:longint);
 begin
 inc(tot);
 e[tot].go:=y;e[tot].next:=head[x];head[x]:=tot;
 end;
function find(x:longint):boolean;
 var i,y:longint;
 begin
 if x<=cur then exit(false);
 i:=head[x];
 while i<>0 do
  begin
  y:=e[i].go;
  if not(v[y]) then
   begin
   v[y]:=true;
   if (py[y]=-1) or (find(py[y])) then
    begin
    py[y]:=x;
    px[x]:=y;
    exit(true);
    end;
   end;
  i:=e[i].next;
  end;
 exit(false);
 end;
procedure init;
 begin
 readln(n);
 for i:=0 to n-1 do
   begin
   read(c[i]);
   a[i]:=(i+c[i]) mod n;
   if n-abs(a[i]-i)<c[i] then a[i]:=-1;
   b[i]:=(i+n-c[i]) mod n;
   if abs(b[i]-i)<c[i] then b[i]:=-1;
   if a[i]>b[i] then swap(a[i],b[i]);
   if a[i]<>-1 then insert(i,a[i]);
   if b[i]<>-1 then insert(i,b[i]);
   writeln(a[i],' ',b[i]);
   end;
 end;
procedure main;
 begin
 for i:=0 to n-1 do py[i]:=-1;
 cur:=-1;
 for i:=0 to n-1 do
  begin
  fillchar(v,sizeof(v),false);
  if find(i) then inc(ans);
  end;
 if ans<>n then begin writeln('No Answer');exit;end;
 fillchar(v,sizeof(v),false);
 for i:=0 to n-1 do
  begin
  if px[i]<>a[i] then
    begin
    cur:=i;
    fillchar(v,sizeof(v),false);
    tmp:=py[a[i]];
    py[a[i]]:=i;
    py[b[i]]:=-1;
    v[a[i]]:=true;
    if find(tmp) then
      begin
       px[i]:=a[i]
      end
    else
      begin
       py[a[i]]:=tmp;
       py[b[i]]:=i;
      end;
    end;
  end;
 write(px[0]);for i:=1 to n-1 do write(' ',px[i]);
 end;
begin
 assign(input,'input.txt');assign(output,'output.txt');
 reset(input);rewrite(output);
 init;
 main;
 close(input);close(output);
end.