2018-2019 ACM-ICPC Brazil Subregional Programming Contest 寒假自训第六场

最新推荐文章于 2022-07-14 23:53:51 发布

长沙橘子猫

最新推荐文章于 2022-07-14 23:53:51 发布

阅读量993

点赞数

分类专栏：比赛----gym 数据结构----数链剖分图论----网络流博弈论

本文链接：https://blog.csdn.net/ccsu_cat/article/details/86760911

版权

图论----网络流同时被 3 个专栏收录

16 篇文章 0 订阅

订阅专栏

比赛----gym

14 篇文章 0 订阅

订阅专栏

数据结构----数链剖分

6 篇文章 0 订阅

订阅专栏

今天终于没爆0了，但是速度慢了，漏了C题

付队的题解

B. Marbles

题意：有n对数（xi, yi），每次操作可以选择一个正数u，使得这对数变成（xi-u, yi）或（xi-u, yi-u）或（xi, yi-u），两个玩家轮流操作一次，谁先把一对数变成（0,0）谁赢。

思路：不是把所有数都变成（0，0）才赢，所以不是nim博弈，因此直接写没法写，但是我们转化一下，如果我把一对数变成了（x, x）或（x,0）或（0, x）那我必输，那么我就假设这三种情况不能存在，谁最后无法操作谁输，这样就变成nim博弈了。

#include<bits/stdc++.h>
using namespace std;
const int maxn=105;
int sg[maxn][maxn],vis[maxn*10];
void init(int n)
{
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++)
	{
		if(i==j)continue;
		memset(vis,0,sizeof(vis));
		for(int k=1;k<i;k++)if(k!=j)vis[sg[k][j]]=1;
		for(int k=1;k<j;k++)if(k!=i)vis[sg[i][k]]=1;
		for(int k=1;k<min(i,j);k++)vis[sg[i-k][j-k]]=1;
		for(int k=0;;k++)if(!vis[k])
		{
			sg[i][j]=k;
			break;
		}
	}
}
int main()
{
	init(100);
	int n,x,y,t=0,flag=0;
	cin>>n;
	while(n--)
	{
		cin>>x>>y;
		t^=sg[x][y];
		if(x==y)flag=1;
	}
	if(flag)puts("Y");
	else if(t)puts("Y");
	else puts("N");
}

C. Pizza Cutter

题意：有一个披萨，你要横着切H刀，竖着切V刀，不存在三刀相交在一个点上，问最后披萨被切了多少块。

思路：首先切一刀下去造成的贡献是与之相交的直线+1，那么横竖相交造成的贡献是（H+1）*(V+1)，加下来考虑横横相交的直线，其实不难发现，把每条横线按照y1排序后，形成的新的y2数组的逆序数就是其造成的贡献，竖竖相交同理。

#include<bits/stdc++.h>
#define ll long long
#define low(x) x&-x
using namespace std;
const int maxn=1e5+10;
struct node
{
	int x,y;
	bool operator<(const node&t)const
	{
		if(x==t.x)
		return y<t.y;
		return x<t.x;
	}
}a[maxn];
int b[maxn],c[maxn];
int qu(int x,int n)
{
	int res=0;
	for(;x<=n;x+=low(x))res+=c[x];
	return res;
}
void up(int x)
{
	for(;x;x-=low(x))c[x]++;
}
ll work(int n)
{
	ll res=0;
	memset(c,0,sizeof(c));
	for(int i=1;i<=n;i++)
	scanf("%d%d",&a[i].x,&a[i].y),b[i]=a[i].y;
	sort(a+1,a+1+n);
	sort(b+1,b+1+n);
	for(int i=1;i<=n;i++)
	{
		int k=lower_bound(b+1,b+1+n,a[i].y)-b;
		res+=qu(k,n);
		up(k);
	}
	return res;
}
int main()
{
	int n,m;
	scanf("%d%d",&n,&m);
	scanf("%d%d",&n,&m);
	ll ans=1ll*(n+1)*(m+1);
	ans+=work(n);
	ans+=work(m);
	printf("%lld\n",ans);
}

G. Gasoline

