[游记&题解]2019暑假中山纪中集训day11&12

最新推荐文章于 2020-06-22 12:42:56 发布

C20190406Panda_hu

最新推荐文章于 2020-06-22 12:42:56 发布

阅读量185

点赞数 1

分类专栏：游记合辑

本文链接：https://blog.csdn.net/qq_34940287/article/details/102934488

版权

游记合辑专栏收录该内容

14 篇文章 0 订阅

订阅专栏

中山纪中NOIP2019提高组模拟赛

A组 day12 solution

Date:2019.8.12

Probelm	A	B	C
题目名称	迷宫	猛汉王	工厂
Difficulty	2300	2600	2700
Finished	Yes	-	-

简单游记

day 11

忘了写了…

day 12

AM
今天是自闭的第8天。
打开先看A题…看起来还挺不好做的…
然后发现很像之前做过的那道跳房子的题，像那道题一样直接维护一个5$\times $5 的 d p 数组就行了 . . . 然后就开始写，中间转移的部分由于脑残写复杂了 . . . 于是调了半天 . . . 然后开专门出了两组比较强的数据，这时候已经过去两个小时了 . . . 然后看 B . . . . 看起来好像只有暴力分 ? ? ? 然后没看多久就调了。然后看 C . . . 好像可以用状压来做，然后试着写了一个$ O(2^nn)$写完之后发现样例都过不了…试着加了一点贪心…还是过不了…
只好试着打更暴力的分数…
结果没有打完，考试就差不多结束了…
PM
回来看成绩(只交了A题)，发现A题又写挂了…
后来检查的时候发现少写了一行…
不过自己写超过100行的代码多半都要写炸…这也不是第一次了…
B题是又是扫描线+线段树???
C题是一个很神仙的状压dp…

迷宫 (maze)

[题解]

算法标签:线段树。

这年头什么东西都能在线段树上维护了。

在线段树上维护dp数组即可。

[实现]

#pragma GCC optimize(2)
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<cstring>
using namespace std;
#define MAXN 200000
#define Pr pair<int,int>
#define X first
#define Y second
#define INF 1000000000
#define RI register int
int read(){
	int x=0,F=1;char c=getchar();
	while(c<'0'||c>'9'){if(c=='-')F=-1;c=getchar();}
	while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
	return x*F;
}
int n,m,q;
int min(int a,int b){return a<b?a:b;}
struct node{
	int dp[5][5],l[5],r[5];
	node(){
		for(RI i=0;i<5;i++)
			for(RI j=0;j<5;j++)
			dp[i][j]=INF;
	}
	node operator +(const node &b)const
	{	node tmp;
		for(RI i=0;i<n;i++)
			for(RI j=0;j<n;j++)
				for(RI k=0;k<n;k++)
				if(r[k]&&b.l[k])tmp.dp[i][j]=min(tmp.dp[i][j],dp[i][k]+b.dp[k][j]);
		for(RI i=0;i<n;i++)tmp.l[i]=l[i];
		for(RI i=0;i<n;i++)tmp.r[i]=b.r[i];
		return tmp;
	}
}p[(MAXN<<2)+5],pre[MAXN+5];
inline void Build(int x,int l,int r){
	if(l==r){
		p[x]=pre[l];
		return ;
	}
	int mid=(l+r)>>1;
	Build(x<<1,l,mid);
	Build(x<<1|1,mid+1,r);
	p[x]=p[x<<1]+p[x<<1|1];
}
inline void modf(int x,int l,int r,int pos){
	if(l==r){
		p[x]=pre[pos];
		return ;
	}
	int mid=(l+r)>>1;
	if(pos<=mid)modf(x<<1,l,mid,pos);
	else modf(x<<1|1,mid+1,r,pos);
	p[x]=p[x<<1]+p[x<<1|1];
}
node Query(int x,int l,int r,int L,int R){
	if(L<=l&&R>=r)return p[x];
	int mid=(l+r)>>1;
	if(R<=mid)return Query(x<<1,l,mid,L,R);
	else if(L>mid)return Query(x<<1|1,mid+1,r,L,R);
	else
	return Query(x<<1,l,mid,L,mid)+Query(x<<1|1,mid+1,r,mid+1,R);
}
void prepare(int id){
	for(RI i=0;i<n;i++)
	{
		if(!pre[id].l[i])continue;
		for(RI j=i-1;j>=0;j--)
		if(pre[id].l[j])pre[id].dp[i][j]=i-j;
		else break;
		for(RI j=i+1;j<n;j++)
		if(pre[id].l[j])pre[id].dp[i][j]=j-i;
		else break;
	}
}
int main()
{
	freopen("maze.in","r",stdin);
	freopen("maze.out","w",stdout);
	n=read(),m=read(),q=read();
	for(RI i=0;i<n;i++)
		for(RI j=1;j<=m;j++)
		{
			RI val=read();
			if(val)pre[j].dp[i][i]=0;
			pre[j].l[i]=pre[j].r[i]=val;
		}
	for(RI j=1;j<=m;j++)
	prepare(j);
	Build(1,1,m);
	for(int i=1;i<=q;i++){
		RI opt=read();
		if(opt==1){
			RI a=read()-1,b=read();
			pre[b].l[a]=pre[b].r[a]=pre[b].l[a]^1;
			for(RI i=0;i<n;i++)
				for(int j=0;j<n;j++)
				pre[b].dp[i][j]=INF;
			for(RI i=0;i<n;i++)
				if(pre[b].l[i])
				pre[b].dp[i][i]=0;
			prepare(b);
			modf(1,1,m,b);
		}else{
			RI a=read()-1,b=read(),c=read()-1,d=read();
			if(d<b){
				printf("-1\n");
				continue;
			}
			node tp=Query(1,1,m,b,d);
			if(tp.dp[a][c]!=INF)printf("%d\n",tp.dp[a][c]+d-b);
			else printf("-1\n");
		}
	}
}

