SOSDP顾名思义就是,救命啊是DP Sum over Subsets(SOS)DP
本文约等于https://codeforces.com/blog/entry/45223的机翻中文版
问题引入
给一个2^n长度的数组A,现在对于任意x要预处理出函数F(x)的返回值。
F(x)的定义:SUM(A[i] | x & i == i)
即 i 的二进制表示被x包含,F(x)返回所有满足条件 i 的A[ i ]总和。
解法
暴力 O(4^n)
//枚举每一种x
for(int mask = 0;mask < (1<<N); ++mask){
//枚举每一个i
for(int i = 0;i < (1<<N); ++i){
if((mask&i) == i){
F[mask] += A[i];
}
}
}
这种解法枚举了很多无用的i,对于给定的x,可以用更高效的方法枚举其包含的二进制数。
优化后的暴力O(3^n)
//依旧枚举所有的x
for (int mask = 0; mask < (1<<n); mask++){
F[mask] = A[0];
// 枚举了所有有效的i
for(int i = mask; i > 0; i = (i-1) & mask){
F[mask] += A[i];
}
}
复杂度证明:枚举i含有k个1位,则有C(n,k)种可能,每种可能有2^k种情况。
暴力的做法最大的浪费在于把x的二进制位混在一起枚举。
DP
要把枚举的x,按照二进制拆分。
设计dp状态:
dp[被数字x包含][且最右边y位和x相同] = 的下标贡献总和
初始状态:
dp[i][0] = A[i];
转移方程:
for(int i = 1 ; i < (1 << maxn) ; i++){