一：CF1463F

Max Correct Set
有一个结论：以x+yx+yx+y为周期排列填充一定是不劣于最后的答案的
令p=x+y,r=n%pp=x+y,r=n\%pp=x+y,r=n%p
⌊np⌋\lfloor\frac{n}{p}\rfloor⌊pn⌋出现次数为ttt次，则rrr的出现次数为t+1t+1t+1，预处理一下即可

设dpi,sdp_{i,s}dpi,s表示一段ppp中的某个点iii为断点，最后max(x,y)max(x,y)max(x,y)位的状态为jjj的答案数

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
int n, x, y;
int dp[2][1 << 22], val[50];int main() {scanf( "%d %d %d", &n, &x, &y );int p = x + y, r = n % p;for( int i = 0;i < p;i ++ )val[i] = n / p + ( i < r );int S = 1 << max( x, y );memset( dp, -1, sizeof( dp ) );dp[1][0] = 0;for( int i = 0;i < p;i ++ ) {int now = i & 1;memset( dp[now], -1, sizeof( dp[now] ) );for( int j = 0;j < S;j ++ )if( ~ dp[now ^ 1][j] ) {dp[now][j << 1 & S - 1] = max( dp[now][j << 1 & S - 1], dp[now ^ 1][j] );if( ! ( 1 << x - 1 & j ) && ! ( 1 << y - 1 & j ) )dp[now][j << 1 & S - 1 | 1] = max( dp[now][j << 1 & S - 1 | 1], dp[now ^ 1][j] + val[i] );}}int ans = 0;for( int i = 0;i < S;i ++ )ans = max( ans, dp[( p - 1 ) & 1][i] );printf( "%d\n", ans );return 0;
}

二：CF1152F2

Neko Rules the Catniverse (Large Version)
简单版：
设dpi,j,sdp_{i,j,s}dpi,j,s：表示前iii个数，从中选了jjj个，离iii最近的mmm个元素是否选择的状态sss
考虑转移到i+1i+1i+1
iii不选，直接转移到dpi+1,j,s>>1dp_{i+1,j,s>>1}dpi+1,j,s>>1，（舍去距离iii为mmm的元素，该元素距离i+1i+1i+1变成了m+1m+1m+1，丢掉
iii选，那么这个元素可以放在最后距离≤m\le m≤m个元素中选了的元素后面，或者放在第一个（目前iii是最大的，一定可以放在最前面

#include <cstdio>
#define int long long
#define mod 1000000007
int n, k, m;
int dp[100005][14][1 << 4];signed main() {scanf( "%lld %lld %lld", &n, &k, &m );dp[0][0][0] = 1;for( int i = 0;i < n;i ++ )for( int j = 0;j <= k;j ++ )for( int s = 0;s < ( 1 << m );s ++ ) {dp[i + 1][j][s >> 1] = ( dp[i + 1][j][s >> 1] + dp[i][j][s] ) % mod;dp[i + 1][j + 1][s >> 1 | 1 << m - 1] = ( dp[i + 1][j + 1][s >> 1 | 1 << m - 1] + dp[i][j][s] * ( __builtin_popcount( s ) + 1 ) % mod ) % mod;}int ans = 0;for( int i = 0;i < ( 1 << m );i ++ )ans = ( ans + dp[n][k][i] ) % mod;printf( "%lld\n", ans );return 0;
}

加强版是建立在发现简单版的第二维和第三维大小一共在210210210左右
是可以用矩阵快速幂优化的，把第二维和第三维压在一起重新编号代替

