【CF1209E】Rotate Columns【状压DP】【位运算】【贪心】

题意:给一个N×MN \times MN×M的矩阵,可以进行任意多次操作将一列轮换,求每一行的最大值之和的最大值。多组数据。

Easy VersionN≤4N \leq 4N4,M≤100M \leq100M100

Hard VersionN≤12N \leq 12N12,M≤2000M \leq2000M2000

看这数据范围显然是个状压

相当于是每一行只能选一个

f(i,S)f(i,S)f(i,S)表示当前到iii,已经选了SSS的最大值

记忆化搜索一波,暴力转一下

复杂度O(4nn2m)O(4^nn^2m)O(4nn2m) 可以通过Easy Version

    #include <iostream>#include <cstdio>#include <cstring>#include <cctype>using namespace std;int n,m;int a[20][105],dp[105][1<<4];int Move(int s){return (s>>1)|((s&1)<<(n-1));}int dfs(int pos,int s){if (dp[pos][s]) return dp[pos][s];if (s==(1<<n)-1) return 0;if (pos==m) return 0;int &ans=dp[pos][s];for (int t=0;t<(1<<n);t++){if (s&t) continue;int res=dfs(pos+1,s|t);int mx=0;for (int k=0,cur=t;k<n;k++){int sum=0;cur=Move(cur);for (int i=0;i<n;i++)if ((1<<i)&cur)sum+=a[i][pos];mx=max(mx,sum);}ans=max(ans,res+mx);}return ans;}void solve(){memset(dp,0,sizeof(dp));scanf("%d%d",&n,&m);for (int i=0;i<n;i++)for (int j=0;j<m;j++)scanf("%d",&a[i][j]);printf("%d\n",dfs(0,0));}int main(){int T;scanf("%d",&T);while (T--) solve();return 0;}

复杂度里有MMM,过不了HV

考虑如何消掉M

我们发现因为只有NNN行,所以最多只有NNN列选了数

能否快速确定这NNN列?

贪心!

我们把列按最大值从大到小排序,如果前面的一列没选数而后面的选了

那我们完全可以改成前面没有选的,一定更优

因为只能选NNN个数,所以只用考虑最大的NNN