猛汉王 (mhw)

[题解]

算法标签:线段树,扫描线,数学推导
GuGu

[实现]

#pragma GCC optimize(3)
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<cstring>
using namespace std;
#define MAXN 100000
#define MAXM 600000
#define LL long long
int n,m,D,pcnt;
LL mx,mn;
int read(){
	int x=0,F=1;char c=getchar();
	while(c<'0'||c>'9'){if(c=='-')F=-1;c=getchar();}
	while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
	return x*F;
}
struct Line{
	LL x,l,r;
	int f,kdf,val,id;
}p[MAXM+5];
int lcnt;
LL l[MAXM+5],cov[2][MAXN+5];
bool cmp(Line s1,Line s2){
	if(s1.x==s2.x)return s1.val>s2.val;
	return s1.x<s2.x;
}
struct Seg_tree{
	int sum[(MAXM<<2)+5];
	Seg_tree(){memset(sum,0,sizeof(sum));}
	void Insert(int x,int l,int r,int L,int R,int val){
		if(r<L||R<l)return ;
		if(l>=L&&R>=r){
			sum[x]+=val;return ;
		}
		int mid=(l+r)>>1;
		Insert(x<<1,l,mid,L,R,val);
		Insert(x<<1|1,mid+1,r,L,R,val);
	}
	int Query(int x,int l,int r,int pos){
		if(l==r)return sum[x];
		int mid=(l+r)>>1;
		if(pos<=mid)return Query(x<<1,l,mid,pos)+sum[x];
		else return Query(x<<1|1,mid+1,r,pos)+sum[x];
	}
}T[2];
int main()
{
	freopen("mhw.in","r",stdin);
	freopen("mhw.out","w",stdout);
	n=read(),m=read(),D=read();
	for(int i=1;i<=n;i++){
		LL a=read(),b=read(),x,y;
		x=a+b,y=a-b;
		p[++lcnt]=(Line){x-D,y-D,y+D,0,0,1,i};
		p[++lcnt]=(Line){x,y,0,0,1,0,i};
		p[++lcnt]=(Line){x+D,y-D,y+D,0,0,-1,i};
		l[++l[0]]=y-D,l[++l[0]]=y,l[++l[0]]=y+D;
	}
	for(int i=1;i<=m;i++){
		LL a=read(),b=read(),x,y;
		x=a+b,y=a-b;
		p[++lcnt]=(Line){x-D,y-D,y+D,1,0,1,i};
		p[++lcnt]=(Line){x,y,0,1,1,0,i};
		p[++lcnt]=(Line){x+D,y-D,y+D,1,0,-1,i};
		l[++l[0]]=y-D,l[++l[0]]=y,l[++l[0]]=y+D;
	}
	sort(p+1,p+lcnt+1,cmp);
	sort(l+1,l+l[0]+1);
	l[0]=unique(l+1,l+l[0]+1)-l-1;
	for(int i=1;i<=lcnt;i++){
		p[i].l=lower_bound(l+1,l+l[0]+1,p[i].l)-l;
		p[i].r=lower_bound(l+1,l+l[0]+1,p[i].r)-l;
		if(p[i].kdf){
			cov[p[i].f][p[i].id]=T[p[i].f^1].Query(1,1,l[0],p[i].l);
		}else{
			T[p[i].f].Insert(1,1,l[0],p[i].l,p[i].r,p[i].val);
		}
	}
	//为了方便统计 cov[1][i]表示白点与黑点i距离大于d的个数
	// cov[0][i]表示黑点与白点i距离小于d的个数
	for(int i=1;i<=m;i++)
	cov[1][i]=n-cov[1][i];
	sort(cov[0]+1,cov[0]+n+1);
	sort(cov[1]+1,cov[1]+m+1);
	for(int i=1;i<=n;i++){
		mn+=1LL*(n-i)*cov[0][i];//贪心算每个cover的贡献
		mx+=1LL*(i-1)*cov[0][i];
		LL val=1LL*(m-cov[0][i])*(m-cov[0][i]-1)/2;//根据公式,这里要取反
		mn-=val,mx-=val;
	}
	for(int i=1;i<=m;i++){
		mn+=1LL*(m-i)*cov[1][i];
		mx+=1LL*(i-1)*cov[1][i];
		LL val=1LL*(n-cov[1][i])*(n-cov[1][i]-1)/2;
		mn-=val,mx-=val;
	}
	printf("%lld %lld",mn,mx);
}

