图论专题qwq

拓扑排序

#include<bits/stdc++.h>
#define N 100005
using namespace std;
int n,m,x,y;
int in[N];
vector<int>v[N];
queue<int>q;
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)v[i].clear();
	memset(in,0,sizeof(in));
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&x,&y);
		in[y]++;
		v[x].push_back(y);
	}
	for(int i=1;i<=n;i++) if(!in[i]) q.push(i);
	while(!q.empty())
	{
		int now=q.front();
		q.pop();
		cout<<now<<" ";
		for(int i=0;i<v[now].size();i++)
		{
			int next=v[now][i];
			in[next]--;
			if(in[next]==0) q.push(next);
		}
	}
	return 0;
}

欧拉回路

#include<bits/stdc++.h>
using namespace std;
const int maxn=100010;
int n,m,u,v;
int du[maxn][2],d[maxn];
stack<int> st;
vector<int> q[maxn];
void dfs(int x)
{
	for(int i=d[x];i<q[x].size();i=d[x])
	{
		d[x]=i+1;
		dfs(q[x][i]);
	}
	st.push(x);
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&u,&v);
		q[u].push_back(v);
		du[u][1]++;
		du[v][0]++;
	}
	for(int i=1;i<=n;i++) sort(q[i].begin(),q[i].end());
	int cnt1=0,cnt2=0,book=1;
	bool flag=1;
	for(int i=1;i<=n;i++)
	{
		if(du[i][1]!=du[i][0]&&i!=1) flag=0;
		if(du[i][1]-du[i][0]==1)
		{
			cnt2++;
			book=i;
		}
		if(du[i][0]-du[i][1]==1)cnt1++;
	}
	if(flag==0&&!(cnt1==cnt2&&cnt1==1))
	{
		printf("No");
		return 0;
	}
	dfs(book);
	while(!st.empty())
	{
		printf("%d ",st.top());
		st.pop();
	}
	return 0;
}

最短路

ＳＰＦＡ：

dis[1]=0;
queue<int>q;///存的是松弛过的顶点
q.push(1);
visit[1]=true;
while(!q.empty())
{
	int u=q.front(); 
	q.pop();
	visit[u]=false; 
	for(int E=head[u]; E; E=e[E].next) 
	{
		int v=e[E].to; 
		if(dis[v]>dis[u]+e[E].w) 
		{
			dis[v]=dis[u]+e[E].w;
			if(!visit[v]) 
			{
				visit[v]=true;
				q.push(v);
			}
		}
	}
}

D i j ：

#include<bits/stdc++.h>
using namespace std;
const int N=100010,M=500010;
struct edge{int to,next,w;} e[M];
int head[N],dis[N];
int n,m,k,cnt;
bool vis[M];
void add_edge(int u,int v,int w)
{
	cnt++;
	e[cnt].w=w;
	e[cnt].to=v;
	e[cnt].next=head[u];
	head[u]=cnt;
}
struct node
{
	int dis,pos;
	bool operator <( const node &x ) const
	{ return x.dis<dis; }
};
priority_queue<node> q;
void dijkstra()
{
	dis[k]=0;
	q.push((node){0,k});
	while(!q.empty())
	{
		node tmp=q.top();
		q.pop();
		int x=tmp.pos;
		if(vis[x]) continue;
		vis[x]=1;
		for(int i=head[x];i;i=e[i].next)
		{
			int y=e[i].to;
			if(dis[y]>dis[x]+e[i].w)
			{
				dis[y]=dis[x]+e[i].w;
				if(!vis[y]) q.push((node){dis[y],y});
			}
		}
	}
}

int main()
{
	scanf("%d%d%d",&n,&m,&k);
	for(int i=1;i<=n;++i) dis[i]=0x7fffffff;
	for(int i=1;i<=m;++i)
	{
		int u,v,w;
		scanf("%d%d%d",&u,&v,&w);
		add_edge(u,v,w);
	}
	dijkstra();
	for(int i=1;i<=n;i++) printf("%d ",dis[i]);
	return 0;
}

