码迷，mamicode.com
首页 > 其他好文 > 详细

csp-s模拟45，50 weight，蔬菜，联盟

时间：2019-09-24 11:56:15      阅读：78      评论：0      收藏：0      [点我收藏+]

标签：put   uil   水过   ems   http   efi   nod   csp   try   

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

蔬菜：

题解说和四维偏序有关，但是看到它没有修改，可以莫队水过

一个二维莫队，定义四个指针，分别表示左上角和右下角的坐标

然后模拟序列上的莫队移动

while(u<ask[i].x) del_h(u++);
while(u>ask[i].x) add_h(--u);
while(d<ask[i].p) add_h(++d);
while(d>ask[i].p) del_h(d--);
while(l<ask[i].y) del_l(l++);
while(l>ask[i].y) add_l(--l);
while(r<ask[i].q) add_l(++r);
while(r>ask[i].q) del_l(r--);

u是up，d是down，l是left，r是right，然后统计答案即可

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
#define int long long
#define re register
using namespace std;
const int MAXN=205;
int n,m,q,a[MAXN][MAXN],cnt=0,tot=0,ans[100005];
int b[MAXN*MAXN],num[MAXN*MAXN],u=1,d=0,l=1,r=0;
struct node{
	int x,y,p,q,id;
	friend bool operator < (node a,node b){
		return (a.x==b.x?(a.y==b.y?(a.p==b.p?(a.q<b.q):a.p<b.p):a.y<b.y):a.x<b.x);
	}
}ask[100005];
void del_h(int ha){
	for(int i=l;i<=r;++i) --num[a[ha][i]];
}
void add_h(int ha){
	for(int i=l;i<=r;++i) ++num[a[ha][i]];
}
void del_l(int li){
	for(int i=u;i<=d;++i) --num[a[i][li]];
}
void add_l(int li){
	for(int i=u;i<=d;++i) ++num[a[i][li]];
}
signed main(){
	//freopen("test.in","r",stdin);
	//freopen("vio.out","w",stdout);
	scanf("%lld%lld%lld",&n,&m,&q);
	for(re int i=1;i<=n;++i){
		for(re int j=1;j<=m;++j){
			scanf("%lld",&a[i][j]);
			b[++tot]=a[i][j];
		}
	}
	sort(b+1,b+tot+1);
	cnt=unique(b+1,b+tot+1)-b-1;
	for(re int i=1;i<=n;++i){
		for(re int j=1;j<=m;++j)
			a[i][j]=lower_bound(b+1,b+cnt+1,a[i][j])-b;
	}
	for(int i=1;i<=q;++i){
		scanf("%lld%lld%lld%lld",&ask[i].x,&ask[i].y,&ask[i].p,&ask[i].q);
		ask[i].id=i;
	}
	sort(ask+1,ask+q+1);
	for(int i=1;i<=q;++i){
		while(u<ask[i].x) del_h(u++);
		while(u>ask[i].x) add_h(--u);
		while(d<ask[i].p) add_h(++d);
		while(d>ask[i].p) del_h(d--);
		while(l<ask[i].y) del_l(l++);
		while(l>ask[i].y) add_l(--l);
		while(r<ask[i].q) add_l(++r);
		while(r>ask[i].q) del_l(r--);
		for(int j=1;j<=cnt;++j){
			ans[ask[i].id]+=num[j]*num[j];
		}
	}
	for(int i=1;i<=q;++i){
		printf("%lld\n",ans[i]);
	}
	return 0;
}

联盟：

这道题和树的直径有关

先求出原树的直径，然后枚举断开的直径上的边（段直径上的边一定最优）

设断开之后两个联通块的直径分别是$l_1$，$l_2$，那么第一问答案就是

$max(l_1,l_2,\lceil\frac{l_1}{2}\rceil+\lceil\frac{l_2}{2}\rceil+1)$

第一问出来了第二问就不是问题了

然后第三问，我们随便找出一个要断的边，连接的点就是两个联通块的直径的中点