工厂 (factory)

[题解]

算法标签:状压dp,二分图

GuGu

[实现]

#pragma GCC optimize(2)
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<cstring>
using namespace std;
#define MAXN 30
int n,decp,INF;
int read(){
	int x=0,F=1;char c=getchar();
	while(c<'0'||c>'9'){if(c=='-')F=-1;c=getchar();}
	while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
	return x*F;
}
int G[MAXN+1][MAXN+1],vis[MAXN*2+5];
int cnt[MAXN+1][MAXN+1],tot;
int f[300000][MAXN+1];
char s[MAXN+1];
struct node{
	int cnt,a,b,base;
}a[MAXN+1];
void xfind(int x,int sd,int &l,int &r){
	if(sd==0)l++;else r++;
	vis[x+sd*n]=1;
	for(int j=1;j<=n;j++){
		if(sd==0&&G[x][j]&&!vis[j+n])xfind(j,sd^1,l,r);
		if(sd==1&&G[j][x]&&!vis[j])xfind(j,sd^1,l,r);
	}
}
bool sum(int S){
	int x=0,y=0;
	for(int i=1;i<=tot;i++){
		int tmp=S%(a[i].cnt+1);
		x+=a[i].a*tmp,y+=a[i].b*tmp;
		S/=(a[i].cnt+1);
	}if(x==y)return 1;
	else return 0;
}
int main()
{
	freopen("factory.in","r",stdin);
	freopen("factory.out","w",stdout);
	n=read();
	for(int i=1;i<=n;i++){
		scanf("%s",s+1);
		for(int j=1;j<=n;j++){
			G[i][j]=s[j]-'0';
			decp+=s[j]-'0';//记录已经存在的边 
		}
	}
	for(int i=1;i<=2*n;i++)
	if(!vis[i]){
		int l=0,r=0;
		xfind((i-1)%n+1,(i-1)/n,l,r);
		cnt[l][r]++;
	}
	for(int i=0;i<=n;i++)
		for(int j=0;j<=n;j++)
		if(cnt[i][j]>0)
		a[++tot]=(node){cnt[i][j],i,j};
	int ms=0,tmp=1;
	for(int i=1;i<=tot;i++){
		ms+=a[i].cnt*tmp;a[i].base=tmp;
		tmp*=(a[i].cnt+1);
	}
	a[tot+1].base=tmp;
	memset(f,127,sizeof(f));
	INF=f[0][0];f[0][0]=0;
	for(int S=0;S<=ms;S++){
		int F=sum(S);
		for(int i=n;i>=0;i--){
			if(f[S][i]==INF)continue;
			for(int j=1;j<=tot;j++){
				if(S/a[j].base%(a[j].cnt+1)!=a[j].cnt&&f[S+a[j].base][i+a[j].a]>f[S][i])
				f[S+a[j].base][i+a[j].a]=f[S][i];
				if(F&&f[S][0]>f[S][i]+i*i)
				 f[S][0]=f[S][i]+i*i;
			}
		}
	}
	printf("%d",f[ms][0]-decp);
}