#include <cstdio>
#include <cstring>
#define int long long
#define mod 1000000007
int n, k, m, cnt;
int id[15][20];struct matrix {int c[210][210];//2^4*13=208matrix() {memset( c, 0, sizeof( c ) );}matrix operator * ( matrix &v ) {matrix ans;for( int i = 1;i <= cnt;i ++ )for( int j = 1;j <= cnt;j ++ )for( int k = 1;k <= cnt;k ++ )ans[i][j] = ( ans[i][j] + c[i][k] * v[k][j] % mod ) % mod;return ans;}int *operator [] ( int x ) { return c[x]; }
}g, f;signed main() {scanf( "%lld %lld %lld", &n, &k, &m );for( int i = 0;i <= k;i ++ )for( int j = 0;j < ( 1 << m );j ++ )id[i][j] = ++ cnt;for( int i = 0;i <= k;i ++ )for( int s = 0;s < ( 1 << m );s ++ ) {//          int t = s << 1 & ( ( 1 << m ) - 1 );
//          ( g[id[i][s]][id[i][t]] += 1 ) %= mod;
//          if( i < k ) ( g[id[i][s]][id[i + 1][t | 1] += __builtin_popcount( s ) + 1 ) %= mod;
//两种写法均是正确的( g[id[i][s]][id[i][s >> 1]] += 1 ) %= mod;if( i < k ) ( g[id[i][s]][id[i + 1][s >> 1 | 1 << m - 1]] += __builtin_popcount( s ) + 1 ) %= mod;}f[1][1] = 1;while( n ) {if( n & 1 ) f = f * g;g = g * g;n >>= 1;}int ret = 0;for( int i = 0;i < ( 1 << m );i ++ )ret = ( ret + f[1][id[k][i]] ) % mod;printf( "%lld\n", ret );return 0;
}

三：CF1342F

Make It Ascending
设dpi,j,sdp_{i,j,s}dpi,j,s：表示划分了iii个集合，最后一个集合的合并点在jjj，最后一个集合选择的点状态为sss
小贪心：一个集合的合并点肯定越靠前越好，集合的和在大于上一个集合的条件下越小越好；这样都是为了增大后面再划分一个集合的概率
预处理集合为sss的aaa值和，然后直接转移就没了o(=·ω·=)m

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
#define INF 0x7f7f7f7f
int T, n;
int a[16], sum[1 << 15], id[16];
pair < int, int > pre[16][16][1 << 15];
int f[16][16][1 << 15];void opt( int x ) {for( int i = x;i < n;i ++) id[i] --;
}void print( int i, int pos, int s ) {if( ! s ) return;int lst = pre[i][pos][s].first, t = lst ^ s;for( int k = 0;k < n;k ++ )if( k != ( pos - 1 ) && ( 1 << k & t ) ) {printf( "%d %d\n", id[k], id[pos - 1] );opt( k );}print( i - 1, pre[i][pos][s].second, lst );
}void solve() {scanf( "%d", &n );for( int i = 0;i < n;i ++ )scanf( "%d", &a[i] );int S = 1 << n;for( int i = 1;i < S;i ++ )sum[i] = a[__builtin_ctz( i )] + sum[i & ( i - 1 )];for( int i = 0;i <= n;i ++ )for( int j = 0;j <= n;j ++ )memset( f[i][j], 0x7f, sizeof( int ) * S );f[0][0][0] = 0;for( int i = 0;i < n;i ++ )for( int s = 0;s < S;s ++ )for( int j = 0;j < n;j ++ ) {if( f[i][j][s] == INF ) continue;int t = ( S - 1 ) ^ s;for( int k = t;k;k = ( k - 1 ) & t ) {if( sum[k] <= f[i][j][s] || ! ( k >> j ) ) continue;int pos = j + __builtin_ctz( k >> j ) + 1;if( f[i + 1][pos][s | k] > sum[k] ) {f[i + 1][pos][s | k] = sum[k];pre[i + 1][pos][s | k] = make_pair( s, j );}}}for( int i = 0;i < n;i ++ ) id[i] = i + 1;for( int i = n;i;i -- )for( int j = 1;j <= n;j ++ )if( f[i][j][S - 1] != INF ) {printf( "%d\n", n - i );print( i, j, S - 1 );return;}
}int main() {scanf( "%d", &T );while( T -- ) solve();return 0;
}

一二三系列之状压DP——Max Correct Set（一）Neko Rules the Catniverse (Large Version)（二）Make It Ascending（三）相关推荐

HDU 4539郑厂长系列故事――排兵布阵（状压DP）
HDU 4539 郑厂长系列故事――排兵布阵基础的状压DP,首先记录先每一行可取的所哟状态(一行里互不冲突的大概160个状态), 直接套了一个4重循环居然没超时我就呵呵了 1 //#pragma ...
CF1463F Max Correct Set（取小样法+状压 DP）
CF1463F Max Correct Set 要求选出集合 \(U=\{1,2,3,\dots,n\}\) 的一个子集 \(S\),满足:如果 \(a \in S\) 并且 \(b \in S\), ...
YYHS-蜀传之单刀赴会（梦回三国系列T2）（最短路+状压dp）
题目描述 [题目背景] 公元215年,刘备取益州,孙权令诸葛瑾找刘备索要荆州.刘备不答应,孙权极为恼恨,便派吕蒙率军取长沙.零陵.桂阳三郡.长沙.桂阳蜀将当即投降.刘备得知后,亲自从成都赶到公安(今湖 ...
POJ 1038 Bugs Integrated Inc （复杂的状压DP）
\(POJ~1038~~*Bugs~Integrated~Inc:\) (复杂的状压DP) \(solution:\) 很纠结的一道题目,写了大半天,就想练练手,结果这手生的.其实根据之前那道炮兵阵地 ...
hdu 4778 Gems Fight! 状压dp
转自wdd :http://blog.csdn.net/u010535824/article/details/38540835 题目链接:hdu 4778 状压DP 用DP[i]表示从i状态选到结束得 ...
P2340 奶牛会展（状压dp）
P2340 奶牛会展题目背景奶牛想证明它们是聪明而风趣的.为此,贝西筹备了一个奶牛博览会,她已经对N 头奶牛进行了面试,确定了每头奶牛的智商和情商. 题目描述贝西有权选择让哪些奶牛参加展览.由 ...
UVa 11825 (状压DP) Hackers' Crackdown
这是我做状压DP的第一道题,状压里面都是用位运算来完成的,只要耐下心来弄明白每次位运算的含义,还是容易理解的. 题意: 有编号为0~n-1的n台服务器,每台都运行着n中服务,每台服务器还和若干台其他服 ...
Codeforces Gym 100676G Training Camp 状压dp
http://codeforces.com/gym/100676 题目大意是告诉你要修n门课,每门课有一个权值w[i], 在第k天修该课程讲获得k*w[i]的学习点数,给出了课程与先修课程的关系,要修 ...
POJ 1185 炮兵阵地【状压DP】
<题目链接> 题目大意: 司令部的将军们打算在N*M的网格地图上部署他们的炮兵部队.一个N*M的地图由N行M列组成,地图的每一格可能是山地(用"H" 表示),也可能是平 ...

一二三系列之状压DP——Max Correct Set（一）Neko Rules the Catniverse (Large Version)（二）Make It Ascending（三）

文章目录

一：CF1463F

二：CF1152F2

三：CF1342F

一二三系列之状压DP——Max Correct Set（一）Neko Rules the Catniverse (Large Version)（二）Make It Ascending（三）相关推荐

最新文章

热门文章