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;
}
题意：求矩阵里有多少个子矩阵全包含零。
思路：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;
}