码迷,mamicode.com
首页 > 其他好文 > 详细

bzoj2436: [Noi2011]Noi嘉年华

时间:2019-01-17 22:35:08      阅读:141      评论:0      收藏:0      [点我收藏+]

标签:unique   高仿   sort   mes   ++   ||   font   ora   ini   

我震惊了,我好菜,我是不是该退役(苦逼)

 

可以先看看代码里的注释

首先我们先考虑一下第一问好了真做起来也就这个能想想了

那么离散化时间是肯定的,看一手范围猜出是二维DP,那对于两个会场,一个放自变量,一个放变量,然后O(n^3)的DP好了

第二问像第一问的做法特判一波就是O(n^4)啦

对于一个嘉年华必选,等价于必选一段区间,我们设f[l][r]为必选l,r放一起,前面一段自己处理,后面一段自己处理的最优解

那么ans=max(f[l][r]) (a[i].l<=l,r<=a[i].r)

可以发现前面一段自己处理我们在第一问已经搞定了。。。后面高仿前面就好。。。

f[l][r]=max(min(s[l-1][x]+h[l][r]+t[r+1][y],x+y)) 

然后这还是个四方的

但是用脑(bai)子(du)想想,t[r+1][y]随y增大减小,把min中两项写成两个一次函数,这个min的图像会是一个v字形(一个下降的直线和一个上升的直线),v的最下方就是决策点

当x增大的时候由于s[l-1][x]也跟着减小,相当于前一个截距变小,而后一个截距变大,决策点左移,所以是有决策单调性的y可以扫一遍完事

 

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn=210;
const int maxT=2*maxn;

struct node{int l,r;}a[maxn];
int lslen,ls[maxT];

int h[maxT][maxT];//时间段内------->有多少嘉年华 
int s[maxT][maxn],t[maxT][maxn];//前/后到i的时间段,给第一个j,另一个最多混到多少 
int f[maxT][maxT];//这个时间段必选,且没有选择和该时间段相交的嘉年华------>较小的最大为多少 
int main()
{
    freopen("a.in","r",stdin);
    freopen("a.out","w",stdout);
    int n;
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d%d",&a[i].l,&a[i].r);a[i].r+=a[i].l;
        a[i].l++;
        ls[++lslen]=a[i].l;
        ls[++lslen]=a[i].r;
    }
    sort(ls+1,ls+lslen+1);
    lslen=unique(ls+1,ls+lslen+1)-ls-1;
    for(int i=1;i<=n;i++)
        a[i].l=lower_bound(ls+1,ls+lslen+1,a[i].l)-ls,
        a[i].r=lower_bound(ls+1,ls+lslen+1,a[i].r)-ls;
        
    for(int l=1;l<=lslen;l++)
        for(int r=l;r<=lslen;r++)
            for(int i=1;i<=n;i++)
                if(l<=a[i].l&&a[i].r<=r)h[l][r]++;
                
    //.......init........
    
    memset(s,-31,sizeof(s));s[0][0]=0;
    for(int i=0;i<lslen;i++)
        for(int j=0;j<=n;j++)
            if(s[i][j]!=-1)
                for(int k=i+1;k<=lslen;k++)
                    s[k][j]=max(s[k][j],s[i][j]+h[i+1][k]),
                    s[k][j+h[i+1][k]]=max(s[k][j+h[i+1][k]],s[i][j]);
    int mx=0;
    for(int j=0;j<=n;j++)
        mx=max(mx,min(j,s[lslen][j]));
    printf("%d\n",mx);

    memset(t,-31,sizeof(t));t[lslen+1][0]=0;
    for(int i=lslen+1;i>1;i--)
        for(int j=0;j<=n;j++)
            if(t[i][j]!=-1)
                for(int k=i-1;k>=1;k--)
                    t[k][j]=max(t[k][j],t[i][j]+h[k][i-1]),
                    t[k][j+h[k][i-1]]=max(t[k][j+h[k][i-1]],t[i][j]);
        
    //......solve1.......
    
    for(int l=1;l<=lslen;l++)
        for(int r=l;r<=lslen;r++)
        {
            int y=n;
            for(int x=0;x<=n;x++)
            {
                if((x+y)<=n)
                    f[l][r]=max(f[l][r],min(s[l-1][x]+h[l][r]+t[r+1][y],x+y));
                while(y>=0&&((x+y)>n||s[l-1][x]+h[l][r]+t[r+1][y]<x+y))
                {
                    y--;
                    f[l][r]=max(f[l][r],min(s[l-1][x]+h[l][r]+t[r+1][y],x+y));
                }
            }
        }
    
    for(int i=1;i<=n;i++)
    {
        int ans=0;
        for(int l=1;l<=a[i].l;l++)
            for(int r=a[i].r;r<=lslen;r++)
                ans=max(ans,f[l][r]);
        printf("%d\n",ans);
    }
    
    return 0;
}

 

bzoj2436: [Noi2011]Noi嘉年华

标签:unique   高仿   sort   mes   ++   ||   font   ora   ini   

原文地址:https://www.cnblogs.com/AKCqhzdy/p/10284820.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!