hdu 6251 Inkopolis [Tarjan+想法]

题意:给你一个基环树,树上的边有各种颜色,每次操作有两种:

①修改一条边的颜色

②询问这个图有多少色块(相同颜色的边公用一个点,则为一块)

题解:(由于一遍一遍的改,代码很不精简=。=)

首先,我们可以考虑修改树上的边那么有7种情况:(假设边原来的颜色为a,修改后为b)

①边的两端的都存在颜色为a的边,修改后边的两端都存在颜色为b的边,则贡献为0。

②边的两端的都存在颜色为a的边,修改后边的两端有一端存在颜色为b的边,则贡献为1。

③边的两端的都存在颜色为a的边,修改后边的两端都不存在颜色为b的边,则贡献为2。

④边的两端有一端存在颜色为a的边,修改后边的两端都存在颜色为b的边,则贡献为-1。

⑤边的两端都不存在颜色为a的边,修改后边的两端都存在颜色为b的边,则贡献为-2。

⑥边的两端都不存在颜色为a的边,修改后边的两端有一端存在颜色为b的边则贡献为-1。

⑦边的两端有一端存在颜色为a的边,修改后边的两端都不存在颜色为b的边,则贡献为1。

其他情况贡献均为0。

然后我们考虑环上的情况

假如环是如下情况:

当加入颜色不是1的边时需要比之前条件的贡献-1

假如环时如下情况

当加入颜色是1的边时需要比之前条件的贡献+1

最后,与点相邻的边就直接加入set维护,Tarjan找基环,mp记录边的情况mp<0代表该边时环边。

AC代码:

#include<stdio.h>
#include<string.h>
#include<set>
#include<queue>
#include<map>
#include<stdlib.h>
#include<algorithm>
#include<vector>
#define N 200005
using namespace std;
inline char nc(){
    static char buf[100000],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline int _read(){
    char ch=nc();int sum=0;
    while(!(ch>='0'&&ch<='9'))ch=nc();
    while(ch>='0'&&ch<='9')sum=sum*10+ch-48,ch=nc();
    return sum;
}

struct edge
{
	int to,next,col;
	edge(){}
	edge(int to,int next,int col)
	{
		this->to=to;
		this->next=next;
		this->col=col;
	}
}ed[N*2];
int head[N],dfn[N],mark[N*2],cir[N],sta[N],low[N];
int lnum,index,tot,top,sum;
void addline(int a,int b,int c)
{
	ed[lnum]=edge(b,head[a],c);
	head[a]=lnum++;
}
void dfs(int u)
{
	dfn[u]=low[u]=++index;
	sta[top++]=u;
	for(int i=head[u];~i;i=ed[i].next)
	{
		int to=ed[i].to;
		if(mark[i])continue;
		mark[i]=mark[i^1]=1;
		if(!dfn[to])
		{
			dfs(to);
			low[u]=min(low[u],low[to]);
		}
		else low[u]=min(low[u],dfn[to]);
	}
	if(low[u]==dfn[u])
	{
		int now=tot,k;
		do
		{
			k=sta[--top];
			cir[tot++]=k;
		}while(k!=u);
		if(tot-now==1)tot=now;
	}	
}
struct node
{
	int col,num;
	node(){}
	node(int col,int num)
	{
		this->col=col;
		this->num=num;
	}
};
multiset<int>st[N];
multiset<int>::iterator it2;
multiset<node>ci;
multiset<node>::iterator it,haha;
map<int,int>mp[N];
bool operator<(node a,node b)
{
	return a.col<b.col;
}
int Get(int a,int col)
{
	it2=st[a].find(col);
	if(it2==st[a].end())return 0;
	it2++;
	if(it2==st[a].end()||*it2!=col)return 0;
	return 1;
}
int judge(int a,int b,int c)
{
	int col=abs(mp[a][b]);
	int f1=Get(a,col);
	int f2=Get(b,col);
	int f3=st[a].find(c)!=st[a].end();
	int f4=st[b].find(c)!=st[b].end();
	st[a].erase(st[a].find(col));
	st[b].erase(st[b].find(col));
	st[a].insert(c);
	st[b].insert(c);
	if(f1&&f2&&f3&&f4)return 0;
	if(f1&&f2&&(f3||f4))return 1;
	if(f1&&f2&&!f3&&!f4)return 2;
	if((f1||f2)&&f3&&f4)return -1;
	if(!f1&&!f2&&f3&&f4)return -2;
	if(!f1&&!f2&&(f3||f4))return -1;
	if((f1||f2)&&!f3&&!f4)return 1;
	return 0;
}
int q1[N],q2[N],q3[N];
int main()
{
	int T;
	T=_read();
	int cas=1;
	while(T--)
	{
		memset(head,-1,sizeof(head));
		memset(mark,0,sizeof(mark));
		memset(dfn,0,sizeof(dfn));
		memset(low,0,sizeof(low));
		lnum=sum=top=tot=index=0;
		ci.clear();
		int n,q;
		n=_read();q=_read();
		for(int i=1;i<=n;i++)st[i].clear(),mp[i].clear();
		for(int i=0;i<n;i++)
		{
			q1[i]=_read();q2[i]=_read();q3[i]=_read();
			addline(q1[i],q2[i],N+1);
			addline(q2[i],q1[i],N+1);
			st[q1[i]].insert(N+1);
			st[q2[i]].insert(N+1);
			mp[q1[i]][q2[i]]=mp[q2[i]][q1[i]]=N+1;
		}
		sta[top++]=1;
		dfs(1);
		for(int i=0;i<tot;i++)
		{
			mp[cir[i]][cir[(i+1)%tot]]*=-1;
			for(int j=head[cir[i]];~j;j=ed[j].next)
			{
				int to=ed[j].to;
				if(to==cir[(i+1)%tot])
				{
					it=ci.find(node(ed[j].col,0));
					if(it==ci.end())
					{
						ci.insert(node(ed[j].col,1));
						continue;
					}
					node k=*it;
					k.num++;
					ci.erase(it);
					ci.insert(k);
				}
			}
		}
		for(int i=0;i<tot;i++)
			mp[cir[i]][cir[(i-1+tot)%tot]]*=-1;
		memset(mark,0,sizeof(mark));	
		printf("Case #%d:\n",cas++);
		sum=1;
		for(int gg=0;gg<q+n;gg++)
		{
			int a,b,c;
			if(gg<n)a=q1[gg],b=q2[gg],c=q3[gg];
			else a=_read(),b=_read(),c=_read();
			if(abs(mp[a][b])==c)
			{
				if(gg>=n)printf("%d\n",sum);
				continue;
			}
			sum+=judge(a,b,c);
			if(mp[a][b]<0)
			{
				it=ci.find(node(abs(mp[a][b]),0));
				node k=*it;
				if(k.num==tot&&k.col!=c)
				{
					sum--;
					k.num--;
					ci.erase(it);
					ci.insert(k);
					ci.insert(node(c,1));
				}
				else if(k.num==1&&ci.size()==2)
				{
					haha=it;
					if(haha==ci.begin())haha++;
					else haha--;
					node g=*haha;
					if(g.col==c)sum++;
					ci.erase(it);
					it=ci.find(node(c,0));
					if(it==ci.end())ci.insert(node(c,1));
					else 
					{
						g=*it;
						g.num++;
						ci.erase(it);
						ci.insert(g);
					}
				}
				else 
				{
					k.num--;
					ci.erase(it);
					ci.insert(k);
					it=ci.find(node(c,0));
					if(it==ci.end())ci.insert(node(c,1));
					else 
					{
						k=*it;
						k.num++;
						ci.erase(it);
						ci.insert(k);
					}
				}
				mp[a][b]=mp[b][a]=-c;
			}
			else mp[a][b]=mp[b][a]=c;
			if(gg>=n)printf("%d\n",sum);
		}
	}
} 

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值