标签:
http://www.lydsy.com/JudgeOnline/problem.php?id=2436
首先离散化,离散化后时间范围为[1,cnt]。
求出H[i][j],表示时间范围在[i,j]的活动有多少个,可以在N^2的时间内解决。
假设场地分别为A和B。
我们容易知道,场地A和场地B的活动安排一定是这样的:
他们的活动安排一定是这样间隔着的。
我们求F[i][j],表示当时间<=i时,A场地有j个活动,B场地最多有多少个活动。
我们这样记并不知道最后一个活动是在A场地还是B场地,但一定是这2种情况的最大值,也就是说:F[i][j]表示“当时间<=i时,A场地有j个活动且最后一个活动在A场地,B场地最多有多少个活动”和当时间<=i时,A场地有j个活动且最后一个活动在B场地,B场地最多有多少个活动"这两种情况的最大值。
从小到大枚举i。
先要倒序循环j,用F[i][j+1]更新F[i][j]:F[i][j]=max(F[i][j],F[i][j+1])。
我们枚举j,向后递推,枚举k,
如果接下来的[i+1,k]的时间区间内都是B场地在进行活动,我们可以转移到F[k][j]:F[k][j]=max(F[k][j],F[i][j]+H[i+1][k])。
如果接下来的[i+1,k]的时间区间内都是A场地在进行活动,我们可以转移到F[k][j+H[i+1][k]]:F[k][j+H[i+1][k]]=max(F[k][j+H[i+1][k]],F[i][j])。
这个可以在N^3的时间内解决。
类似的,我们求G[i][j],表示当时间>=i时,A场地有j个活动,B场地最多有多少个活动,也是N^3的时间内解决。
然后我们求T[i][j],表示[i,j]内的所有活动都在A场地时,活动相对较少的嘉年华的活动数量的最大值。
很直接的一个想法是:
$T[i][j]=max\{min(x+y+H[i][j],F[i-1][x]+G[j+1][y])\} (0\leqslant x\leqslant N,0\leqslant y\leqslant N)$
但是这样算肯定是N^4的。
我们发现,当x变大的时候,F[i-1][x]跟着变小:
$min(x↑+y+H[i][j],F[i-1][x]↓+G[j+1][y])$
如果y也跟着变大,G[j+1][y]跟着变小:
$min(x↑+y↑+H[i][j],F[i-1][x]↓+G[j+1][y]↓)$
由于我们取的是min,这样并没有什么卵用。
所以y只能变小,G[j+1][y]跟着变大:
$min(x↑+y↓+H[i][j],F[i-1][x]↓+G[j+1][y]↑)$
所以随着x的递增,y递减。
并且容易知道,这一定是个单峰的。
所以可以N^3求出ans[i][j]。
第1行的输出就是$max\{min(i,F[cnt][i])\}(0\leqslant i\leqslant N)$
第i+1个输出就是$max(T[i][j])(1\leqslant i\leqslant l,r\leqslant j\leqslant cnt)$
#include<cstdio> #include<cstdlib> #include<iostream> #include<fstream> #include<algorithm> #include<cstring> #include<string> #include<cmath> #include<queue> #include<stack> #include<map> #include<utility> #include<set> #include<bitset> #include<vector> #include<functional> #include<deque> #include<cctype> #include<climits> #include<complex> //#include<bits/stdc++.h>适用于CF,UOJ,但不适用于poj using namespace std; typedef long long LL; typedef double DB; typedef pair<int,int> PII; typedef complex<DB> CP; #define mmst(a,v) memset(a,v,sizeof(a)) #define mmcy(a,b) memcpy(a,b,sizeof(a)) #define re(i,a,b) for(i=a;i<=b;i++) #define red(i,a,b) for(i=a;i>=b;i--) #define fi first #define se second #define m_p(a,b) make_pair(a,b) #define SF scanf #define PF printf #define two(k) (1<<(k)) template<class T>inline T sqr(T x){return x*x;} template<class T>inline void upmin(T &t,T tmp){if(t>tmp)t=tmp;} template<class T>inline void upmax(T &t,T tmp){if(t<tmp)t=tmp;} const DB EPS=1e-9; inline int sgn(DB x){if(abs(x)<EPS)return 0;return(x>0)?1:-1;} const DB Pi=acos(-1.0); inline int gint() { int res=0;bool neg=0;char z; for(z=getchar();z!=EOF && z!=‘-‘ && !isdigit(z);z=getchar()); if(z==EOF)return 0; if(z==‘-‘){neg=1;z=getchar();} for(;z!=EOF && isdigit(z);res=res*10+z-‘0‘,z=getchar()); return (neg)?-res:res; } inline LL gll() { LL res=0;bool neg=0;char z; for(z=getchar();z!=EOF && z!=‘-‘ && !isdigit(z);z=getchar()); if(z==EOF)return 0; if(z==‘-‘){neg=1;z=getchar();} for(;z!=EOF && isdigit(z);res=res*10+z-‘0‘,z=getchar()); return (neg)?-res:res; } const int maxN=200; const int maxcnt=2*maxN; int N; struct Tdata{int l,r,id;}a[maxN+10]; int cnt,bak[2*maxN+100]; int t[maxcnt+10]; int H[maxcnt+10][maxcnt+10]; inline bool cmpr(Tdata x,Tdata y){return x.r<y.r;} int F[maxcnt+10][maxN+10],G[maxcnt+10][maxN+10]; int T[maxcnt+10][maxcnt+10]; int ans,res[maxN+10]; int main() { /*freopen("show.in","r",stdin); freopen("show.out","w",stdout);*/ int i,j,k; N=gint(); re(i,1,N)a[i].l=gint(),a[i].r=a[i].l+gint()-1,a[i].id=i; re(i,1,N)bak[++cnt]=a[i].l,bak[++cnt]=a[i].r; sort(bak+1,bak+cnt+1); cnt=unique(bak+1,bak+cnt+1)-bak-1; re(i,1,N)a[i].l=lower_bound(bak+1,bak+cnt+1,a[i].l)-bak,a[i].r=lower_bound(bak+1,bak+cnt+1,a[i].r)-bak; sort(a+1,a+N+1,cmpr); int head=1; re(i,1,cnt) { while(head<=N && a[head].r<=i)t[a[head].l]++,head++; int sum=0; red(j,i,1)sum+=t[j],H[j][i]=sum; } mmst(F,-1); F[0][0]=0; re(i,0,cnt-1) { red(j,N-1,0)upmax(F[i][j],F[i][j+1]); re(j,0,N) { if(F[i][j]!=-1)re(k,i+1,cnt)upmax(F[k][j],F[i][j]+H[i+1][k]); if(F[i][j]!=-1)re(k,i+1,cnt)upmax(F[k][j+H[i+1][k]],F[i][j]); } } mmst(G,-1); G[cnt+1][0]=0; red(i,cnt+1,2) { red(j,N-1,0)upmax(G[i][j],G[i][j+1]); re(j,0,N) { if(G[i][j]!=-1)re(k,1,i-1)upmax(G[k][j],G[i][j]+H[k][i-1]); if(G[i][j]!=-1)re(k,1,i-1)upmax(G[k][j+H[k][i-1]],G[i][j]); } } ans=0; re(j,0,N)upmax(ans,min(j,F[cnt][j])); cout<<ans<<endl; re(i,1,cnt)re(j,i,cnt) { T[i][j]=0; /*int x,y; re(x,0,N)re(y,0,N) if(F[i-1][x]!=-1 && G[j+1][y]!=-1) upmax(T[i][j],min(x+y+H[i][j],F[i-1][x]+G[j+1][y]));*/ int x,y=N; while(y-1>=0 && G[j+1][y]==-1)y--; re(x,0,N) { if(F[i-1][x]==-1)continue; while(y-1>=0 && min(x+y+H[i][j],F[i-1][x]+G[j+1][y])<=min(x+y-1+H[i][j],F[i-1][x]+G[j+1][y-1]))y--; upmax(T[i][j],min(x+y+H[i][j],F[i-1][x]+G[j+1][y])); } } re(i,1,N) { int l=a[i].l,r=a[i].r,id=a[i].id; res[id]=0; re(j,1,l)re(k,r,cnt)upmax(res[id],T[j][k]); } re(i,1,N)cout<<res[i]<<endl; return 0; }
标签:
原文地址:http://www.cnblogs.com/maijing/p/4700986.html