Codeforces Round #622 (Div. 2) D. Happy New Year 状压dp

传送门

文章目录

  • 题意:
  • 思路:

题意:

在这里插入图片描述
n≤1e5,m≤1e9,k≤8.n\le 1e5,m\le 1e9,k\le 8.n1e5,m1e9,k8.

思路:

注意到题目中保证了每个孩子至多收到kkk个,且k≤8k\le 8k8,注意到这是题目保证的,并不是说这个孩子能收到很多,但是收到kkk个就接受不了别的了,我理解错题意想了一晚上。
看到kkk这么小,自然的想到状压了,复杂度O(n∗2k)O(n*2^k)O(n2k)很完美。
发现区间很大,自然想到离散化。通过离散化, 我们可以将其分成若干个左闭右开的区间。这个时候就可以定义f[i][j]f[i][j]f[i][j]表示选到了第iii个区间,且该区间的状态为jjjjjj是一个长度为kkk的二进制,某一位iii111表示第iii条线覆盖了这个区间。由于区间是左闭右开的区间,所以左端点才有信息存储,下面分情况讨论:
设当前状态为jjj,其含有的111的个数为cntcntcnt,当前这段区间在第kkk个。
当前需要加上一段区间:
(1)(1)(1)如果jjj的第k−1k-1k1位已经存在,那么应该从前一个状态不含这一位转移过来,即 f[i][j]=f[i−1][jxor1<<k]+len∗(cntmod2)f[i][j]=f[i-1][j \ \ xor \ \ 1<<k]+len*(cnt\bmod 2)f[i][j]=f[i1][j  xor  1<<k]+len(cntmod2)
(2)(2)(2)如果jjj的第k−1k-1k1位不存在,那么应该从上一个jjj的状态转移,即f[i][j]=f[i−1][j]+len∗(cntmod2)f[i][j]=f[i-1][j]+len*(cnt\bmod 2)f[i][j]=f[i1][j]+len(cntmod2)
当前需要减去一段区间:
(1)(1)(1)如果jjjk−1k-1k1位已经存在,由于我们是要删掉他,所以应该将他置为−INF-INFINF,即f[i][j]=−INFf[i][j]=-INFf[i][j]=INF
(2)(2)(2)如果jjjk−1k-1k1位不存在,那么从上一个状态的jjjjxor1<<kj\ \ xor\ \ 1<<kj  xor  1<<k二者大的一个转移过来,再加上len∗(cntmod2)len*(cnt\bmod 2)len(cntmod2)即可,即f[i][j]=max(f[i−1][j],f[i−1][jxor1<<k])+len∗(cntmod2)f[i][j]=max(f[i-1][j],f[i-1][j \ \ xor \ \ 1<<k])+len*(cnt\bmod 2)f[i][j]=max(f[i1][j],f[i1][j  xor  1<<k])+len(cntmod2)

第一维直接滚动一下就好了,或者直接开一维,让后注意一下遍历的顺序就好啦,具体看代码吧。

// Problem: D. Happy New Year
// Contest: Codeforces - Codeforces Round #622 (Div. 2)
// URL: https://codeforces.com/contest/1313/problem/D
// Memory Limit: 512 MB
// Time Limit: 2000 ms
// 
// Powered by CP Editor (https://cpeditor.org)//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid (tr[u].l+tr[u].r>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;const int N=400010,mod=1e9+7,INF=0x3f3f3f3f;
const double eps=1e-6;int n,m,k;
PII p[N];
vector<int>v;
int st[N];
LL f[2][(1<<8)+10];int find(int x)
{return lower_bound(v.begin(),v.end(),x)-v.begin();
}int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);int tot=0,op=0;scanf("%d%d%d",&n,&m,&k);for(int i=1;i<=n;i++){int l,r; scanf("%d%d",&l,&r);v.pb(l); v.pb(r+1);p[++tot]={l,i}; p[++tot]={r+1,-i};}//for(int i=1;i<(1<<k);i++) f[0][i]=-INF;memset(f,-0x3f,sizeof(f));f[op][0]=0;sort(p+1,p+1+tot);LL ans=0;for(int i=1;i<=tot;i++) {op^=1;int len=p[i+1].X-p[i].X;if(i==tot) len=0;int id=p[i].Y,pos;if(id>0) {for(int j=0;j<k;j++) if(st[j]==0) { pos=j; st[j]=id; break; }for(int j=0;j<(1<<k);j++)if(j>>pos&1) f[op][j]=f[op^1][j^1<<pos]+len*__builtin_parity(j);else f[op][j]=f[op^1][j]+len*__builtin_parity(j);}else {for(int j=0;j<k;j++) if(st[j]==-id) { pos=j; st[j]=0; break; }for(int j=0;j<(1<<k);j++)if(j>>pos&1) f[op][j]=-INF;else f[op][j]=max(f[op^1][j],f[op^1][j^1<<pos])+len*__builtin_parity(j);}}cout<<f[op][0]<<endl;/*for(int i=1;i<=tot;i++){int len;if(i!=tot) len=p[i+1].X-p[i].X;else len=0;int id=p[i].Y,pos;if(id>0){for(int i=0;i<k;i++) if(st[i]==0) { pos=i; st[i]=id; break; }for(int i=(1<<k)-1;i>=0;i--) if(i>>pos&1) f[i]=f[i^1<<pos]+len*__builtin_parity(i);else f[i]=f[i]+len*__builtin_parity(i);}else {for(int i=0;i<k;i++) if(st[i]==-id) { pos=i; st[i]=0; break; }for(int i=0;i<(1<<k);i++) if(i>>pos&1) f[i]=INF;else f[i]=max(f[i],f[i^1<<pos])+len*__builtin_parity(i);}} printf("%lld\n",f[0]);*/return 0;
}
/**/

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

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