题意：有p个加油站，r个油田，给出每个加油站需要的油以及每个油田库存的油，有C条加油站到油田的带权路，求每个加油站都能加满油时的最长的运输时间的最小值。

思路：二分答案，然后根据答案重新建图跑最大流即可。

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn=2010,inf=1e9;
struct Edge
{
	int from,to;
	ll cap,flow;
};
struct Dinic
{
	int n,m,s,t;
	vector<Edge>edges;
	vector<int>G[maxn];
	bool vis[maxn];
	int d[maxn],cur[maxn];
	void init(int n)
	{
		this->n=n;
		for(int i=0;i<n;i++)
		G[i].clear();
		edges.clear();
	}
	void AddEdge(int from,int to,ll cap)
	{
		edges.push_back((Edge){from,to,cap,0});
		edges.push_back((Edge){to,from,0,0});
		m=edges.size();
		G[from].push_back(m-2);
		G[to].push_back(m-1);
	}
	bool bfs()
	{
		memset(vis,0,sizeof(vis));
		queue<int>Q;
		Q.push(s);
		d[s]=0;
		vis[s]=1;
		while(!Q.empty())
		{
			int x=Q.front();Q.pop();
			for(int i=0;i<G[x].size();i++)
			{
				Edge& e=edges[G[x][i]];
				if(!vis[e.to]&&e.cap>e.flow)
				{
					vis[e.to]=1;
					d[e.to]=d[x]+1;
					Q.push(e.to);
				}
			}
		}
		return vis[t];
	}
	ll dfs(int x,ll a)
	{
		if(x==t||a==0)return a;
		ll flow=0,f;
		for(int& i=cur[x];i<G[x].size();i++)
		{
			Edge& e=edges[G[x][i]];
			if(d[x]+1==d[e.to]&&(f=dfs(e.to,min(a,e.cap-e.flow)))>0)
			{
				e.flow+=f;
				edges[G[x][i]^1].flow-=f;
				flow+=f;
				a-=f;
				if(a==0)break;
			}
		}
		return flow;
	}
	ll Maxflow(int s,int t)
	{
		this->s=s,this->t=t;
		ll flow=0;
		while(bfs())
		{
			memset(cur,0,sizeof(cur));
			flow+=dfs(s,inf);
		}
		return flow;
	}
}solve;
int u[20005],v[20005],w[20005],a[1005],b[1005];
int main()
{
	int p,r,c,sum=0;
	scanf("%d%d%d",&p,&r,&c);
	int S=0,T=p+r+1;
	for(int i=1;i<=p;i++)scanf("%d",&a[i]),sum+=a[i];
	for(int i=1;i<=r;i++)scanf("%d",&b[i]);
	for(int i=1;i<=c;i++)
	scanf("%d%d%d",&u[i],&v[i],&w[i]);
	int L=0,R=1e7,mid,all=1e7;
	while(L<R)
	{
		mid=(L+R)/2;
		solve.init(T+1);
		for(int i=1;i<=p;i++)
		solve.AddEdge(S,i,a[i]);
		for(int i=1;i<=r;i++)
		solve.AddEdge(p+i,T,b[i]);
		for(int i=1;i<=c;i++)
		if(w[i]<=mid)
		solve.AddEdge(u[i],p+v[i],inf);
		ll res=solve.Maxflow(S,T);
		if(res==sum)
		R=mid;
		else L=mid+1;
	}
	if(L==all)puts("-1");
	else printf("%d\n",L);
}

L. Subway Lines

题意：给一棵树，询问两条路径上公共点的个数。

思路：树剖，将第一条路径的所有点用线段树标记，然后查询第二条路径在线段树对应区间的和即可。

