本题是运筹学的经典模型——流水线调度作业模型的 Johnson 不等式的运用。其他的常见模型可参考
运营管理学习笔记-第十章。
Johnson 不等式是 S. M. Johnson 在 1954 年提出的用于解决两台机器流水车间调度问题以最小化最大完工时间的最优算法。它是一种贪心算法,步骤如下:
将所有 N N N 架飞船分成两个集合。集合 S 1 S_1 S 1 包含所有满足 U ( i ) ≤ V ( i ) U(i) \le V(i) U ( i ) ≤ V ( i ) 的飞船 i i i 。集合 S 2 S_2 S 2 包含所有满足 U ( i ) > V ( i ) U(i) > V(i) U ( i ) > V ( i ) 的飞船 i i i 。
对集合 S 1 S_1 S 1 中的飞船,按照它们的 U ( i ) U(i) U ( i ) 值进行非递减 (升序) 排序。对集合 S 2 S_2 S 2 中的飞船,按照它们的 V ( i ) V(i) V ( i ) 值进行非递增 (降序) 排序。
最终的最优处理顺序是先按照排序后的顺序处理 S 1 S_1 S 1 中的所有飞船,然后紧接着按照排序后的顺序处理 S 2 S_2 S 2 中的所有飞船。
这也就是其他题解中所写的:
U U U 值小的尽量往前放,
V V V 值大的尽量往后放。如果
U ≤ V U \le V U ≤ V ,就看
U U U 值,越小越靠前;如果
U > V U > V U > V ,就看
V V V 值,越大越靠后。这样进行排序统计答案,即可通过本题。
这样做确实能完成本题。但是我不认为这个贪心是非常直观显然的。而这一个贪心的证明也并不简单,甚至可以说是十分繁琐。证明 Johnson 不等式可以使用邻项交换法完成。下面给出一个简要证明。
假设存在一个飞船处理序列
π \pi π ,其不完全符合 Johnson 不等式,也即存在一对相邻的飞船
( i , j ) (i,j) ( i , j ) ,它们的顺序不符合 Johnson 不等式。我们需要证明的是,将
( i , j ) (i,j) ( i , j ) 交换可以让总花费时间不增加的同时,其符合 Johnson 不等式的排序思路,成为最优方案。为方便叙述,我们定义
C k , 1 ( π ) C_{k,1}(\pi) C k , 1 ( π ) 表示在一个飞船处理序列
π \pi π 中,第
k k k 艘飞船到达火星的时间,
C k , 2 ( π ) C_{k,2}(\pi) C k , 2 ( π ) 表示第
k k k 艘飞船到达天王星的时间。
假设
i i i 是序列中第
k k k 个元素,那么,原本不交换的飞船处理序列是
π = ( π 1 , π 2 , … , π k − 1 , i , j , π k + 2 , … , π N ) \pi=(\pi_1,\pi_2,\dots,\pi_{k-1},i,j,\pi_{k+2},\dots,\pi_N) π = ( π 1 , π 2 , … , π k − 1 , i , j , π k + 2 , … , π N ) ,而交换后的序列是
π ′ = ( π 1 , π 2 , … , π k − 1 , j , i , π k + 2 , … , π N ) \pi'=(\pi_1,\pi_2,\dots,\pi_{k-1},j,i,\pi_{k+2},\dots,\pi_N) π ′ = ( π 1 , π 2 , … , π k − 1 , j , i , π k + 2 , … , π N ) 。实际上我们只需要比较的是
C k + 1 , 2 ( π ) C_{k+1,2}(\pi) C k + 1 , 2 ( π ) 和
C k + 1 , 2 ( π ′ ) C_{k+1,2}(\pi') C k + 1 , 2 ( π ′ ) 的大小,因为之后的元素不会再产生差异。接下来,为了表述方便,我们简记
A = C k − 1 , 1 ( π ) A=C_{k-1,1}(\pi) A = C k − 1 , 1 ( π ) ,即第
k − 1 k-1 k − 1 艘飞船到达火星的时间,以及
B = C k − 1 , 2 ( π ) B=C_{k-1,2}(\pi) B = C k − 1 , 2 ( π ) ,即第
k − 1 k-1 k − 1 艘飞船到达天王星的时间。根据这些定义和变量,我们可以得出以下的几个数值:
C k , 1 ( π ) = A + U ( i ) C_{k,1}(\pi) = A + U(i) C k , 1 ( π ) = A + U ( i ) ;
C k + 1 , 1 ( π ) = C k , 1 ( π ) + U ( j ) = A + U ( i ) + U ( j ) C_{k+1,1}(\pi) = C_{k,1}(\pi) + U(j) = A + U(i) + U(j) C k + 1 , 1 ( π ) = C k , 1 ( π ) + U ( j ) = A + U ( i ) + U ( j ) ;
C k , 2 ( π ) = max ( C k , 1 ( π ) , B ) + V ( i ) = max ( A + U ( i ) , B ) + V ( i ) C_{k,2}(\pi) = \max(C_{k,1}(\pi), B) + V(i) = \max(A + U(i), B) + V(i) C k , 2 ( π ) = max ( C k , 1 ( π ) , B ) + V ( i ) = max ( A + U ( i ) , B ) + V ( i ) 。这里的 max \max max 函数是因为需要等待前一艘飞船到达天王星后,后一艘飞船才能出发;
C k + 1 , 2 ( π ) = max ( C k + 1 , 1 ( π ) , C k , 2 ( π ) ) + V ( j ) C_{k+1,2}(\pi) = \max(C_{k+1,1}(\pi), C_{k,2}(\pi)) + V(j) C k + 1 , 2 ( π ) = max ( C k + 1 , 1 ( π ) , C k , 2 ( π )) + V ( j ) 。这里的 max \max max 函数指,第 k + 1 k+1 k + 1 艘飞船已经到达了火星待飞天王星,同时前一艘飞船已经到达了天王星后,才可起飞。
代入上述数值可以得出:
C k + 1 , 2 ( π ) = max ( A + U ( i ) + U ( j ) , max ( A + U ( i ) , B ) + V ( i ) ) + V ( j ) C_{k+1,2}(\pi) = \max(A + U(i) + U(j), \max(A + U(i), B) + V(i)) + V(j) C k + 1 , 2 ( π ) = max ( A + U ( i ) + U ( j ) , max ( A + U ( i ) , B ) + V ( i )) + V ( j ) ;
对于飞船处理序列
π ′ \pi' π ′ ,类似地我们有:
C k , 1 ( π ′ ) = A + U ( j ) C_{k,1}(\pi') = A + U(j) C k , 1 ( π ′ ) = A + U ( j ) ;
C k + 1 , 1 ( π ′ ) = C k , 1 ( π ′ ) + U ( i ) = A + U ( j ) + U ( i ) C_{k+1,1}(\pi') = C_{k,1}(\pi') + U(i) = A + U(j) + U(i) C k + 1 , 1 ( π ′ ) = C k , 1 ( π ′ ) + U ( i ) = A + U ( j ) + U ( i ) 。这一数值与 C k + 1 , 1 ( π ) C_{k+1,1}(\pi) C k + 1 , 1 ( π ) 相同;
C k , 2 ( π ′ ) = max ( C k , 1 ( π ′ ) , B ) + V ( j ) = max ( A + U ( j ) , B ) + V ( j ) C_{k,2}(\pi') = \max(C_{k,1}(\pi'), B) + V(j) = \max(A + U(j), B) + V(j) C k , 2 ( π ′ ) = max ( C k , 1 ( π ′ ) , B ) + V ( j ) = max ( A + U ( j ) , B ) + V ( j ) ;
C k + 1 , 2 ( π ′ ) = max ( C k + 1 , 1 ( π ′ ) , C k 2 ( π ′ ) ) + V ( i ) C_{k+1,2}(\pi') = \max(C_{k+1,1}(\pi'), C_{k2}(\pi')) + V(i) C k + 1 , 2 ( π ′ ) = max ( C k + 1 , 1 ( π ′ ) , C k 2 ( π ′ )) + V ( i ) ;
代入上述数值可以得出:
C k + 1 , 2 ( π ′ ) = max ( A + U ( j ) + U ( i ) , max ( A + U ( j ) , B ) + V ( j ) ) + V ( i ) C_{k+1,2}(\pi') = \max(A + U(j) + U(i), \max(A + U(j), B) + V(j)) + V(i) C k + 1 , 2 ( π ′ ) = max ( A + U ( j ) + U ( i ) , max ( A + U ( j ) , B ) + V ( j )) + V ( i )
下面我们证明:如果飞船
i i i 在飞船
j j j 之前,但它们满足
min ( U ( i ) , V ( j ) ) > min ( U ( j ) , V ( i ) ) \min(U(i), V(j)) > \min(U(j), V(i)) min ( U ( i ) , V ( j )) > min ( U ( j ) , V ( i )) ,那么将它们的顺序交换为
( j , i ) (j, i) ( j , i ) ,至少不会增加第
k + 1 k+1 k + 1 个位置飞船的完工时间。
根据
max \max max 函数的运算性质:
max ( x , y ) + z = max ( x + z , y + z ) \max(x,y)+z=\max(x+z,y+z) max ( x , y ) + z = max ( x + z , y + z ) ,展开表达式可得:
C k + 1 , 2 ( π ) = max ( A + U ( i ) + U ( j ) + V ( j ) , A + U ( i ) + V ( i ) + V ( j ) , B + V ( i ) + V ( j ) ) C_{k+1, 2}(\pi) = \max(A+U(i)+U(j)+V(j), A+U(i)+V(i)+V(j), B+V(i)+V(j)) C k + 1 , 2 ( π ) = max ( A + U ( i ) + U ( j ) + V ( j ) , A + U ( i ) + V ( i ) + V ( j ) , B + V ( i ) + V ( j ))
C k + 1 , 2 ( π ′ ) = max ( A + U ( j ) + U ( i ) + V ( i ) , A + U ( j ) + V ( j ) + V ( i ) , B + V ( j ) + V ( i ) ) C_{k+1, 2}(\pi') = \max(A+U(j)+U(i)+V(i), A+U(j)+V(j)+V(i), B+V(j)+V(i)) C k + 1 , 2 ( π ′ ) = max ( A + U ( j ) + U ( i ) + V ( i ) , A + U ( j ) + V ( j ) + V ( i ) , B + V ( j ) + V ( i ))
由于
C k + 1 , 2 ( π ) C_{k+1, 2}(\pi) C k + 1 , 2 ( π ) 和
C k + 1 , 2 ( π ′ ) C_{k+1, 2}(\pi') C k + 1 , 2 ( π ′ ) 的
max \max max 函数中都含有
B + V ( i ) + V ( j ) B+V(i)+V(j) B + V ( i ) + V ( j ) ,我们可以忽略它,而只在意
max \max max 函数的前两项。具体地,我们设
T 1 = A + U ( i ) + U ( j ) + V ( j ) T_1=A+U(i)+U(j)+V(j) T 1 = A + U ( i ) + U ( j ) + V ( j ) ,
T 2 = A + U ( i ) + V ( i ) + V ( j ) T_2=A+U(i)+V(i)+V(j) T 2 = A + U ( i ) + V ( i ) + V ( j ) 。对应地,
T 1 ′ = A + U ( j ) + U ( i ) + V ( i ) T_1'=A+U(j)+U(i)+V(i) T 1 ′ = A + U ( j ) + U ( i ) + V ( i ) ,
T 2 ′ = A + U ( j ) + V ( j ) + V ( i ) T_2'=A+U(j)+V(j)+V(i) T 2 ′ = A + U ( j ) + V ( j ) + V ( i ) 。此时我们有:
T 1 − T 1 ′ = [ A + U ( i ) + U ( j ) + V ( j ) ] − [ A + U ( j ) + U ( i ) + V ( i ) ] = V ( j ) − V ( i ) T_1-T_1'=[A+U(i)+U(j)+V(j)]-[A+U(j)+U(i)+V(i)]=V(j)-V(i) T 1 − T 1 ′ = [ A + U ( i ) + U ( j ) + V ( j )] − [ A + U ( j ) + U ( i ) + V ( i )] = V ( j ) − V ( i ) ;
T 2 − T 2 ′ = [ A + U ( i ) + V ( i ) + V ( j ) ] − [ A + U ( j ) + V ( j ) + V ( i ) ] = U ( i ) − U ( j ) T_2-T_2'=[A+U(i)+V(i)+V(j)]-[A+U(j)+V(j)+V(i)]=U(i)-U(j) T 2 − T 2 ′ = [ A + U ( i ) + V ( i ) + V ( j )] − [ A + U ( j ) + V ( j ) + V ( i )] = U ( i ) − U ( j )
由于
min ( U ( i ) , V ( j ) ) > min ( U ( j ) , V ( i ) ) \min(U(i), V(j)) > \min(U(j), V(i)) min ( U ( i ) , V ( j )) > min ( U ( j ) , V ( i )) ,我们可以进行分类讨论,讨论获取
U ( i ) , U ( j ) , V ( i ) , V ( j ) U(i),U(j),V(i),V(j) U ( i ) , U ( j ) , V ( i ) , V ( j ) 的大小关系:
U ( i ) ≤ V ( j ) U(i)\leq V(j) U ( i ) ≤ V ( j ) :
情况 1:U ( j ) ≤ V ( i ) U(j)\leq V(i) U ( j ) ≤ V ( i ) 。此时 U ( i ) > U ( j ) U(i)>U(j) U ( i ) > U ( j ) 成立,即 T 2 > T 2 ′ T_2>T_2' T 2 > T 2 ′ 。而至于 T 1 − T 1 ′ = V ( j ) − V ( i ) T_1-T_1'=V(j)-V(i) T 1 − T 1 ′ = V ( j ) − V ( i ) ,若 V ( j ) ≥ V ( i ) V(j)\geq V(i) V ( j ) ≥ V ( i ) 则直接得证;若 V ( j ) < V ( i ) V(j)<V(i) V ( j ) < V ( i ) ,我们希望此时的 T 1 ′ T_1' T 1 ′ 不会成为交换后较大的候选值,因此应当比较 T 2 T_2 T 2 和 T 1 ′ T_1' T 1 ′ 。由于 T 2 − T 1 ′ = [ A + U ( i ) + V ( i ) + V ( j ) ] − [ A + U ( j ) + U ( i ) + V ( i ) ] = V ( j ) − U ( j ) T_2-T_1'=[A+U(i)+V(i)+V(j)]-[A+U(j)+U(i)+V(i)]=V(j)-U(j) T 2 − T 1 ′ = [ A + U ( i ) + V ( i ) + V ( j )] − [ A + U ( j ) + U ( i ) + V ( i )] = V ( j ) − U ( j ) ,由于 U ( i ) ≤ V ( j ) U(i)\leq V(j) U ( i ) ≤ V ( j ) ,U ( i ) > U ( j ) U(i)>U(j) U ( i ) > U ( j ) 可以得出 T 2 − T 1 ′ > 0 T_2-T_1'>0 T 2 − T 1 ′ > 0 ,即 T 2 > T 1 ′ T_2>T_1' T 2 > T 1 ′ ,且又有 T 2 > T 2 ′ T_2>T_2' T 2 > T 2 ′ ,从而 max ( T 1 ′ , T 2 ′ ) ≤ max ( T 1 , T 2 ) \max(T_1',T_2')\leq \max(T_1,T_2) max ( T 1 ′ , T 2 ′ ) ≤ max ( T 1 , T 2 ) 。
情况 2:V ( i ) < U ( j ) V(i)<U(j) V ( i ) < U ( j ) 。此时 U ( i ) > V ( i ) U(i)>V(i) U ( i ) > V ( i ) 成立。由于 U ( i ) ≤ V ( j ) U(i)\leq V(j) U ( i ) ≤ V ( j ) ,因此有 V ( j ) > V ( i ) V(j)>V(i) V ( j ) > V ( i ) ,可得 T 1 > T 1 ′ T_1>T_1' T 1 > T 1 ′ 。至于 T 2 − T 2 ′ = U ( i ) − U ( j ) T_2-T_2'=U(i)-U(j) T 2 − T 2 ′ = U ( i ) − U ( j ) ,若 U ( i ) ≥ U ( j ) U(i)\geq U(j) U ( i ) ≥ U ( j ) 则直接得证;若 U ( i ) < U ( j ) U(i)<U(j) U ( i ) < U ( j ) 则考虑 T 1 − T 2 ′ = [ A + U ( i ) + U ( j ) + V ( j ) ] − [ A + U ( j ) + V ( j ) + V ( i ) ] = U ( i ) − V ( i ) T_1-T_2'=[A+U(i)+U(j)+V(j)]-[A+U(j)+V(j)+V(i)]=U(i)-V(i) T 1 − T 2 ′ = [ A + U ( i ) + U ( j ) + V ( j )] − [ A + U ( j ) + V ( j ) + V ( i )] = U ( i ) − V ( i ) 。此时因为 U ( i ) > V ( i ) U(i)>V(i) U ( i ) > V ( i ) ,因此 T 1 > T 2 ′ T_1>T_2' T 1 > T 2 ′ 。又有 T 2 > T 2 ′ T_2>T_2' T 2 > T 2 ′ ,从而 max ( T 1 ′ , T 2 ′ ) ≤ max ( T 1 , T 2 ) \max(T_1',T_2')\leq \max(T_1,T_2) max ( T 1 ′ , T 2 ′ ) ≤ max ( T 1 , T 2 ) 。
V ( j ) < U ( i ) V(j)<U(i) V ( j ) < U ( i ) :
情况 1:U ( j ) ≤ V ( i ) U(j)\leq V(i) U ( j ) ≤ V ( i ) 。证明与上方一致,从略。
情况 2:V ( i ) < U ( j ) V(i)<U(j) V ( i ) < U ( j ) 。证明与上方一致,从略。
由于两个序列中第三个项完全相同,所以可得:
C k + 1 , 2 ( π ′ ) = max { T 1 ′ , T 2 ′ , B + V ( i ) + V ( j ) } ≤ max { T 1 , T 2 , B + V ( i ) + V ( j ) } = C k + 1 , 2 ( π ) C_{k+1,2}(\pi') = \max\{T'_1,\,T'_2,\,B+V(i)+V(j)\} \le \max\{T_1,\,T_2,\,B+V(i)+V(j)\} = C_{k+1,2}(\pi) C k + 1 , 2 ( π ′ ) = max { T 1 ′ , T 2 ′ , B + V ( i ) + V ( j )} ≤ max { T 1 , T 2 , B + V ( i ) + V ( j )} = C k + 1 , 2 ( π )
也即排序方式成立。因此,通过逐步交换所有违反 Johnson 不等式的相邻对
( i , j ) (i,j) ( i , j ) ,调整得到最优的飞船处理序列
π \pi π ,计算该序列所需的时间,即为答案。
接下来从编写代码的角度出发,如果直接根据
min ( U ( i ) , V ( j ) ) ≤ min ( U ( j ) , V ( i ) ) \min(U(i), V(j)) \leq \min(U(j), V(i)) min ( U ( i ) , V ( j )) ≤ min ( U ( j ) , V ( i )) 编写排序的方式是错误的。在这一篇
题解 里面说的非常清楚明白。总而言之,使用本文开头所讲述的方式:
将所有 N N N 架飞船分成两个集合。集合 S 1 S_1 S 1 包含所有满足 U ( i ) ≤ V ( i ) U(i) \le V(i) U ( i ) ≤ V ( i ) 的飞船 i i i 。集合 S 2 S_2 S 2 包含所有满足 U ( i ) > V ( i ) U(i) > V(i) U ( i ) > V ( i ) 的飞船 i i i 。
对集合 S 1 S_1 S 1 中的飞船,按照它们的 U ( i ) U(i) U ( i ) 值进行非递减 (升序) 排序。对集合 S 2 S_2 S 2 中的飞船,按照它们的 V ( i ) V(i) V ( i ) 值进行非递增 (降序) 排序。
即可完成对飞船的排序,从而计算得出结论。从证明贪心方式到错误的排序,都体现出本题难度很高,并非部分用户所认为的只有黄题难度。
参考代码:
CPP #include <bits/stdc++.h>
using namespace std;
int main () {
ios_base::sync_with_stdio (false );
cin.tie (nullptr );
int n;
cin >> n;
vector <pair<int , int >> s (n);
for (int i = 0 ; i<n; i++)
cin >> s[i].first >> s[i].second;
vector <pair<int , int >> g1, g2;
for (int i = 0 ; i<n; i++) {
if (s[i].first <= s[i].second)
g1. push_back (s[i]);
else
g2. push_back (s[i]);
}
sort (g1. begin (), g1. end (), [](auto a, auto b) {
return a.first<b.first;
});
sort (g2. begin (), g2. end (), [](auto a, auto b) {
return a.second>b.second;
});
vector <pair<int , int >> ans;
for (auto s: g1)
ans.push_back (s);
for (auto s: g2)
ans.push_back (s);
long long t1 = 0 , t2 = 0 ;
for (auto s: ans) {
t1 += s.first;
t2 = max (t1, t2) + s.second;
}
cout << t2 << endl;
return 0 ;
}