bzoj4418&&bzoj4419&&bzoj4420:SHOI2013Day2题解

这三题截止现在（2016.3.11）窝居然都是跑的最快的……可啪……

T1 bzoj4418

这题叫做扇形面积并，看到这个名字我就方了，因为我不会计算几何啊QAQ

一看题目，发现是傻逼题……（雾）

又是一个树状数组上二分（辣鸡SHOI）

用BIT维护每个半径出现的次数，然后二分找到第k大的，如果不足k个就是0，然后扫一遍就完了

另外我没看懂m<=10^6的意义……因为这个条件完全用不到……

另外这题我手里的考场版本（众所周知窝是大弱省SH的选手）是说a1<a2……然后数据又有a1>a2……辣鸡出题人颓我时间！辣鸡SHOI吃枣药丸！（大雾

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#define ll long long
#define N 400005
 
using namespace std;
inline int read(){
	int ret=0;char ch=getchar();
	bool flag=0;
	while (ch<'0'||ch>'9'){
		flag=ch=='-';
		ch=getchar();
	}
	while ('0'<=ch&&ch<='9'){
		ret=ret*10-48+ch;
		ch=getchar();
	}
	return flag?-ret:ret;
}
 
struct BIT{
	int c[N];
	void clear(){memset(c,0,sizeof(c));}
	int lowbit(int x){return x&-x;}
	void modify(int pos,int delta){
		for (int i=pos;i<=100000;i+=lowbit(i))c[i]+=delta;
	}
	int solve(int value){
		int x=0,now=0;
		for (int i=(1<<16);i;i=(i>>1))
			if (x+i<=100000&&now+c[x+i]<=value)
				now+=c[x+=i];
		return x;
	}
} bit;
 
struct node{
	int pos,r,claris;
	node(){}
	node(int _pos,int _r,int _claris):pos(_pos),r(_r),claris(_claris){}
} a[N];
inline bool operator <(const node &x,const node &y){
	return x.pos<y.pos;
}
int n,K;
 
int main(){
	int n0=read(),m0=read();K=-read();n=0;
	while (n0--){
		int x=read(),a1=read(),a2=read();
		a[++n]=node(a1,x,1);
		a[++n]=node(a2,x,-1);
		if (a1>a2){
			a[++n]=node(-m0,x,1);
			a[++n]=node(m0,x,-1);
		}
	}
	sort(a+1,a+n+1);
	bit.clear();
	ll ans=0;
	for (int i=1;i<n;++i){
		K+=a[i].claris;
		bit.modify(a[i].r,a[i].claris);
		if (a[i].pos==a[i+1].pos||K<0) continue;
		ll tmp=bit.solve(K)+1;
		ans+=tmp*tmp*(a[i+1].pos-a[i].pos);
	}
	printf("%lld\n",ans);
	return 0;
}

T2 bzoj4419

这题叫做发微博……有普及组水平的肯定能秒了这个题啦！

我们计算每次加/删好友对答案的贡献，我们倒过来做（顺过来肯定可以，我懒得想了而已），用一个数组维护每个人发的微博数。

对于每一条加好友信息u+v，对u的答案有当前v发过的微博数的贡献，反之亦然。

对于删好友，取负即可。

由于每次加/删好友都是合法的，所以这个算法的正确性是显然的。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#define N 500005
 
using namespace std;
inline int read(){
	int ret=0;char ch=getchar();
	while (ch<'0'||ch>'9') ch=getchar();
	while ('0'<=ch&&ch<='9'){
		ret=ret*10-48+ch;
		ch=getchar();
	}
	return ret;
}
 
int op[N],u[N],v[N];
int sum[N],ans[N];
int n,m;
 
int main(){
	n=read();m=read();
	for (int i=1;i<=m;u[i++]=read()){
		char ch=getchar();
		while (ch!='!'&&ch!='-'&&ch!='+')ch=getchar();
		op[i]=ch!='-'?ch=='+':-1;
		if (op[i]) v[i]=read();
	}
	memset(sum,0,sizeof(sum));
	memset(ans,0,sizeof(ans));
	for (int i=m;i;--i)
		if (op[i]){
			ans[u[i]]+=op[i]*sum[v[i]];
			ans[v[i]]+=op[i]*sum[u[i]];
		}
		else sum[u[i]]++;
	for (int i=1;i<n;++i) printf("%d ",ans[i]);
	printf("%d\n",ans[n]);
	return 0;
}

