BZOJ 5161 最長上公升子串行 狀壓DP

2022-05-09 16:26:39 字數 1101 閱讀 9324

題意

求乙個\(1\sim n\)的排列lis的期望長度,\(n\leq 28\)

題解

考慮樸素的lis:\(f[i] = min(f[j]) + 1\)

記\(mx[i]\)為\(f\)的字首最大值,那麼可以得到乙個性質\(mx[i + 1] \in [mx[i], mx[i] + 1]\)

對\(mx\)陣列進行差分,則差分陣列只有\(01\),可以狀壓

由於\(mx[1] - mx[0]=1\),從第二位開始狀壓

然後考慮從\(1\sim i\)的排列推到\(1\sim i+1\)的排列,\(1\sim i\)的差分陣列為\(s\)。把i插入第\(j\)(\(1\sim i+1\))個數前面:

若\(j=1\):新狀態為\(s << 1\),即\(mx[2]\)為\(0\),而其餘不改

若\(j>1\):\(mx[1]\) 到 \(mx[j - 1]\)不變,\(mx[j]\)到\(mx[i]\)整體右移一位,新的\(mx[j]\)為1,但要把新的\(mx[j]\)後面第乙個\(1\)刪去(因為\(j\)處答案變大了,所以該點和新答案一樣大,差分陣列對應\(0\))

注意實現的時候第\(k\)個位置對應\(1<<(k-2)\)

打表程式:

static int dp[2][134217728], ans;

for(int n = 1; n <= 28; n ++) }}

for(int i = 0; i < (1 << (n - 1)); i ++)

upd(ans, 1ll * dp[r][i] * (__builtin_popcount(i) + 1) % mo);

int fac = 1;

for(int i = 2; i <= n; i ++) fac = 1ll * fac * i % mo;

printf("%d, ", 1ll * ans * qpow(fac, mo - 2) % mo);

}

提交程式:

#include const int qwq = ;

int main()

最長上公升子串行

問題描述 乙個數的序列bi,當b1 b2 bs的時候,我們稱這個序列是上公升的。對於給定的乙個序列 a1,a2,an 我們可以得到一些上公升的子串行 ai1,ai2,aik 這裡1 i1 i2 ik n。比如,對於序列 1,7,3,5,9,4,8 有它的一些上公升子串行,如 1,7 3,4,8 等等...

最長上公升子串行

最長上公升子串行問題是各類資訊學競賽中的常見題型,也常常用來做介紹動態規劃演算法的引例,筆者接下來將會對poj上出現過的這類題目做乙個總結,並介紹解決lis問題的兩個常用 演算法 n 2 和 nlogn 問題描述 給出乙個序列a1,a2,a3,a4,a5,a6,a7.an,求它的乙個子串行 設為s1...

最長上公升子串行

最長上公升子串行問題 給出乙個由n個數組成的序列x 1.n 找出它的最長單調上公升子串行。即求最大的m和a1,a2 am,使得a1動態規劃求解思路分析 o n 2 經典的o n 2 的動態規劃演算法,設a i 表示序列中的第i個數,f i 表示從1到i這一段中以i結尾的最長上公升子串行的長度,初始時...