【训练题】航线设计 | 使用最长上升子序列(LIS)长度的O(nlogn)算法优化
【問題描述】
有一個國家被一條河劃分為南北兩部分,在南岸和北岸總共有N對城鎮,每一城鎮在對岸都有唯一的友好城鎮。任何兩個城鎮都沒有相同的友好城鎮。每一對友好城鎮都希望有一條航線來往。于是他們向政府提出了申請。由于河終年有霧。政府決定不允許有任兩條航線交叉(如果兩條航線交叉,將有很大機會撞船)。
你的任務是寫一個程序來幫政府官員決定他們應撥款興建哪些航線以使得沒有出現交叉的航線最多。
【輸入格式】
第一行一個整數N,表示分布在河兩岸的城鎮對數。接下來的N行每行有兩個由空格分隔的正數C,D (C、D<=10^9),描述每一對友好城鎮沿著河岸與西邊境線的距離,C表示北岸城鎮的距離而D表示南岸城鎮的距離。在河的同一邊,任何兩個城鎮的位置都是不同的。
【輸出格式】
在安全條件下能夠開通的最大航線數目。
【輸入樣例】
7 22 4 2 6 10 3 15 12 9 8 17 17 4 2【輸出樣例】
4【數據范圍】
1<=N<=500000
題解:
O(n2)算法代碼:
優化:
O(nlogn)算法:
最長上升子序列nlogn算法 - mengxm
假設存在一個序列d[1..9] = 2 1 5 3 6 4 8 9 7,可以看出來它的LIS長度為5。
下面一步一步試著找出它。
我們定義一個序列B,然后令 i = 1 to 9 逐個考察這個序列。
此外,我們用一個變量Len來記錄現在最長算到多少了
首先,把d[1]有序地放到B里,令B[1] = 2,就是說當只有1一個數字2的時候,長度為1的LIS的最小末尾是2。這時Len=1
然后,把d[2]有序地放到B里,令B[1] = 1,就是說長度為1的LIS的最小末尾是1,d[1]=2已經沒用了,很容易理解吧。這時Len=1
接著,d[3] = 5,d[3]>B[1],所以令B[1+1]=B[2]=d[3]=5,就是說長度為2的LIS的最小末尾是5,很容易理解吧。這時候B[1..2] = 1, 5,Len=2
再來,d[4] = 3,它正好加在1,5之間,放在1的位置顯然不合適,因為1小于3,長度為1的LIS最小末尾應該是1,這樣很容易推知,長度為2的LIS最小末尾是3,于是可以把5淘汰掉,這時候B[1..2] = 1, 3,Len = 2
繼續,d[5] = 6,它在3后面,因為B[2] = 3, 而6在3后面,于是很容易可以推知B[3] = 6, 這時B[1..3] = 1, 3, 6,還是很容易理解吧? Len = 3 了噢。
第6個, d[6] = 4,你看它在3和6之間,于是我們就可以把6替換掉,得到B[3] = 4。B[1..3] = 1, 3, 4, Len繼續等于3
第7個, d[7] = 8,它很大,比4大,嗯。于是B[4] = 8。Len變成4了
第8個, d[8] = 9,得到B[5] = 9,嗯。Len繼續增大,到5了。
最后一個, d[9] = 7,它在B[3] = 4和B[4] = 8之間,所以我們知道,最新的B[4] =7,B[1..5] = 1, 3, 4, 7, 9,Len = 5。
于是我們知道了LIS的長度為5。
!!!!! 注意。這個1,3,4,7,9不是LIS,它只是存儲的對應長度LIS的最小末尾。有了這個末尾,我們就可以一個一個地插入數據。雖然最后一個d[9] = 7更新進去對于這組數據沒有什么意義,但是如果后面再出現兩個數字 8 和 9,那么就可以把8更新到d[5], 9更新到d[6],得出LIS的長度為6。
然后應該發現一件事情了:在B中插入數據是有序的,而且是進行替換而不需要挪動——也就是說,我們可以使用二分查找,將每一個數字的插入時間優化到O(logN)~于是算法的時間復雜度就降低到了O(NlogN)~!
代碼:
#include<cstdio> #include<algorithm> using namespace std;int n,s[500005]; //s[i]:存儲長度為i的上升子序列的末尾元素的最小值struct info {int A,B; }a[500005];bool cmp(info i,info j) {return i.A<j.A; }int main() {//freopen("in.txt","r",stdin);scanf("%d",&n);for(int i=1;i<=n;i++){scanf("%d%d",&a[i].A,&a[i].B);}sort(a+1,a+1+n,cmp);int lenth=0,t; //lenth為結果 s[++lenth]=a[1].B;for(int i=2;i<=n;i++){if(a[i].B>s[lenth]){s[++lenth]=a[i].B;}else{t=upper_bound(s+1,s+1+lenth,a[i].B)-s;s[t]=a[i].B;}}printf("%d\n",lenth);return 0; }總結
以上是生活随笔為你收集整理的【训练题】航线设计 | 使用最长上升子序列(LIS)长度的O(nlogn)算法优化的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 一些数据格式化-Eval( )和Dat
- 下一篇: wpf-调用简单的js代码