相关文章

2021牛客第四场-E Tree Xor-线段树区间异或

https://ac.nowcoder.com/acm/contest/11255/E 题目大意&#xff1a;有n个点&#xff0c;每一点有一个区间[l,r]范围的权值可能&#xff0c;有n-1条边&#xff0c;边上的权值为w,ww(u)^w(v),问你满足权值范围和边权值的方案数。 思路&#xff1a;我们把树根的值设为0&#xff…

【SPOJ2666】QTree4【链分治】

传送门 题意&#xff1a;给一棵带边权的树&#xff0c;每个点开始时为白色&#xff0c;维护两种操作&#xff1a; 1.改变一个点的颜色&#xff08;白变黑&#xff0c;黑变白&#xff09; 2.询问最远的两个白点之间的距离 树分治国集论文 链分治的本质其实就是树链剖分。它们…

2021牛客第五场 I.Interval Queries-回滚莫队

如图&#xff1a;一道很裸的回滚莫队&#xff0c;注意加入的操作和回滚的操作就好了。 #include <cstdio> #include <cstring> #include <string> #include <cmath> #include <iostream> #include <algorithm> #include <queue> #i…

Linux下搭建asp.net运行环境

最近有个项目&#xff0c;是在Windows平台下开发的&#xff0c;需要把 asp.net web应用移植到 CentOS下&#xff0c;甚是头疼&#xff1b;翻阅资料&#xff0c;发现Jexus是个可行的方案&#xff0c;下面是官方对Jexus的定义&#xff1a;什么是JexusJexus是一款Linux平台上的高性…

Codeforces Round #736 (Div. 2)E. The Three Little Pigs-长除法求多项式系数

https://codeforces.com/contest/1549/problem/E 题目大意自己悟吧&#xff0c;不解释了。 这题在赛时我想的是fft&#xff0c;但是题解说不用fft就能求出多项式的系数&#xff0c;fft在这题时间复杂度有点高。 使用多项式长除法就能求出系数&#xff0c;k^x的系数就是抓x个p…

在 Asp.Net Core WebAPI 中防御跨站请求伪造攻击

什么是跨站请求伪造跨站请求伪造&#xff08;英语&#xff1a;Cross-site request forgery&#xff09;&#xff0c;也被称为 one-click attack 或者 session riding&#xff0c;通常缩写为 CSRF 或者 XSRF&#xff0c; 是一种挟制用户在当前已登录的Web应用程序上执行非本意的…

HDU2021多校第五天 1009 Array-树状数组求二阶前缀和

https://acm.hdu.edu.cn/showproblem.php?pid7020 题目大意&#xff1a;求出区间众数严格大于区间一半大小的子区间的数量。 思路&#xff1a;先对每一个数求出对应的位置排列&#xff0c;S[j]代表1-有多少个i&#xff0c;易得s[r]-s[l]>r-l-(s[r]-s[l])的时候子区间{l1&a…

C# 8 新增小功能

尽管 C# 8 应该会在今年发布&#xff0c;并且 C# 8.x 和 9 的路线图也开始形成&#xff0c;但是&#xff0c;微软正在继续审查下一个版本的特性。目标类型表达式假设我们有 A、B 和 C 三种类型&#xff0c;其中类型 B 和 C 是类型 A 的子类型。如果我们使用如下所示的语句&…

2021HDU多校第五场1004 Another String-双指针

https://acm.hdu.edu.cn/showproblem.php?pid7015 题意&#xff1a;对于每次以i来分割字符串&#xff0c;计算两个分割串的子串的k-匹配的数量。 思路&#xff1a;当我们计算到ans[i]是&#xff0c;我们可以当成从ans[i-1]加上【1-i】字符串和【i1&#xff0c;n】的k-匹配数量…

