csp-s模拟测试53u,v,w题解

csp-s模拟测试53u,v,w题解,第1张

概述题面:https://www.cnblogs.com/Juve/articles/11602450.html u: 用差分优化修改 二维差分:给(x1,y1),(x2,y2)加上s: $d[x1][y1]+=s,d[x1][y2+1]-=s,d[x2+1][y1]-=s,d[x2+1][y2+1]+=s$ 定义2个差分数组d1,d2,分别记录竖列和斜边的差分 $d1[r][c]+=s,d1[r+l

题面:https://www.cnblogs.com/Juve/articles/11602450.html

u:

用差分优化修改

二维差分:给(x1,y1),(x2,y2)加上s:

$d[x1][y1]+=s,d[x1][y2+1]-=s,d[x2+1][y1]-=s,d[x2+1][y2+1]+=s$

定义2个差分数组d1,d2,分别记录竖列和斜边的差分

$d1[r][c]+=s,d1[r+l-1][c]-=s,d2[r][c+1]-=s,d2[r+l-1][c+l-1]+=s$

统计是求2个前缀和

 

#include<iostream>#include<cstdio>#include<algorithm>#include<cstring>#define int long longusing namespace std;const int MAXN=1e3+5;int n,q,a[MAXN][MAXN],b[MAXN][MAXN],c[MAXN][MAXN],d[MAXN][MAXN],ans=0;int d1[MAXN][MAXN],d2[MAXN][MAXN];signed main(){	scanf("%lld%lld",&n,&q);	while(q--){		int r,c,l,s;		scanf("%lld%lld%lld%lld",&r,&c,&l,&s);		int N=min(n,r+l-1),M=min(n,c+l-1);		d1[r][c]+=s,d1[N+1][c]-=s;		d2[r][c+1]-=s,d2[N+1][M+2]+=s;	}	for(int i=1;i<=n;++i){		for(int j=1;j<=n;++j){			a[i][j]=(a[i-1][j]+d1[i][j]);			b[i][j]=(b[i-1][j-1]+d2[i][j]);			c[i][j]=a[i][j]+b[i][j];			d[i][j]=d[i][j-1]+c[i][j];		}	}	for(int i=1;i<=n;++i){		for(int j=1;j<=n;++j)			ans^=d[i][j];	}	printf("%lld\n",ans);	return 0;}

 

v:

 

搜索每一个情况,对于重复的用map记忆化一下

注意判断长度为奇数的情况,map超时可以用hash表

 

#include<iostream>#include<cstdio>#include<cstring>#include<algorithm>#define re registerusing namespace std;int n,k,sta=0;char ch[35];struct hash_map{	int pre[40000000],cnt;	struct node{		int nxt,to,w;		double v;	}e[6000000];	double &operator[](int sta){		int key=1LL*sta%30000019;		for(int i=pre[key];i;i=e[i].nxt)			if(e[i].to==sta) return e[i].v;		e[++cnt].nxt=pre[key];e[cnt].to=sta;;e[cnt].v=-1.0;		pre[key]=cnt;		return e[cnt].v;	}}mp;inline double dfs(re int x,re int st){	if(x==n-k) return 0.0;	if(mp[st]!=-1) return mp[st];	mp[st]=0;	re double sum=0.0;	for(re int i=1;i<=x>>1;i++){        re int con1=st>>i-1&1,con2=st>>x-i&1;        re int to1=st>>1&~((1<<i-1)-1)|st&(1<<i-1)-1;        re int to2=st>>1&~((1<<x-i)-1)|st&((1<<x-i)-1);        sum+=2.0*max(dfs(x-1,to1)+con1,dfs(x-1,to2)+con2)/x;    }    if(x&1){        re int i=x+1>>1;        re int to=st>>1&~((1<<i-1)-1)|st&(1<<i-1)-1;        re int con=st>>i-1&1;        sum+=(dfs(x-1,to)+con)/x;    }    return mp[st]=sum;}signed main(){	scanf("%d%d",&k);	scanf("%s",ch+1);	for(re int i=1;i<=n;i++)        sta|=(ch[i]==‘W‘)<<n-i;    sta|=1<<n;	printf("%0.7lf\n",dfs(n,sta));	return 0;}

 

w:

如果我们选出一条边,就给边两端的点度数加一,那么第一个答案就是度数为奇数的点个数除以2

我们定义dp[i][0/1]表示i节点,没有/有向父亲的反转i边的两个答案

