BZOJ4426 : [Nwerc2015]Better Productivity最大生产率
2016-03-18 17:14
459 查看
如果一个区间包含另一个区间,那么这两个区间是否在一起的生产率是一样的。
将所有这种包含了其他区间的区间放入数组$b$,其余的放入数组$c$,有多个相同的时候则从$b$移一个到$c$。
那么$c$里所有区间左端点递增,右端点也递增,设$f[i][j]$为$c$中前$j$个区间划分成$i$组的最大收益,直接DP即可,决策具有单调性。
然后把$p$分配给$b$和$c$,求出$b$和$c$组合取来的最大收益即可。
时间复杂度$O(n^2\log n)$。
将所有这种包含了其他区间的区间放入数组$b$,其余的放入数组$c$,有多个相同的时候则从$b$移一个到$c$。
那么$c$里所有区间左端点递增,右端点也递增,设$f[i][j]$为$c$中前$j$个区间划分成$i$组的最大收益,直接DP即可,决策具有单调性。
然后把$p$分配给$b$和$c$,求出$b$和$c$组合取来的最大收益即可。
时间复杂度$O(n^2\log n)$。
#include<cstdio> #include<algorithm> #define N 210 using namespace std; int n,m,i,j,flag,cb,cc,s ,f ,ans=-2147400000; struct P{int x,y;}a ,b ,c ; bool cmpb(P a,P b){return a.y-a.x>b.y-b.x;} bool cmpc(P a,P b){return a.x<b.x;} void dp(int p,int l,int r,int dl,int dr){ int m=(l+r)>>1,dm=dl,t=ans; for(int i=min(m-1,dr);i>=dl;i--){ if(c[i+1].y<=c[m].x)break; int now=f[p-1][i]+c[i+1].y-c[m].x; if(now>=t)t=now,dm=i; } f[p][m]=t; if(l<m)dp(p,l,m-1,dl,dm); if(r>m)dp(p,m+1,r,dm,dr); } int main(){ scanf("%d%d",&n,&m); for(i=1;i<=n;i++)scanf("%d%d",&a[i].x,&a[i].y); for(i=1;i<=n;i++){ for(flag=0,j=1;j<=n;j++)if(a[i].x<=a[j].x&&a[j].y<=a[i].y&&(a[i].x!=a[j].x||a[i].y!=a[j].y||i<j)){flag=1;break;} if(flag)b[++cb]=a[i];else c[++cc]=a[i]; } sort(b+1,b+cb+1,cmpb); for(i=1;i<=cb;i++)s[i]=s[i-1]+b[i].y-b[i].x; sort(c+1,c+cc+1,cmpc); for(i=1;i<=cc;i++)f[0][i]=ans; for(i=1;i<=m;i++)f[i][0]=ans; for(i=1;i<=m;i++)dp(i,1,cc,0,cc); for(i=1;i<=m;i++)if(m-i<=cb&&f[i][cc]>=0)ans=max(ans,f[i][cc]+s[m-i]); return printf("%d",ans),0; }
相关文章推荐
- jquery.mobile新手说明与示例
- 详解spl_autoload_register()函数
- 在windows xp x86中配置python开发环境
- error C4235: 使用了非标准扩展: 不支持在此结构上使用“_asm”关键字
- python3.5之string
- @RequestMapping("/myway")
- iOS中.pch文件如何使用
- 程序员也需要这种魄力
- MyBatis对于Java对象里的枚举类型处理
- SpringMvc+Mybatis常用插件
- super()函数
- 使用Genymotion调试出现错误INSTALL_FAILED_CPU_ABI_INCOMPATIBLE解决办法
- TLCL学习笔记
- no debuggable application
- mysql中二进制形式保存文件数据
- 蓝桥杯模拟试题—找数字
- Android 打造底部导航控件
- 光源,摄像机,简单的工厂
- c++与c不太相同的一些地方2
- 备份mysql显示mysqldump: Got errno 28 on write