tarjan:

tarjan2_xiyuping24的博客-CSDN博客

2-SAT:【模板】2-SAT 问题 

逻辑上的关系（a^b）<==>（非a=>b）并且（非b=>a）把a,b,非a,非b等抽象成点，推出关系抽象成有向边，经典的不成立情况就是a推出非a。如果要赋值的话，就随机找点开始DFS随便赋值

#include<bits/stdc++.h>
using namespace std;
#define int long long
int read()
{
	int x=0,f=1;
	char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
int n,m,tot,head[2000010];
struct edge{int to,next;}e[2000010];
void add_edge(int u,int v)
{
	e[++tot].to=v;
	e[tot].next=head[u];
	head[u]=tot;
}
const int N=2000010;
int dfn[N],st[N],low[N],col[N],ans[N];
int top,scc,dfs_clock;
bool vis[N];
void tarjan(int u)
{
	vis[u]=1;
	st[++top]=u;
	dfn[u]=low[u]=++dfs_clock;
	for(int i=head[u];i;i=e[i].next)
	{
		int v=e[i].to;
		if(!dfn[v])
		{
			tarjan(v);
			low[u]=min(low[u],low[v]);
		}
		else if(vis[v]) low[u]=min(low[u],dfn[v]);
	}
	if(dfn[u]==low[u])
	{
		scc++;int s; 
		do
		{
			s=st[top--];
			vis[s]=0,col[s]=scc;
		}while(s!=u);
	}
}
signed main()
{
	n=read(),m=read();
	while(m--)
	{
		int xi=read(),a=read();
		int xj=read(),b=read();
		add_edge(2*xi+!a,2*xj+b);
		add_edge(2*xj+!b,2*xi+a);
	}
	for(int i=1;i<=2*n;i++) if(!dfn[i]) tarjan(i);
	for(int i=1;i<=n;i++)
	{
		if(col[i*2]==col[i*2+1])
		{
			puts("IMPOSSIBLE");
			return 0;
		}
		else ans[i]=col[i*2]>col[i*2+1];
	}
	puts("POSSIBLE");
	for(int i=1;i<=n;i++) printf("%lld ",ans[i]);
	return 0;
}

二分图最大匹配：

二分图最大匹配_xiyuping24的博客-CSDN博客

LCA：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int N=500002;
int n,m,s,k=0;
int head[N],dep[N],f[N][21];
struct edge{ int v,nxt; }e[N*2];
void add(int u,int v)
{
	e[k].v=v;
	e[k].nxt=head[u];
	head[u]=k++;
}
void dfs(int u,int fa)
{
	dep[u]=dep[fa]+1;
	f[u][0]=fa;
	for(int i=1;(1<<i)<=dep[u];i++)
	{
		f[u][i]=f[f[u][i-1]][i-1];
	}
	for(int i=head[u];i!=-1;i=e[i].nxt)
	{
		if(e[i].v!=fa) dfs(e[i].v,u);
	}
}
int LCA(int a,int b)
{
	if(dep[a]>dep[b]) swap(a,b);
	for(int i=20;i>=0;i--) if(dep[a]<=dep[b]-(1<<i)) b=f[b][i];
	if(a==b) return a;
	for(int i=20;i>=0;i--)
	{
		if(f[a][i]!=f[b][i])
		{
			a=f[a][i];
			b=f[b][i];
		}
	}
	return f[a][0];
}
int main()
{
	memset(head,-1,sizeof(head));
	int a,b;
	scanf("%d%d%d",&n,&m,&s);
	for(int i=1;i<n;i++)
	{
		scanf("%d%d",&a,&b);
		add(a,b);
		add(b,a);
	}
	dfs(s,0);
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&a,&b);
		printf("%d\n",LCA(a,b));
	}
	return 0;
}