Codeforces Round #219 (Div. 1)

这场一开始A题用set乱贪一发就交了，然后WA第六组。。。。然后二分乱搞过的。B题真的是想晕了，4秒的时限就是骗人的嘛。。。最后用1600*1600复杂度过的。。。思维不行，还需要好好练习。。。

A：

题意：

      把一个袋鼠装进另一个袋鼠。每个袋鼠只能装一个袋鼠且必须体积小于等于自己的一半。不能装 装了袋鼠的袋鼠。求最少剩       多少只袋鼠。

思路：二分有多少只袋鼠被装。

code：

#include <algorithm>#include <iostream>#include <string.h>#include <stdlib.h>#include <stdio.h>#include <string>#include <math.h>#include <vector>#include <queue>#include <stack>#include <cmath>#include <list>#include <set>#include <map>using namespace std;#define N  500010#define ALL(x)     x.begin(),x.end()#define CLR(x,a)   memset(x,a,sizeof(x))typedef long long    ll;typedef pair<int,int> PI;const int INF    = 0x3fffffff;const int MOD    = 100000007;const double EPS = 1e-7;int n;int a[N];bool check(int x){for(int i=0;i<x;i++){if(2*a[i]>a[n-x+i]) return false;}return true;}int main(){scanf("%d",&n);for(int i=0;i<n;i++)scanf("%d",&a[i]);sort(a,a+n);int l=0,r=n/2+1;while(l<r){int mid=(l+r)>>1;if(check(mid)) l=mid+1;else r=mid;}printf("%d\n",n-r+1);    return 0;}

B：

题意：求矩阵里有多少个子矩阵全包含零。

思路：dp[x0][y0][x1][y1] 表示矩阵{x0<=x<=x1,y0<=y<=y1}有多少个包含1的矩阵，tot[x0][y0][x1][y1]则表示有多少个子矩阵，对减下就是答案，递推的方法看代码，类似容斥。

比赛的时候没有去递推tot，一直以为是算一下组合数就是了，然后狂卡样例。。。。真是弱。。代码略翔，递推太长了。

code：

#include <algorithm>#include <iostream>#include <string.h>#include <stdlib.h>#include <stdio.h>#include <string>#include <math.h>#include <vector>#include <queue>#include <stack>#include <cmath>#include <list>#include <set>#include <map>using namespace std;#define N  500010#define ALL(x)     x.begin(),x.end()#define CLR(x,a)   memset(x,a,sizeof(x))typedef long long    ll;typedef pair<int,int> PI;const int INF    = 0x3fffffff;const int MOD    = 100000007;const double EPS = 1e-7;int n,m;int sum[64][64][64][64];int tot[64][64][64][64];int main(){int x0,x1,y0,y1;scanf("%d%d",&n,&m);int Q;scanf("%d",&Q);for(int i=1;i<=n;i++)for(int j=1;j<=m;j++){scanf("%1d",&sum[i][j][i][j]);tot[i][j][i][j]=1;}for(int h=1;h<=n;h++){for(int w=1;w<=m;w++){if(h==1 && w==1) continue;for(int i=1;i<=n;i++){for(int j=1;j<=m;j++){int k=i+h-1,q=j+w-1;if(k>n || q>m) break;sum[i][j][k][q]+=sum[i+1][j][k][q]+sum[i][j+1][k][q]+sum[i][j][k-1][q]+sum[i][j][k][q-1];sum[i][j][k][q]-=sum[i+1][j+1][k][q]+sum[i+1][j][k-1][q]+sum[i+1][j][k][q-1]+sum[i][j+1][k-1][q]+sum[i][j+1][k][q-1]+sum[i][j][k-1][q-1];sum[i][j][k][q]+=sum[i+1][j+1][k-1][q]+sum[i+1][j+1][k][q-1]+sum[i+1][j][k-1][q-1]+sum[i][j+1][k-1][q-1];sum[i][j][k][q]-=sum[i+1][j+1][k-1][q-1];sum[i][j][k][q]+=(sum[i][j][k][q]>0);tot[i][j][k][q]+=tot[i+1][j][k][q]+tot[i][j+1][k][q]+tot[i][j][k-1][q]+tot[i][j][k][q-1];tot[i][j][k][q]-=tot[i+1][j+1][k][q]+tot[i+1][j][k-1][q]+tot[i+1][j][k][q-1]+tot[i][j+1][k-1][q]+tot[i][j+1][k][q-1]+tot[i][j][k-1][q-1];tot[i][j][k][q]+=tot[i+1][j+1][k-1][q]+tot[i+1][j+1][k][q-1]+tot[i+1][j][k-1][q-1]+tot[i][j+1][k-1][q-1];tot[i][j][k][q]-=tot[i+1][j+1][k-1][q-1];tot[i][j][k][q]++;}}}}while(Q--){scanf("%d%d%d%d",&x0,&y0,&x1,&y1);printf("%d\n",tot[x0][y0][x1][y1]-sum[x0][y0][x1][y1]);}    return 0;}

C:

题意：

        有n个看台排成一条直线。现在要放m次烟花。每次放烟花，你会获得b-|a-x|的收益，a表示烟花在哪个看台放，x是你所在的位置。每单位时间你最多能移动距离d，给出了每组烟花燃放的时间，求最大收益。

思路：

        dp[i][j] 表示燃放第i组烟花，此时人在第j个看台时候的最大收益。

        用dis表示从i-1到i最多能移动的距离，那么转移就是dp[i][j]=max{dp[i-1][k] | abs(j-k)<=dis }。用单调队列优化下，就能过了。

code：

#include <algorithm>#include <iostream>#include <string.h>#include <stdlib.h>#include <stdio.h>#include <string>#include <math.h>#include <vector>#include <queue>#include <stack>#include <cmath>#include <list>#include <set>#include <map>using namespace std;#define N  200010#define ALL(x)     x.begin(),x.end()#define CLR(x,a)   memset(x,a,sizeof(x))typedef long long    ll;typedef pair<int,int> PI;const int INF    = 0x3fffffff;const int MOD    = 100000007;const double EPS = 1e-7;deque<int> que;ll dp[2][N];int main(){int n,m,d,pret=0,pre=1,now=0;int a,b,t;scanf("%d%d%d",&n,&m,&d);for(int i=0;i<m;i++){scanf("%d%d%d",&a,&b,&t);ll dis=1ll*d*(t-pret);pret=t;que.clear();for(int j=1;j<=n;j++){while(!que.empty() && dp[pre][j]>dp[pre][que.back()]) que.pop_back();que.push_back(j);while(dis<abs(que.front()-j)) que.pop_front();dp[now][j]=dp[pre][que.front()];}que.clear();for(int j=n;j>=1;j--){while(!que.empty() && dp[pre][j]>dp[pre][que.back()]) que.pop_back();que.push_back(j);while(dis<abs(que.front()-j)) que.pop_front();dp[now][j]=max(dp[now][j],dp[pre][que.front()]);}for(int j=1;j<=n;j++) dp[now][j]+=b-abs(j-a);swap(pre,now);}printf("%I64d\n",*max_element(dp[pre]+1,dp[pre]+n+1));    return 0;}