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

文章目录

  • 一:CF1463F
  • 二:CF1152F2
  • 三:CF1342F

一: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}\rfloorpn出现次数为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,(舍去距离iiimmm的元素,该元素距离i+1i+1i+1变成了m+1m+1m+1,丢掉
iii选,那么这个元素可以放在最后距离≤m\le mm个元素中选了的元素后面,或者放在第一个(目前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
小贪心:一个集合的合并点肯定越靠前越好,集合的和在大于上一个集合的条件下越小越好;这样都是为了增大后面再划分一个集合的概率
预处理集合为sssaaa值和,然后直接转移就没了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;
} 

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/317493.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

C - Log Calculator FZU - 2036

C - Log Calculator FZU - 2036 题意&#xff1a; 已知a和b&#xff0c;满足alog2x&#xff0c;blog2y 求log2(xy) 题解&#xff1a; 因为alog2x&#xff0c;blog2y 所以x2a, yab log2(xy)log2(2a2b)log2(2b(12a-b))blog2(2a-b1) x>32时&#xff0c;blog2(2x1)xb 代码&…

实用帖 | 使用Visual Studio开发.NET Core推荐设置

Visual Studio 一直是.NET程序员的好伙伴&#xff0c;如今虽然.NET Core的开发可以用VSCode等其他环境&#xff0c;但VS依旧有强大的优势。今天来介绍下如何让VS用得更顺手的技能Here We Go!⛷基本配置首先要搞清楚的一点是&#xff0c;并不是所有版本的Visual Studio都支持.NE…

CF1654F-Minimal String Xoration【倍增】

正题 题目链接:https://www.luogu.com.cn/problem/CF1654F 题目大意 给出一个长度为2n2^n2n的字符串sss&#xff08;下标为0∼2n−10\sim 2^n-10∼2n−1&#xff09; 你要找到一个xxx满足tisixorxt_{i}s_{i\ xor\ x}ti​si xor x​&#xff0c;并且ttt的字典序最小。 1≤n≤…

P3642 [APIO2016]烟火表演(左偏树、函数)

解析 感觉是左偏树的神题了. 首先有一个比较显然的结论&#xff0c;一个合法的方案中&#xff0c;两个叶子到它们 lca\text{lca}lca 的距离必须相等. 考虑设计 dp\text{dp}dp &#xff1a; fi,xf_{i,x}fi,x​ 表示 iii 的子树中&#xff0c;所有叶子到它的距离为 xxx 的最小…

2011年全国大学生程序设计邀请赛(福州)

题号题目考点难度APassword table模拟签到BAxial symmetry计算几何代码题CLog Calculatorjava大数DMaximum Value Problem思维题EAnother Postman Problem思维题FPets二分图匹配模板题GTilingHChecker思维题角度很妙&#xff0c;值得多练练IThe Mad MathematicianJSocial Netwo…

概率与期望游戏地图

文章目录前言拾得宝典的玄灵大陆P2911P4316P6154P1297P6862参悟升级的三千界P1365P1850[CF]Crypto Lights打野刷野怪副本[AT-abc215G]Colorful Candies 2前言 其实从去年开始初次接触概率期望后&#xff0c;后面的考试几乎场场见到&#xff1b;而每次都只能敲最原始的暴力&…

.NET Core IdentityServer4实战 第一章-入门与API添加客户端凭据

内容&#xff1a;本文带大家使用IdentityServer4进行对API授权保护的基本策略作者&#xff1a;zara(张子浩) 欢迎分享&#xff0c;但需在文章鲜明处留下原文地址。本文将要讲述如何使用IdentityServer4对API授权保护以及如何携带Token进行访问受保护的API&#xff0c;通过HttpC…

Loj#3077-「2019 集训队互测 Day 4」绝目编诗【结论,虚树,鸽笼原理】

正题 题目链接:https://loj.ac/p/3077 题目大意 给出nnn个点mmm条边的一张简单无向图&#xff0c;求是否存在两个长度相等的简单环。 1≤n≤104,1≤m≤1061\leq n\leq 10^4,1\leq m\leq 10^61≤n≤104,1≤m≤106 解题思路 先考虑一个暴力的做法&#xff0c;我们暴力搜索图上…

蒋金楠:200行代码7个对象《ASP.NET Core框架揭密》苏州.NET俱乐部课程分享

【课程名称】《ASP.NET Core框架揭密》【老师介绍】蒋金楠&#xff0c;同程艺龙机票事业部技术专家&#xff0c;微软最有价值专家&#xff08;MVP&#xff0c;连续12&#xff09;&#xff0c;多本.NET专著作者。博客园Artech&#xff0c;公众号“大内老A”。【课程简介】讲述AS…

CF1137C:Museums Tour(缩点、分层图)

解析 比起CF&#xff0c;更像是主题库的一道题 也是很好的一道题 由于可以在环里转圈&#xff0c;容易想到缩点 但是出现了一个问题&#xff1a;缩点是无法统计一个强连通分量里有哪些环及其大小的 而本题中显然环的大小和d的关系至关重要 考虑转换思路 把每个城市x割成d个点…

H - Checker FZU - 2041

H - Checker FZU - 2041 题意&#xff1a; 一个长度为n的01串&#xff0c;现在能将里面的1移动m次&#xff0c;问最长的连续0是多长 题解&#xff1a; 没想出来&#xff0c;看了其他人代码&#xff0c;就是对于每个0空间进行扩充&#xff0c;然后记录每次扩充的最大值 关键…

2021-4-1 多校省选模拟赛

文章目录考试复盘nmd考试复盘 T1T1T1 我可太喜欢这种不要脑子的莽试题了&#xff01;&#xff01; 考场上猜的结论d>3d>3d>3无解 d1,d2d1,d2d1,d2的填法也试出来了 d3d3d3只试出来了n,mn,mn,m是333的倍数的填法 但是另外的情况就GGGGGG了 最后交上去竟然还敲挂了(…

CF838C-Future Failure【dp,子集卷积】

正题 题目链接:https://www.luogu.com.cn/problem/CF838C 题目大意 一个字符串sss&#xff0c;两个人轮流操作&#xff0c;每次每个人可以选择删掉一个字符或者重排列这个字符串&#xff0c;但是不能出现之前出现过的字符串&#xff0c;不能操作者输。 求有多少个长度为nnn且…

Reactive Extensions 相见恨晚的Rx.Net

何为Reactive Extensions&#xff08;Rx&#xff09;Rx是一个遵循函数式编程的类库&#xff0c;它引用观察者以及迭代器设计模式对可观察对象产生的数据进行异步消费。使用Rx&#xff0c;开发人员将使用LINQ运算符操作异步数据流&#xff0c;并使用调度程序参数化异步数据流中的…

CF819E:Mister B and Flight to the Moon(构造、归纳法)

解析 本题也算看了一半题解吧 看到“数学归纳法”退出来自己推的 这题想到归纳法后面也就简单了 首先&#xff0c;n3和n4的时候显然有解&#xff0c;可以打表 然后考虑在获得n-2的答案时&#xff0c;如何获得n的答案 如果n为奇数&#xff0c;我们可以把(1,n-1,2,n)、(3,n-1,4…

AcWing 1303. 斐波那契前 n 项和

AcWing 1303. 斐波那契前 n 项和 题意: 题解&#xff1a; 这个题解分析的太到位了 代码&#xff1a; #include<iostream> #include<cstring>#define ll long longusing namespace std;const int N3;int n,m; ll f1[N][N]{{1,1,1}}; ll A[N][N]{{0,1,0},{1,1,1…

2021-4-4 省选模拟赛(灯,十字路口,密室逃脱)

文章目录考试复盘A&#xff1a;灯(light)B&#xff1a;十字路口(crossing)C&#xff1a;密室逃脱(escape)考试复盘 第一题分块虽然明显&#xff0c;但是说实话自己没怎么做过分块的题 就不会做大块的处理。。。(;_) 今天听H老说分块可以成替代数据结构的骗分暴力对拍神器 这么…

AGC056E-Cheese【dp】

前言 奶酪可能会长腿&#xff0c;但绝对不会变质 _ _&#xff0c;_ _ _ _ _ _ _ 正题 题目链接:https://atcoder.jp/contests/agc056/tasks/agc056_e 题目大意 有一个长度为nnn的环&#xff0c;第i0.5(0≤i<n)i0.5(0\leq i<n)i0.5(0≤i<n)位置上各有一只老鼠。 然…

张高兴的 .NET Core IoT 入门指南:(二)GPIO 的使用

什么是 GPIOGPIO 是 General Purpose Input Output 的缩写&#xff0c;即“通用输入输出”。 Raspberry Pi 有两行 GPIO 引脚&#xff0c; Raspberry Pi 通过这两行引脚进行一些硬件上的扩展&#xff0c;与传感器进行交互等等。Raspberry Pi B/2B/3B/3B/Zero 引脚图简单的讲&am…

acwing1304. 佳佳的斐波那契

1304. 佳佳的斐波那契 题意&#xff1a; S(n)表示Fibonacci的前n项和mod m T(n)(F12F23F3…nFn)mod m 给n和m&#xff0c;求T(n) 题解&#xff1a; 矩阵快速幂 参考题解 关键在于构造矩阵相乘的形式 代码&#xff1a; #include<bits/stdc.h> #define debug(a,b) pri…