【bzoj1044】[HAOI2008]木棍分割 二分+dp
題目描述
有n根木棍, 第i根木棍的長度為Li,n根木棍依次連結了一起, 總共有n-1個連接處. 現在允許你最多砍斷m個連接處, 砍完后n根木棍被分成了很多段,要求滿足總長度最大的一段長度最小, 并且輸出有多少種砍的方法使得總長度最大的一段長度最小. 并將結果mod 10007。。。
輸入
輸入文件第一行有2個數n,m.接下來n行每行一個正整數Li,表示第i根木棍的長度.n<=50000,0<=m<=min(n-1,1000),1<=Li<=1000.
輸出
輸出有2個數, 第一個數是總長度最大的一段的長度最小值, 第二個數是有多少種砍的方法使得滿足條件.
樣例輸入
3 2
1
1
10
樣例輸出
10 2
題解
二分+dp
第一問即 noip2015跳石頭 。。。一眼二分,然后看不滿足條件時就切一刀,判斷是否小于m。
第二問求方案數,很顯然是個dp。
設$f[i][j]$表示前$i$個分了$j$段的方案數,那么狀態轉移方程應該為$f[i][j]=\sum\limits_{len(t+1,i)\le ans1}f[t][j-1]$,邊界條件$f[0][0]=1$,其中$len(a,b)表示$[a,b]$所有木棍的長度總和。
可以發現$t$的取值范圍是一段連續的單調的區間,因此可以用類似雙指針的方法掃出$t$的取值左端點。然后$\sum$又可以使用前綴和維護,這樣時間復雜度就降為了$O(nm)$。
然而這樣還會炸空間。。。
因此使用滾動數組就好了,顯然第二維是可以滾動的,因此先枚舉第二維,滾動一下就好了。
#include <cstdio> #include <algorithm> #define N 50010 #define mod 10007 using namespace std; int n , m , a[N] , sl[N] , f[2][N] , sum[2][N]; bool judge(int mid) {int i , now = 0 , cnt = 0;for(i = 1 ; i <= n ; i ++ ){if(now + a[i] > mid) now = 0 , cnt ++ ;now += a[i];}return cnt <= m; } int main() {int i , j , l = 0 , r = 0 , mid , ans = -1 , p = 0 , ret = 0 , d;scanf("%d%d" , &n , &m);for(i = 1 ; i <= n ; i ++ ) scanf("%d" , &a[i]) , l = max(l , a[i]) , r += a[i] , sl[i] = sl[i - 1] + a[i];while(l <= r){mid = (l + r) >> 1;if(judge(mid)) ans = mid , r = mid - 1;else l = mid + 1;}printf("%d " , ans);for(i = 0 ; i <= n ; i ++ ) sum[0][i] = 1;for(i = d = 1 ; i <= m + 1 ; i ++ , d ^= 1){sum[d][0] = p = 0;for(j = 1 ; j <= n ; j ++ ){while(sl[j] - sl[p] > ans) p ++ ;f[d][j] = sum[d ^ 1][j - 1];if(p) f[d][j] = (f[d][j] - sum[d ^ 1][p - 1] + mod) % mod;sum[d][j] = (sum[d][j - 1] + f[d][j]) % mod;}ret = (ret + f[d][n]) % mod;}printf("%d\n" , ret);return 0; }?
?
轉載于:https://www.cnblogs.com/GXZlegend/p/7491510.html
總結
以上是生活随笔為你收集整理的【bzoj1044】[HAOI2008]木棍分割 二分+dp的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 创建基于Web的实时系统
- 下一篇: 【自动化__持续集成】___java__