T3 bzoj4420

这题叫做二重镇……看完题我就想吐槽，这不就是个消除游戏吗！辣鸡出题人尽我食粮！（雾

这题看规模就知道是状压DP（咦…怎么又考状压…辣鸡SHOI吃枣药丸！）

然后听说这题是码农题，我觉得顶多算个普及组级别的码农题（？）NOIP2015的斗地主都比这个题难写啊……

先dfs找出所有状态，以及各个状态加各物品能延伸到的状态和得分

这里有个小优化，把所有能消除的状态全部剪掉，大概能剪掉2/5吧

其实这里还有一个小优化，不过我没有做（反正测速已经rank1了？），就是把对称的状态全部剪掉，这样又能剪掉将近一半

然后，我们用dp[i][k][mask]表示已经处理了前i个物品，仓库里是物品k，当前状态是mask的最大得分，很显然能脑补转移，然后这题就做完了

（不能脑补转移的，看我代码里的int DP(){}里有很详细的转移姿势）

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#define N 8
#define L 105
#define MASK 33333
 
using namespace std;
inline void Max(int &x,int y){if (y>x) x=y;}
 
int bin[N];
int id[666669],cnt;
int n,len;
int getcode(int *a,int lastpos,int &point){
	point=0;
	for (;a[lastpos];){
		int l=lastpos,r=lastpos,tmp=a[lastpos];
		while (l>1&&a[l]==a[l-1]) --l;
		while (r<n&&a[r]==a[r+1]) ++r;
		if (l==r) break;
		point+=(r-l+1)*bin[tmp];
		for (int i=l;i<=r;++i) a[i]=0;
		a[lastpos]=(tmp+1)%6;
	}
	int tmp=0;
	for (int i=1;i<=n;++i) (tmp*=10)+=a[i];
	if (!id[tmp]) id[tmp]=cnt++;
	return id[tmp];
}
 
int a[N],b[N];
int dgr[MASK],edge[MASK][N][N],pnt[MASK][N][N];
void dfs(int u){
	if (u>n){
		int now=getcode(a,0,b[0]);
		for (int i=1;i<=n;++i)if (!a[i]){
			++dgr[now];
			for (int j=1;j<=5;++j){
				for (int k=1;k<=n;++k) b[k]=a[k];
				b[i]=j;
				edge[now][j][dgr[now]]=getcode(b,i,pnt[now][j][dgr[now]]);
			}
		}
		return;
	}
	for (int i=0;i<=5;++i)
		if (u==1||!a[u-1]||a[u-1]!=i)
			a[u]=i,dfs(u+1);
}
 
int dp[L][N][MASK],gd[L];
int DP(){
	memset(dp,-1,sizeof(dp));
	int ans=dp[0][0][0]=0;
	for (int i=0;i<=len;++i)
		for (int k=5;k>=0;--k)
			for (int mask=0;mask<cnt;++mask)if (dp[i][k][mask]>=0){
				if (i<len) for (int j=1;j<=dgr[mask];++j)
					Max(dp[i+1][k][edge[mask][gd[i+1]][j]],dp[i][k][mask]+pnt[mask][gd[i+1]][j]);
				if (k) for (int j=1;j<=dgr[mask];++j)
					Max(dp[i][0][edge[mask][k][j]],dp[i][k][mask]+pnt[mask][k][j]);
				else Max(dp[i+1][gd[i+1]][mask],dp[i][k][mask]);
				Max(ans,dp[i][k][mask]);
			}
	return ans;
}
 
int main(){
	scanf("%d%d",&n,&len);
	for (int i=1;i<=len;++i){
		char ch=getchar();
		while (ch<'0'||ch>'9') ch=getchar();
		gd[i]=ch-48;
	}
	for (int i=bin[0]=1;i<=5;++i) bin[i]=bin[i-1]<<1;
	memset(dgr,0,sizeof(dgr));cnt=0;
	dfs(1);
	printf("%d\n",DP());
	return 0;
}

ps：听说这场lxlwdgy大爷AK了%%%

pss：这篇文章的标签我直接复制上一篇的再删一个就行了……