C# 默认接口方法更新完成,很多细节问题尚待解决

随着对默认接口方法的支持越来越接近完成&#xff0c;一些潜在的问题被提了出来。虽然已经完成了很多工作&#xff0c;但这是一个复杂的特性&#xff0c;许多细节问题还没有解决。但首先&#xff0c;这里有一些已解决的问题。接口允许使用 static 和 const 字段了。除 和! 之外…

2021HDU多校第九场1008HDU7073 Integers Have Friends 2. 随机化

HUD地址&#xff1a;https://acm.hdu.edu.cn/showproblem.php?pid7073 题目大意&#xff1a;选择数组中最多的数&#xff0c;使得他们模m同余&#xff08;m>2&#xff09;。求最大的数量。 思路&#xff1a;然后我们全选奇数或者全选偶数的话&#xff0c;那么我们的答案是…

一个引发程序员们干架的问题

这里是Z哥的个人公众号每周五早8点 按时送达当然了&#xff0c;也会时不时加个餐&#xff5e;我的第「77」篇原创敬上在一个分布式系统的开发团队中&#xff0c;有一些问题是很容易产生程序员之间矛盾的。其中之一就是「业务归属」&#xff0c;就是当新加/修改一个业务的时候&a…

微服务探索与实践—服务注册与发现

前言微服务从大规模使用到现在已经有很多年了&#xff0c;从之前的探索到一步步的不断完善与成熟&#xff0c;微服务已经成为众多架构选择中所必须面对的一个选项。服务注册与发现是相辅相成的&#xff0c;所以一般会合起来思索。其依托组件有很多&#xff0c;比如Zookeeper,Co…

微服务探索与实践—总述

背景软件开发是一个不断发展的过程&#xff0c;从当初的面向过程为主到如今的面向对象的开发&#xff0c;软件开发者不断探索与实践更加符合时代发展要求的开发模式与架构思想&#xff0c;而这&#xff0c;也在极大程度上提高了软件开发的效率。微服务是一种架构模式或者说是架…

这周,全球首个IT技术全中文免费学习平台诞生了!太惊艳!

本周三的Developer Wednesday为大家揭秘了Microsoft Learn网站的全新学习模式一时间Microsoft Learn圈粉无数如果你正在为错过直播而暴风哭泣等一等&#xff0c;看这里&#xff01;全程高清回播将为你持续解密这一次你一定不可以再错过了哦&#xff01;什么是Microsoft Learn&…

十分钟了解Kubernetes

何为Kubernetes?最简单的一句话来概括Kubernetes。它就是一套成熟的商用服务编排解决方案。Kubernetes定位在Paas层&#xff0c;重点解决了微服务大规模部署时的服务编排问题。Kubernetes组件介绍了解Kubernetes都是从Pod开始的。Pod是Kubernetes最小的调度单元&#xff0c;所…

【CF1194E】Count The Rectangles【类扫描线】【单调性】【树状数组】

传送门 题意&#xff1a;给定NNN条与坐标轴平行的线段&#xff0c;保证不垂直的线段没有交点&#xff0c;求一共构成多少个矩形&#xff08;以线段交点为顶点&#xff09;。 1≤N≤50001\leq N\leq50001≤N≤5000 显然是个数据结构乱搞题。 直觉告诉我们先枚举一条线段。 假…

.NET Core 的过去、现在和未来

在最新的微软 Build 大会主题网站上&#xff0c;微软.NET 程序经理 Scott Hunter 发表了一篇文章&#xff0c;指出.NET Core 是.NET 的未来。去年&#xff0c;微软就已经确认 Visual Studio Live 中的.NET 框架将会被.NET Core 取代。将.NET Core 作为开源开发技术栈在 2014 年…

Visual Basic 兴衰记

这是一篇暴露年龄的文章。为什么这么说&#xff1f;因为现在年轻的程序员可能没有接触过当年红极一时的 Visual Basic。28 年前的 1991 年 4 月&#xff0c;Microsoft 发布了 Visual Basic 1.0 for Windows&#xff0c;次年 9 月&#xff0c;发布了 Visual Basic 1.0 for DOS。…

函数式编程里的Materialization应该翻译成什么?

Materialization是函数式编程里的一个专业术语, 用于特指函数式编程中查询被实际执行并生成结果的这一过程.首先, 搜了一下中文资料, 暂时没有对该词的中文翻译, CSDN\博客园\阿里云上所有关于它的文档都没有做中文翻译, 直接沿用这个英文单词.难道轮到我来创造这个词的中文翻…