复杂度O(4nn3)O(4^nn^3)O(4nn3)

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
#include <algorithm>
using namespace std;
int n,m;
int a[20][2005],dp[2005][1<<12],ord[2005];
int Move(int s){return (s>>1)|((s&1)<<(n-1));}
int dfs(int pos,int s)
{if (dp[pos][s]) return dp[pos][s];if (s==(1<<n)-1) return 0;if (pos==m) return 0;int &ans=dp[pos][s];for (int t=0;t<(1<<n);t++){if (s&t) continue;int res=dfs(pos+1,s|t);int mx=0;
//		printf("%d",t);for (int k=0,cur=t;k<n;k++){int sum=0;cur=Move(cur);
//			printf("->%d",cur);for (int i=0;i<n;i++)if ((1<<i)&cur)sum+=a[i][ord[pos]];mx=max(mx,sum);}
//		puts("");ans=max(ans,res+mx);}return ans;
}
int mx[2005];
inline bool cmp(const int& a,const int& b){return mx[a]>mx[b];}
void solve()
{memset(dp,0,sizeof(dp));memset(mx,0,sizeof(mx));scanf("%d%d",&n,&m);for (int i=0;i<n;i++)for (int j=0;j<m;j++)scanf("%d",&a[i][j]),mx[j]=max(mx[j],a[i][j]);for (int i=0;i<m;i++) ord[i]=i;sort(ord,ord+m,cmp);m=min(n,m);printf("%d\n",dfs(0,0));
}
int main()
{int T;scanf("%d",&T);while (T--) solve();return 0;
}

然而发现慢得飞起

我们发现:在dp的时候,会花费O(n)O(n)O(n)找出最佳的旋转方案,而这个O(n)O(n)O(n)是在O(4n)O(4^n)O(4n)的基础上的

为什么不预处理呢

然后得到了O(2nn3+4nn)O(2^nn^3+4^nn)O(2nn3+4nn)

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
#include <algorithm>
using namespace std;
int n,m;
int a[20][2005],mem[2005][1<<12],dp[2005][1<<12],ord[2005];
int dfs(int pos,int s)
{if (dp[pos][s]) return dp[pos][s];if (s==(1<<n)-1) return 0;if (pos==m) return 0;int &ans=dp[pos][s];for (int t=0;t<(1<<n);t++)if (!(s&t))ans=max(ans,dfs(pos+1,s|t)+mem[ord[pos]][t]);return ans;
}
int mx[2005];
inline bool cmp(const int& a,const int& b){return mx[a]>mx[b];}
void solve()
{memset(dp,0,sizeof(dp));memset(mx,0,sizeof(mx));memset(mem,0,sizeof(mem));scanf("%d%d",&n,&m);for (int i=0;i<n;i++)for (int j=0;j<m;j++)scanf("%d",&a[i][j]),mx[j]=max(mx[j],a[i][j]);for (int i=0;i<m;i++) ord[i]=i;sort(ord,ord+m,cmp);m=min(n,m);for (int pos=0;pos<m;pos++)for (int s=0;s<(1<<n);s++)for (int k=0;k<n;k++){int sum=0;for (int i=0;i<n;i++)if (s&(1<<i))sum+=a[(i+k)%n][ord[pos]];mem[ord[pos]][s]=max(mem[ord[pos]][s],sum);				}		printf("%d\n",dfs(0,0));
}
int main()
{int T;scanf("%d",&T);while (T--) solve();return 0;
}

然而还是慢得飞起

这玩意还有多组数据

我们发现这一句

		if (!(s&t))

是求和sss没有交集的ttt

能否快速得到这玩意呢?

还真不行

但我们可以换一种实现方式

直接递推实现dp

这样从上一维转移过来只用枚举子集

可以用这句

			for (int t=s;t;t=((t-1)&s))

这样就只枚举了sss的子集

复杂度?

加上枚举sss

一共是

∑i=0nCni2i\sum_{i=0}^{n}C_n^i2^ii=0nCni2i

=∑i=0nCnn−i2i=\sum_{i=0}^{n}C_n^{n-i}2^i=i=0nCnni2i

=(1+2)n=3n=(1+2)^n=3^n=(1+2)n=3n

所以复杂度是O(2nn3+3nn)O(2^nn^3+3^nn)O(2nn3+3nn)

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
#include <algorithm>
using namespace std;
int n,m;
int a[12][2000],mx[2000],pos[2000],mem[2000][1<<12],dp[2000][1<<12];
inline bool cmp(const int& a,const int& b){return mx[a]>mx[b];}
void solve()
{scanf("%d%d",&n,&m);memset(mx,0,sizeof(mx));for (int i=0;i<n;i++)for (int j=0;j<m;j++)scanf("%d",&a[i][j]),mx[j]=max(mx[j],a[i][j]);for (int i=0;i<m;i++) pos[i]=i;sort(pos,pos+m,cmp);m=min(n,m);for (int p=0;p<m;p++)for (int s=0;s<(1<<n);s++){mem[p][s]=0;for (int k=0;k<n;k++){int sum=0;for (int i=0;i<n;i++)if ((1<<i)&s)sum+=a[(i+k)%n][pos[p]];mem[p][s]=max(mem[p][s],sum);}}	memset(dp,0,sizeof(dp));for (int s=0;s<(1<<n);s++) dp[0][s]=mem[0][s];for (int p=1;p<m;p++)for (int s=0;s<(1<<n);s++){for (int t=s;t;t=((t-1)&s))dp[p][s]=max(dp[p][s],dp[p-1][s^t]+mem[p][t]);dp[p][s]=max(dp[p][s],dp[p-1][s]);			}printf("%d\n",dp[m-1][(1<<n)-1]);
}
int main()
{int T;scanf("%d",&T);while (T--) solve();return 0;
}

过于毒瘤

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

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

相关文章

Exceptionless(二) - 使用进阶

作者&#xff1a;markjiang7m2原文地址&#xff1a;https://www.cnblogs.com/markjiang7m2/p/11100563.html官网地址&#xff1a;http://letyouknow.net在上一篇文章Exceptionless - .Net Core开源日志框架中就说到如何对Exceptionless进行本地化部署&#xff0c;不过我也跟大家…

Ozon Tech Challenge 2020 (Div.1 + Div.2) F. Kuroni and the Punishment 随机化

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你nnn个数&#xff0c;每次操作可以选择将某个数1,−11,-11,−1&#xff0c;求最少进行多少次操作使得所有数都为正数且gcd>1gcd>1gcd>1。 思路&#xff1a; 考虑gcd2gcd2gcd2的情况&#xff0…

【洛谷P4841】城市规划【指数型生成函数】【麦克劳林级数】【多项式对数】

传送门 题意&#xff1a;求NNN个点的带标号无向连通简单图的个数。 N≤130000N \leq 130000N≤130000 这个问题的主要矛盾在于连通 这个并不好表示&#xff0c;但可以用这个表示出不要求连通的方案数 由于带标号&#xff0c;先构造答案的EGF f(x)∑i0∞aii!xif(x)\sum_{i0}…

基于 Kong 和 Kubernetes 的 WebApi 多版本解决方案

前言大家好&#xff0c;很久没有写博客了&#xff0c;最近半年也是比较的忙&#xff0c;所以给关注我的粉丝们道个歉。去年和朱永光大哥聊的时候提了一下我们的这个方案&#xff0c;他说让我有空写篇博客讲一下&#xff0c;之前是非常的忙&#xff0c;所以这次趁着有些时间就写…

CodeCraft-20 (Div. 2) D. Nash Matrix 构造 + dfs

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给定一个n∗nn*nn∗n的矩阵&#xff0c;每个点上面都有一个xi,yix_i,y_ixi​,yi​&#xff0c;表示这个点到xi,yix_i,y_ixi​,yi​这个点停下&#xff0c;当xi−1,yi−1x_i-1,y_i-1xi​−1,yi​−1的时候代表…

【HAOI2018】染色【反向二项式反演】【NTT卷积】

传送门 题意&#xff1a;NNN个位置染MMM种颜色&#xff0c;恰好出现SSS次的颜色数量恰好为kkk时的愉悦度为wkw_kwk​,求所有方案的愉悦度之和。对100453580910045358091004535809取模。 N≤1e7N \leq 1e7N≤1e7,M≤1e5M \leq 1e5M≤1e5,S≤150S \leq 150S≤150 本题的恶心之处…

从严治码-别人在项目中下毒,我该怎么治?

01 从软考说起从4月份开始&#xff0c;由于备考《系统集成项目管理工程师》的原因&#xff0c;博客没有持续更新&#xff0c;在上半年考试结束之后&#xff0c;又对项目进行了一些收尾的工作。下面就这段时间的学习作一个记录和总结吧。在学习的过程中&#xff0c;提炼了一些自…

CodeCraft-20 (Div. 2) C. Primitive Primes 思维 + 数论

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你两个长度分别为n,mn,mn,m的多项式&#xff0c;将他们乘起来&#xff0c;问系数modp0\bmod p 0modp0的项的指数是多少&#xff0c;两个多项式所有项的系数gcd1gcd1gcd1。 n,m<1e6n,m<1e6n,m<1e6…

【LOJ166】拉格朗日插值2【拉格朗日插值】【NTT卷积】

传送门 题意&#xff1a;给定n,m,f(0),f(1),......,f(n)n,m,f(0),f(1),... ...,f(n)n,m,f(0),f(1),......,f(n)&#xff0c;求f(m),f(m1),......f(mn)f(m),f(m1),... ...f(mn)f(m),f(m1),......f(mn) 模998244353998244353998244353 n≤100000,m≤1e8,n<mn \leq 100000,m \…

学习MVVM设计模式后第一次用于生产

WPF的MVVM设计模式从winform转变到WPF的过程&#xff0c;难点主要还是在MVVM的设计模式。当然&#xff0c;如果依然采用winform的涉及方式&#xff0c;在每个控件背后绑定事件的方式运用在wpf中&#xff0c;依然可行&#xff0c;但是假如GUI改版&#xff0c;其背后绑定的特别为…

剑指 Offer 14- II. 剪绳子 II

给你一根长度为 n 的绳子&#xff0c;请把绳子剪成整数长度的 m 段&#xff08;m、n都是整数&#xff0c;n>1并且m>1&#xff09;&#xff0c;每段绳子的长度记为 k[0],k[1]...k[m - 1] 。请问 k[0]*k[1]*...*k[m - 1] 可能的最大乘积是多少&#xff1f;例如&#xff0c;…

Educational Codeforces Round 108 (Rated for Div. 2) D. Maximum Sum of Products 思维 + dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你两个长度为nnn的数组a,ba,ba,b&#xff0c;你可以至多反转一段连续区间&#xff0c;求∑i1nai∗bi\sum _{i1}^n a_i*b_i∑i1n​ai​∗bi​最大是多少。 n<5e3n<5e3n<5e3 思路&#xff1a; 首…

【CF1215E】Marbles【状压DP】

传送门 题意&#xff1a;给一个长为NNN的序列aaa&#xff0c;每次操作交换两个相邻位置&#xff0c;求最少操作次数使得所有相同的值连成一片。 N≤400000N \leq 400000N≤400000,ai≤20a_i \leq20ai​≤20 我们发现aia_iai​很小&#xff0c;盲猜单独考虑 我们重新确认一个…

netcore mvc快速开发系统(菜单,角色,权限[精确到按钮])开源

基于netcore2.0 mvc 开发的 快速搭建具有如下特色的后台管理系统用户管理菜单管理角色管理权限管理[精确到按钮]&#xff09;代码生成器代码克隆到本地 用vs2017或以上版本 打开工程。项目结构如下&#xff1a;找到DbModel下面的初始化db脚本里面包含4张表的schema和初始化数据…

剑指 Offer 25. 合并两个排序的链表

输入两个递增排序的链表&#xff0c;合并这两个链表并使新链表中的节点仍然是递增排序的。 示例1&#xff1a; 输入&#xff1a;1->2->4, 1->3->4 输出&#xff1a;1->1->2->3->4->4 限制&#xff1a; 0 < 链表长度 < 1000 思路&#xff1a;和…

Codeforces Round #720 (Div. 2) C. Nastia and a Hidden Permutation 交互

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个序列ppp长度nnn&#xff0c;每次可以执行两个种询问&#xff1a; t1max(min(x,pi),min(x1,pj))t1\ \ max(min(x,p_i),min(x1,p_j))t1 max(min(x,pi​),min(x1,pj​)) t2min(max(x,pi),max(x1,pj))t…

【NOIP2018】赛道修建【二分】【树形dp】【multiset】【贪心】

传送门 题意&#xff1a;给一棵带边权的树&#xff0c;求MMM条没有公共边的路径的最小长度的最大值。 N≤50000N \leq50000N≤50000 抛开NOIP不谈&#xff0c;其实这题本身出的很好 显然二分 问题转化成了“最多可以找多少条长度不小于kkk的路径” 递归处理 对于每个结点…

asp.net core 系列之Startup

这篇文章简单记录 ASP.NET Core中 &#xff0c;startup类的一些使用。一.前言在 Startup类中&#xff0c;一般有两个方法&#xff1a;ConfigureServices 方法: 用来配置应用的 service 。 Configure 方法&#xff1a;创建应用的请求处理管道它们都在应用启动时&#xff0c;被AS…

可持久化Splay 学习笔记

可持久化Splay是怎么回事呢&#xff1f;Splay相信大家都很熟悉&#xff0c;但是可持久化Splay是怎么回事呢&#xff0c;下面就让小编带大家一起了解吧。   可持久化Splay&#xff0c;其实就是将Splay持久化一下&#xff0c;大家可能会很惊讶Splay怎么可以持久化呢&#xff1f…

译 | .NET Core 基础架构进化之路(一)

原文&#xff1a;Matt Mitchell翻译&#xff1a;Edi Wang随着 .NET Core 3.0 Preview 6 的推出&#xff0c;我们认为简要了解一下我们基础设施系统的历史以及过去一年左右所做的重大改进会很有用。如果您对构建基础结构感兴趣&#xff0c;或者想要了解我们如何构建与 .NET Core…