接下来我们考虑更新:
在更新时我们使用两个参量:p和q,作为更新dp的中间步骤,用p代表不以i为端点做链,q代表以i为端点做链,设i的某个子节点为to,于是有:
p=min(p+dp[to][0],q+dp[to][1])
q=min(p+dp[to][1],q+dp[to][0])
其中p初始化为(0,0),q初始化为(INF,INF)
解释一下:这里pair的add就是对应元素相加(手写!非内置!)
而min *** 作表示先按pair第一关键字比较,再按第二关键字比较
那么这一步就是一个合并的过程:
首先,i不作为链的端点:分两类来合并:如果子节点与i的边翻转了,那么就要累在以i为端点的链里(因为i与父亲的边不翻转,那么i就将是个端点),如果子节点与i的边没有翻转,那么仅针对这棵子树而言,i并没有作为路径的端点,所以更新没有以i为端点链的代价
如果i作为链的端点,同样分两类来合并:如果子节点与i的边翻转了,那么i显然可以成为链的端点,前提是在此之前i并不是链的端点,所以用之前i不是链的端点的代价来更新;反之,如果子节点与i的边没有翻转,那么就此而言i并不是端点,可要求i是链的一个端点,这样就必须用i原先就是链的端点的代价来更新
遍历根节点所有子节点后,更新dp:
如果i与父亲的边没有翻转。即状态dp[i][0]:
首先,i不作为链的端点肯定是一种可能性,直接比较
接着,如果i是链的一个端点,i和父节点的边还没有翻转,那么说明在这个状态下i是真正的奇度点,所以将状态q的first+1后更新
如果i与父亲的边翻转了,同样分两类更新:
首先,i本身作为了链的端点,而由于i本身就是奇度点,所以仅需将q.second+1来更新即可
还有,如果i本身在下面并没有作为链的端点,而i却与父节点的边发生了翻转,所以i就成为了新的奇度点,同时链长还增加了,所以p.first,p.second均增加即可
最后答案即为dp[1][0].first/2,dp[1][0].second

 

#include<iostream>#include<cstdio>#include<cstring>#include<vector>#include<algorithm>using namespace std;const int MAXN=5e5+5;const int inf=0x3f3f3f3f;int n;int to[MAXN<<1],nxt[MAXN<<1],pre[MAXN],flag[MAXN<<1],cnt=0;voID add(int u,int v,int val){	++cnt,to[cnt]=v,nxt[cnt]=pre[u],pre[u]=cnt,flag[cnt]=val;}struct node{	int tim,len;	frIEnd node operator + (node a,node b){		return (node){a.tim+b.tim,a.len+b.len};	}	frIEnd bool operator < (node a,node b){		return a.tim==b.tim?a.len<b.len:a.tim<b.tim;	}}dp[MAXN][2];node min(node a,node b){	return a<b?a:b;}voID dfs(int x,int fa,int fl){	node p={0,0},q={inf,inf};	for(int i=pre[x];i;i=nxt[i]){		int y=to[i];		if(y==fa) continue;		dfs(y,x,flag[i]);		node a=p,b=q;		p=min(a+dp[y][0],b+dp[y][1]);		q=min(b+dp[y][0],a+dp[y][1]);	}	if(fl==2){		dp[x][0]=min(p,q+(node){1,0});		dp[x][1]=min(p+(node){1,1},q+(node){0,1});	}	if(fl==1){		dp[x][0]=(node){inf,inf};		dp[x][1]=min(p+(node){1,1});	}	if(fl==0){		dp[x][0]=min(p,0});		dp[x][1]=(node){inf,inf};	}}signed main(){	scanf("%d",&n);	for(int i=1,a,b,d;i<n;++i){		scanf("%d%d%d%d",&a,&b,&d);		if(d==2){			add(a,2),add(b,2);		}else{			if(c==d) add(a,0),0);			else add(a,1),1);		}	}	dfs(1,0);	printf("%d %d\n",dp[1][0].tim/2,dp[1][0].len);	return 0;}
总结

以上是内存溢出为你收集整理的csp-s模拟测试53u,v,w题解全部内容,希望文章能够帮你解决csp-s模拟测试53u,v,w题解所遇到的程序开发问题。

如果觉得内存溢出网站内容还不错,欢迎将内存溢出网站推荐给程序员好友。

欢迎分享,转载请注明来源:内存溢出

原文地址: http://outofmemory.cn/web/1032557.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
上一篇 2022-05-24
下一篇 2022-05-24

发表评论

登录后才能评论

评论列表(0条)

保存