Noi 2009 day1 第一题 变换序列

本文解析了一道来自NOI2009的竞赛题——变换序列,介绍了题目的背景与要求,并详细阐述了如何利用二分图的最大匹配理论(匈牙利算法)来解决该问题。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题目:

 变换序列

来源:

 Noi 2009 day1 第一题

题目大意:

    对于N个整数0,1,...,N-1,一个变换序列T可以将i变成Ti,其中Ti∈{0,1,...,N-1}且 。  U……(图片无法显示,请连接    题目来源)∀;x,y∈{0,1,...,N-1},定义x和y之间的  距离D(x,y)=min{|x-y|,N-|x-y|} 。给定每个i和  Ti之间的距离D(i,Ti),你需要求出一个满  足要求的变换序列T。如果有多个满足条件  的序列,输出其中字典序最小的一个。

数据范围:

 N≤10000

样例:

 5

 1 1 2 2 1

1 2 4 0 3

做题思路:

 刚看到题,从理解题以上就让我多次丧失信心,终于找到突破,询问了同学匈牙利算  法,然后在网上通过stall4(Usaco)习得此算法,发现这个题就是一个二分图求最  大匹配,好了,大头出来了就好做了

知识点:

 二分图求最大匹配(匈牙利算法)

var
 a:array[0..10000,0..10000]of boolean;
 b:array[0..10000]of boolean;
 c:array[0..10000]of longint;
 watch:array[0..10000]of longint;
 n,m,y,i,k,j,ans:longint;
function path(x:longint):boolean;
var
 i:longint;
begin
 for i:=1 to n do
 if (a[x,i])and(not b[i]) then
  begin
    b[i]:=true;
    if (c[i]=0)or(path(c[i])) then
     begin
      c[i]:=x;
      watch[x]:=i;{<记录对应序列的点,不断更新的,不要把它与c[]混了昂,c每次都要初始化的>}
      exit(true);
     end;
  end;
 exit(false);
end;
procedure hungary;
var
 i:longint;
begin
 fillchar(c,sizeof(c),0);
 for i:=n downto 1 do
 begin
  fillchar(b,sizeof(b),false);
  if not path(i) then
    begin
     writeln('No Answer');
     halt;
    end;
 end;
end;
begin
 assign(input,'transform.in');reset(input);
 assign(output,'transform.out');rewrite(output);
 fillchar(a,sizeof(a),false);
 readln(n);
 for i:=1 to n do
 begin
  read(y);
  a[i,(i+y-1)mod n+1]:=true; {<向右距y的>}
  a[i,(i-y+n-1)mod n+1]:=true;{<向左距y的>}
 end;
 ans:=0;
 hungary;
 for i:=1 to n-1 do write(watch[i]-1,' ');
 writeln(watch[n]-1);{<貌似最好一个数字后不要空格,所以我谨慎的这么写了>}
 close(input);close(output);
end.
题目来源: http://www.rqnoj.cn/Problem_522.html

目描述 给出一个由 $n$ 个数排成的序列,你可以进行如下操作: 交换相邻的两个数。 将该序列中的第 $k$ 个数改为 $x$。 你的任务是将该序列变成一个等差数列。 输入格式 第一行包含两个整数 $n$ 和 $m$,表示序列长度和操作次数。 第二行包含 $n$ 个整数,表示初始序列。 接下来 $m$ 行,每行描述一个操作,格式如下: 1. swap x y:表示交换原序列中下标为 $x$ 和 $y$ 的数。 2. set x y:表示将原序列中下标为 $x$ 的数修改为 $y$。 输出格式 如果无法将序列变成等差数列,输出 Impossible。 否则,第一行输出 Yes。 第二行输出变换后的序列。 数据范围 $1 \leq n \leq 5000$ $1 \leq m \leq 10^5$ $1 \leq |a_i|, |x|, |y| \leq 10^9$ 输入样例1: 5 3 5 4 3 2 1 set 2 6 swap 1 2 set 4 1 输出样例1: Yes 6 4 3 1 2 输入样例2: 5 3 5 4 3 2 1 set 2 6 swap 1 2 set 4 10 输出样例2: Impossible 算法1 (二分图匹配) $O(n^3)$ 由于每次可以交换相邻的两个数,可以通过这个操作来使序列中的某个数字移动到另一个位置,因此可以考虑每个数字在序列中的位置之间建立一条边,如果两个数字之间可以通过交换相邻数变成等差数列,则在它们之间连一条边。 具体而言,如果有两个数字 $a,b$,假设它们当前的位置分别为 $i,j$,那么它们之间就可以连一条边,如果它们能够使得序列变成等差数列。 对于 $set$ 操作,可以考虑将该位置对应的点与其他点之间的边全部删除,然后重新建立连接。 最后,判断是否存在一个匹配,使得所有点都与另一个点进行匹配,如果存在,则说明可以通过交换相邻数的方式使得序列变成等差数列。 时间复杂度 - 建图 $O(n^2)$ - 匈牙利算法 $O(n^3)$ C++ 代码
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值