#include<bits/stdc++.h>
using namespace std;
const int maxn=2e5+10;
int id[maxn],dep[maxn],cnt;
int top[maxn],size[maxn],f[maxn],son[maxn];
vector<int>G[maxn];
int sum[maxn*4],tag[maxn*4];
void dfs1(int u,int fa,int deep)
{
	f[u]=fa;
	size[u]=1;
	dep[u]=deep;
	for(int i=0;i<G[u].size();i++)
	{
		int v=G[u][i];
		if(v==fa)continue;
		dfs1(v,u,deep+1);
		size[u]+=size[v];
		if(size[v]>size[son[u]])son[u]=v;
	}
}
void dfs2(int u,int root)
{
	top[u]=root;
	id[u]=++cnt;
	if(son[u])dfs2(son[u],root);
	for(int i=0;i<G[u].size();i++)
	{
		int v=G[u][i];
		if(v==f[u]||v==son[u])continue;
		dfs2(v,v);
	}
}
void pushdown(int o,int ls,int rs,int l,int m,int r)
{
	if(tag[o]!=-1&&l!=r)
	{
		sum[ls]=(m-l+1)*tag[o];
		sum[rs]=(r-l)*tag[o];
		tag[ls]=tag[rs]=tag[o];
		tag[o]=-1;
	}	
}
void update(int o,int l,int r,int ql,int qr,int v)
{	
	int m=(l+r)/2,ls=o*2,rs=o*2+1;
	pushdown(o,ls,rs,l,m,r);
	if(l>=ql&&r<=qr)
	{
		sum[o]=v*(r-l+1);
		tag[o]=v;
		return;
	}
	if(ql<=m)update(ls,l,m,ql,qr,v);
	if(qr>m)update(rs,m+1,r,ql,qr,v);
	sum[o]=sum[ls]+sum[rs];
}
int query(int o,int l,int r,int ql,int qr)
{
	int m=(l+r)/2,ls=o*2,rs=o*2+1,res=0;
	pushdown(o,ls,rs,l,m,r);
	if(sum[o]==r-l+1)
	return min(r,qr)-max(l,ql)+1;	
	if(sum[o]==0)return 0;
	if(l>=ql&&r<=qr)
	return sum[o];	
	if(ql<=m)res+=query(ls,l,m,ql,qr);
	if(qr>m)res+=query(rs,m+1,r,ql,qr);
	return res;
}
void up_path(int x,int y,int z)
{
	while(top[x]!=top[y])
	{
		if(dep[top[x]]<dep[top[y]])swap(x,y);
		update(1,1,cnt,id[top[x]],id[x],z);
		x=f[top[x]];
	}
	if(id[x]>id[y])swap(x,y);
	update(1,1,cnt,id[x],id[y],z);
}
int qu_path(int x,int y)
{
	int ans=0;
	while(top[x]!=top[y])
	{
		if(dep[top[x]]<dep[top[y]])swap(x,y);
		ans+=query(1,1,cnt,id[top[x]],id[x]);
		x=f[top[x]];
	}
	if(id[x]>id[y])swap(x,y);
	ans+=query(1,1,cnt,id[x],id[y]);
	return ans;
}
int main()
{
	int n,u,v,a,b,c,d,q;
	scanf("%d%d",&n,&q);
	for(int i=1;i<n;i++)
	{
		scanf("%d%d",&u,&v);
		G[u].push_back(v);
		G[v].push_back(u);
	}
	memset(tag,-1,sizeof(tag));
	dfs1(1,0,1);
	dfs2(1,1);
	while(q--)
	{
		scanf("%d%d%d%d",&a,&b,&c,&d);
		up_path(a,b,1);
		printf("%d\n",qu_path(c,d));
		up_path(a,b,0);
	}
}

长沙橘子猫

关注

0
点赞
踩
1

收藏

觉得还不错? 一键收藏
打赏
0
评论
2018-2019 ACM-ICPC Brazil Subregional Programming Contest 寒假自训第六场

今天终于没爆0了，但是速度慢了，漏了C题付队的题解B. Marbles题意：有n对数（xi, yi），每次操作可以选择一个正数u，使得这对数变成（xi-u, yi）或（xi-u, yi-u）或（xi, yi-u），两个玩家轮流操作一次，谁先把一对数变成（0,0）谁赢。思路：不是把所有数都变成（0，0）才赢，所以不是nim博弈，因此直接写没法写，但是我们转化一下，如果我把一对数变成了...
复制链接

扫一扫