打了4个dfs和4个bfs，码长4.0k，好像我的做法很麻烦？

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#define re register
using namespace std;
const int MAXN=3e5+5;
int n,fmx=0,smx=0,ans=0x3f3f3f3f;
int to[MAXN<<1],nxt[MAXN<<1],pre[MAXN],cnt=1,id[MAXN<<1],fr[MAXN<<1];
inline void add(re int u,re int v,re int rk){
	++cnt,to[cnt]=v,nxt[cnt]=pre[u],pre[u]=cnt,id[cnt]=rk,fr[cnt]=u;
}
int dp[2][MAXN];
int DFs(int x,int fa,int flag){
	int maxx=0;
	dp[flag][x]=0;
	for(int i=pre[x];i;i=nxt[i]){
		int y=to[i];
		if(y==fa) continue;
		int p=DFs(y,x,flag);
		dp[flag][x]=max(dp[flag][x],max(maxx+p,dp[flag][y]));
		maxx=max(maxx,p);
	}
	return maxx+1;
}
int d[MAXN],p[MAXN],qd,zd,dis[MAXN];
int stac[MAXN],topc=0;
int bfs(int s){
	topc=0;
	memset(d,0x3f,sizeof(d));
	d[s]=0;
	p[s]=0;
	queue<int>q;
	q.push(s);
	while(!q.empty()){
		int x=q.front();
		q.pop();
		for(int i=pre[x];i;i=nxt[i]){
			int y=to[i];
			if(d[y]==0x3f3f3f3f){
				d[y]=d[x]+1;
				p[y]=i;
				q.push(y);
			}
		}
	}
	int y=1;
	for(int i=1;i<=n;++i){
		if(dis[i]==0x3f3f3f3f) continue;
		if(d[i]>d[y]) y=i,topc=0;
		//if(d[i]==d[y]) stac[++topc]=i;
	}
	return y;
}
int BFS(int s,int ed){
	if(s==ed) return 0;
	memset(dis,0x3f,sizeof(dis));
	dis[s]=0;
	queue<int>q;
	q.push(s);
	while(!q.empty()){
		int x=q.front();
		q.pop();
		for(int i=pre[x];i;i=nxt[i]){
			int y=to[i];
			if(y==ed) continue;
			if(dis[y]==0x3f3f3f3f){
				dis[y]=dis[x]+1;
				q.push(y);
			}
		}
	}
	int res=0;
	for(int i=1;i<=n;++i){
		if(dis[i]==0x3f3f3f3f) continue;
		res=max(res,dis[i]);
	}
	return res;
}
int sta[MAXN],top=0;
/*void work(int x){
	if(x==0) return ;
	int i=p[x];
	int p=BFS(zd,fr[i]),q=BFS(qd,x);
	//cout<<x<<‘ ‘<<fr[i]<<‘ ‘<<p<<‘ ‘<<q<<‘ ‘<<‘ ‘<<(p+1)/2+(q+1)/2+1<<‘ ‘<<i<<endl;
	int t=max(max(p,q),(p+1)/2+(q+1)/2+1);
	if(ans==t) sta[++top]=i;
	if(ans>t) top=0,sta[++top]=i;
	ans=min(ans,t);
	work(fr[i]);
}*/
void work(int x){
	if(x==0) return ;
	int i=p[x];
	int p=dp[1][fr[i]],q=dp[0][x];
	//cout<<x<<‘ ‘<<fr[i]<<‘ ‘<<p<<‘ ‘<<q<<‘ ‘<<‘ ‘<<(p+1)/2+(q+1)/2+1<<‘ ‘<<i<<endl;
	int t=max(max(p,q),(p+1)/2+(q+1)/2+1);
	if(ans==t) sta[++top]=i;
	if(ans>t) top=0,sta[++top]=i;
	ans=min(ans,t);
	work(fr[i]);
}
int Bfs(int s){
	memset(d,0x3f,sizeof(d));
	d[s]=0;
	p[s]=0;
	queue<int>q;
	q.push(s);
	while(!q.empty()){
		int x=q.front();
		q.pop();
		for(int i=pre[x];i;i=nxt[i]){
			int y=to[i];
			if(d[y]==0x3f3f3f3f){
				d[y]=d[x]+1;
				p[y]=i;
				q.push(y);
			}
		}
	}
	int y=1;
	for(int i=1;i<=n;++i){
		if(dis[i]==0x3f3f3f3f) continue;
		if(d[i]>d[y]) y=i;
	}
	return y;
}
int edge,frr,too,len;
int BFs(int s,int ed){
	if(s==ed) return 0;
	memset(dis,0x3f,sizeof(dis));
	dis[s]=0;
	queue<int>q;
	q.push(s);
	while(!q.empty()){
		int x=q.front();
		q.pop();
		for(int i=pre[x];i;i=nxt[i]){
			int y=to[i];
			if(y==ed) continue;
			if(dis[y]==0x3f3f3f3f){
				dis[y]=dis[x]+1;
				p[y]=i;
				q.push(y);
			}
		}
	}
	int y=s;
	for(int i=1;i<=n;++i){
		if(dis[i]==0x3f3f3f3f) continue;
		if(dis[i]>dis[y]) y=i;
	}
	return y;
}
int ans1=0,ans2=0;
void dfs(int x,int l){
	if(x==0) return ;
	if(l==len){
		ans1=x;
		return ;
	}
	int i=p[x];
	dfs(fr[i],l+1);
}
void DFS(int x,int l){
	if(x==0) return ;
	if(l==len){
		ans2=x;
		return ;
	}
	int i=p[x];
	DFS(fr[i],l+1);
}
signed main(){
	//freopen("dt.in","r",stdin);
	//freopen("my.out","w",stdout);
	scanf("%d",&n);
	for(int i=1,u,v;i<n;++i){
		scanf("%d%d",&u,&v);
		add(u,v,i),add(v,u,i);
	}
	qd=bfs(1);
	zd=bfs(qd);
	DFs(qd,0,0);
	DFs(zd,0,1);
	work(zd);
	printf("%d\n",ans);
	sort(sta+1,sta+top+1);
	top=unique(sta+1,sta+top+1)-sta-1;
	printf("%d ",top);
	for(int i=1;i<=top;++i){
		printf("%d ",id[sta[i]]);
	}
	puts("");
	edge=sta[top];
	frr=fr[edge],too=to[edge];
	qd=BFs(frr,too);
	zd=BFs(qd,too);
	len=dis[zd]/2;
	//cout<<len<<endl;
	dfs(zd,0);
	qd=BFs(too,frr);
	zd=BFs(qd,frr);
	len=dis[zd]/2;
	DFS(zd,0);
	printf("%d %d %d %d\n",frr,too,ans1,ans2);
	return 0;
}

weight：

先kruskal跑出最小生成树，然后对于树边和非树边进行讨论

对于一条非树边,我们至少要将它的权值调整到树上这两个端点对应路径边权最大值 -1 才可以,

否则我们一定可以不选这条边。显然,我们调整到这么大也足够了。

对于一条树边,我们关心的显然是两个端点对应的简单路径经过这条树边的那些边,

我们最大的可能选择是那些边中权值最小的边的权值 -1, (否则我们可以选那条最小边而不选这条树边),

而我们如果将这条树边的边权调整成那个值,它也一定还会在最小生成树中。

然后我们对于最小生成树进行树剖来实现操作

 

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int MAXN=1e6+5;
const int inf=0x3f3f3f3f;
int n,m,xkl,ans[MAXN];
struct EDGE{
	int fr,to,nxt,val,id;
	bool flag;
	friend bool operator < (EDGE a,EDGE b){
		return a.val<b.val;
	}
}ed[MAXN];
int cnt=0,pre[MAXN],head[MAXN],tot=0;
void add(int u,int v,int w,int id){
	ed[++cnt]=(EDGE){u,v,pre[u],w,id,0};
	pre[u]=cnt;
}
struct node{
	int fr,to,nxt,val,id;
}e[MAXN];
void ADD(int u,int v,int val,int id){
	e[++tot]=(node){u,v,head[u],val,id};
	head[u]=tot;
}
int fa[MAXN];
int find(int x){
	return fa[x]=(fa[x]==x?x:find(fa[x]));
}
void kruskal(){
	for(int i=1;i<=n;++i) fa[i]=i;
	sort(ed+1,ed+cnt+1);
	int sum=0;
	for(int i=1;i<=cnt;++i){
		int x=find(ed[i].fr),y=find(ed[i].to);
		if(x!=y){
			fa[x]=y;
			++sum;
			ed[i].flag=1;
			ADD(ed[i].fr,ed[i].to,ed[i].val,ed[i].id);
			ADD(ed[i].to,ed[i].fr,ed[i].val,ed[i].id);
			if(sum==n-1) break;
		}
	}
}
int deep[MAXN],size[MAXN],son[MAXN],val[MAXN],rcd[MAXN];
void dfs(int x,int father){
	size[x]=1;
    for(int i=head[x];i;i=e[i].nxt){
        int y=e[i].to;
        if(y==father) continue;
        deep[y]=deep[x]+1;
        fa[y]=x;
        val[y]=e[i].val;
        rcd[y]=e[i].id;
        dfs(y,x);
        size[x]+=size[y];
        if(size[y]>size[son[x]]) son[x]=y;
    }
}
int top[MAXN],id[MAXN],rk[MAXN],dfs_order=0;
void DFS(int x,int chain){
	top[x]=chain;
    id[x]=++dfs_order;
    rk[dfs_order]=x;
    if(son[x]) DFS(son[x],chain);
    for(int i=head[x];i;i=e[i].nxt){
        int y=e[i].to;
        if(y==fa[x]||y==son[x]) continue;
        DFS(y,y);
    }
}
struct TREE{
	int l,r,laz,mx,mn;
}tr[MAXN<<2];
void build(int k,int l,int r){
	tr[k].l=l,tr[k].r=r,tr[k].laz=inf;
	if(l==r){
		tr[k].mn=inf;
		tr[k].mx=val[rk[l]];
		return ;
	}
	int mid=(l+r)>>1;
	build(k<<1,l,mid),build(k<<1|1,mid+1,r);
	tr[k].mx=max(tr[k<<1].mx,tr[k<<1|1].mx);
}
void down(int k){
    tr[k<<1].laz=min(tr[k<<1].laz,tr[k].laz);
    tr[k<<1|1].laz=min(tr[k<<1|1].laz,tr[k].laz);
    tr[k<<1].mn=min(tr[k].laz,tr[k<<1].mn);
    tr[k<<1|1].mn=min(tr[k].laz,tr[k<<1|1].mn);
    tr[k].laz=inf;
}
int query(int k,int opl,int opr){
	int l=tr[k].l,r=tr[k].r;
    if(opl<=l&&r<=opr){
        return tr[k].mx;
    }
    if(tr[k].laz!=inf) down(k);
    int mid=(l+r)>>1,res=0;
    if(opl<=mid) res=max(res,query(k<<1,opl,opr));
    if(opr>mid) res=max(res,query(k<<1|1,opl,opr));
	return res;
}
void update(int k,int opl,int opr,int val){
	int l=tr[k].l,r=tr[k].r;
	if(opl<=l&&r<=opr){
		tr[k].mn=min(tr[k].mn,val);
		tr[k].laz=min(tr[k].laz,val);
		return ;
	}
	if(tr[k].laz!=inf) down(k);
	int mid=(l+r)>>1;
	if(opl<=mid) update(k<<1,opl,opr,val);
    if(opr>mid) update(k<<1|1,opl,opr,val);
	tr[k].mn=min(tr[k<<1].mn,tr[k<<1|1].mn);
}
int query_path(int x,int y){
	int res=0;
    while(top[x]!=top[y]){
        if(deep[top[x]]<deep[top[y]]) swap(x,y);
        res=max(res,query(1,id[top[x]],id[x]));
        x=fa[top[x]];
    }
    if(deep[x]>deep[y]) swap(x,y);
    res=max(res,query(1,id[x]+1,id[y]));
    return res;
}
void update_path(int x,int y,int val){
	while(top[x]!=top[y]){
        if(deep[top[x]]<deep[top[y]]) swap(x,y);
        update(1,id[top[x]],id[x],val);
        x=fa[top[x]];
    }
    if(deep[x]>deep[y]) swap(x,y);
    update(1,id[x]+1,id[y],val);
}
void ask_ans(int k){
	if(tr[k].l==tr[k].r){
        ans[rcd[rk[tr[k].l]]]=tr[k].mn-1;
        if(ans[rcd[rk[tr[k].l]]]==inf-1){
            ans[rcd[rk[tr[k].l]]]=-1;
        }
        return ;
    }
	if(tr[k].laz!=inf) down(k);
	ask_ans(k<<1),ask_ans(k<<1|1);
}
signed main(){
	//freopen("test.in","r",stdin);
	//freopen("vio.out","w",stdout);
	scanf("%d%d%d",&n,&m,&xkl);
	for(int i=1,u,v,w;i<=m;++i){
		scanf("%d%d%d",&u,&v,&w);
		add(u,v,w,i);
		//add(v,u,w,i);
	}
	kruskal();
	memset(fa,0,sizeof(fa));
	deep[1]=1;
	dfs(1,0);DFS(1,1);
	build(1,1,n);
	for(int i=1;i<=cnt;++i){
		if(ed[i].flag==0){
			//cout<<ed[i].id<<endl;
			ans[ed[i].id]=query_path(ed[i].fr,ed[i].to)-1;
			//cout<<ans[ed[i].id]<<endl;
			//cout<<ed[i].fr<<‘ ‘<<ed[i].to<<endl;
			update_path(ed[i].fr,ed[i].to,ed[i].val);
		}
	}
	ask_ans(1);
	for(int i=1;i<=m;++i)
		printf("%d ",ans[i]);
	puts("");
	return 0;
}

 

csp-s模拟45，50 weight，蔬菜，联盟

标签：put   uil   水过   ems   http   efi   nod   csp   try   

原文地址：https://www.cnblogs.com/Juve/p/11576165.html
踩
(0)
赞
(0)
   
举报
评论 一句话评论（0）
登录后才能评论！
分享档案
更多>
2021年07月29日 (22)
2021年07月28日 (40)
2021年07月27日 (32)
2021年07月26日 (79)
2021年07月23日 (29)
2021年07月22日 (30)
2021年07月21日 (42)
2021年07月20日 (16)
2021年07月19日 (90)
2021年07月16日 (35)
周排行
mamicode.com排行更多图片
更多
友情链接
兰亭集智   国之画   百度统计   站长统计   阿里云   chrome插件   新版天听网
关于我们 - 联系我们 - 留言反馈
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码！