图论板子搞搞

参考博客：
本校潘武杰大佬
本校秦狗
求欧拉回路

#include<cstdio>
const int N=101108,M=201108;
struct Side{
    int v,ne,ok,id;
}S[M<<1];
int t,n,m,sn,cnt,head[N],cur[N],in[N],out[N],fa[N],ans[M<<1];
void init(){
    sn=0;
    for(int i=1;i<=n;i++){
        fa[i]=i;
        head[i]=-1;
        in[i]=out[i]=0;
    }
}
void add(int u,int v,int id){
    S[sn].ok=1;
    S[sn].id=id;
    S[sn].v=v;
    S[sn].ne=head[u];
    head[u]=sn++;
}
int find(int x){
    return fa[x]==x ? x : fa[x]=find(fa[x]); 
}
void bing(int x,int y){
    int fx=find(x),fy=find(y);
    if(fx!=fy) fa[fx]=fy;
    return ;
}
void dfs(int u){
    for(int &i=cur[u];~i;i=S[i].ne){
        if(S[i].ok){
            S[i].ok=0;
            int temp=i;
            if(t==1) S[i^1].ok=0;
            dfs(S[i].v);
            ans[cnt++]=S[temp].id;
            if(i==-1) break;
        }
    }
}
int main(){
    int u,v;
    while(~scanf("%d",&t)){
        scanf("%d%d",&n,&m);
        init();
        for(int i=1;i<=m;i++){
            scanf("%d%d",&u,&v);
            bing(u,v);
            if(t==1){
                in[u]++;
                in[v]++;
                add(u,v,i);
                add(v,u,-i);
            }else{
                in[v]++;
                out[u]++;
                add(u,v,i);
            }
        }
        int beg=1,flag=1,num=0;
        for(int i=1;i<=n;i++){
            cur[i]=head[i];
            if(fa[i]==i&&in[i]) num++;
            if(t==1){
                if(in[i]&1) flag=0;
                else if(in[i]) beg=i;
            }else{
                if(in[i]!=out[i]) flag=0;
                else if(in[i]) beg=i;
            }
            if(!flag||num>1) break;
        }
        if(flag&&num<=1){
            printf("YES\n");
            if(num){
                cnt=0;
                dfs(beg);
                for(int i=cnt-1;i>=0;i--) printf("%d%c",ans[i]," \n"[i==0]);
            }
        }else printf("NO\n");
    }
    return 0;
}

最短哈密顿路径

#include <bits/stdc++.h>
using namespace std;
int n,f[1<<20][21],i,j,k;
int weight[21][21];
int main()
{
    ios::sync_with_stdio(false);//加快cin的读入速度，但是scanf将会不能用。
    memset(f,0x3f,sizeof(f));//初始化最大值
    cin>>n;
    for (i=0; i<n; i++)
        for (j=0; j<n; j++)
            cin>>weight[i][j];
    f[1][0]=0;//第一个点是不需要任何费用的
    for (i=1; i<(1<<n); i++)//i代表着是一个方案集合，其中每一个位置1和0，代表着这个点经过还是没有经过
        for (j=0; j<n; j++)//枚举当前到了哪一个点
            if ((i>>j & 1))//如果i集合中第j位是1，也就是到达过这个点
                for (k=0; k<n; k++)//枚举到达j的点k
                    if ((i^(1<<j)) >> k & 1)//重点，判断k和j的条件，具体在上面解说
                        f[i][j]=min(f[i][j],f[i^(1<<j)][k]+weight[k][j]);//选择最小值，也就是判断，k点到j点最优，还是以前的方案最优
    cout<<f[(1<<n)-1][n-1];//输出最后的最优值
    return 0;
}

有限制条件（走某个点之前必须已经走过某个点）的最短哈密顿路径

#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
using namespace std;
const int INF=3000000;
int d[25][25],dp[2500000][25],f[2500000];
int main(){                                     //就是旅行商问题变形
    int n,m,i,j,u,v,s,ans;                      //dp[s][i]表示状态是s终点是i
    while(scanf("%d%d",&n,&m)!=EOF){
        for(i=0;i<n;i++)
        for(j=0;j<n;j++)
        scanf("%d",&d[i][j]);
        memset(f,0,sizeof(f));
        for(i=0;i<m;i++){                       //将每个点之前要走的点进行
            scanf("%d%d",&u,&v);                //状态压缩
            f[v]|=(1<<u);
        }
        for(i=0;i<(1<<n);i++)
        for(j=0;j<n;j++)
        dp[i][j]=INF;
        dp[1][0]=0;
        for(s=0;s<(1<<n);s++){
            for(i=0;i<n;i++){
                if(dp[s][i]!=INF)
                for(j=0;j<n;j++){
                    if((d[i][j]==-1||(!(s&(1<<i)))||(s&(1<<j))||f[j]!=(s&f[j])))
                    continue;                   //主要判断j点之前的点是否都走过
                    dp[s|(1<<j)][j]=min(dp[s|(1<<j)][j],dp[s][i]+d[i][j]);
                }
            }
        }
        ans=INF;
        for(i=0;i<n;i++)
        ans=min(ans,dp[(1<<n)-1][i]);
        if(ans==INF)
        puts("-1");
        else
        printf("%d\n",ans);
    }
    return 0;
}

LCA倍增

#include<iostream>
#include<algorithm>
#include<string.h>
#include<map>
#include<queue>
#include<cmath>
#include<cstdio>
#include<stack>
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=5e4+5;
const int DEG=22;
struct node{
	int to,next;
	ll val;
}edge[maxn];
int cnt,head[maxn];
ll dis[maxn];
void add(int u,int v,ll val)
{
	edge[cnt].to=v;
	edge[cnt].next=head[u];
	edge[cnt].val=val;
	head[u]=cnt++;
}
void init()
{
	cnt=0;
	memset(head,-1,sizeof head);
	memset(dis,0,sizeof dis);
}
int fa[maxn][DEG];
int deg[maxn];
void BFS(int root)
{
	queue<int> q;
	deg[root]=0;
	fa[root][0]=root;
	q.push(root);
	while(!q.empty())
	{
		int temp=q.front();
		q.pop();
		for(int i=1;i<DEG;i++)
			fa[temp][i]=fa[fa[temp][i-1]][i-1];
		for(int i=head[temp];i!=-1;i=edge[i].next)
		{
			int v=edge[i].to;
			if(v==fa[temp][0])
				continue;
			dis[v]=dis[temp]+edge[i].val;
			deg[v]=deg[temp]+1;
			fa[v][0]=temp;
			q.push(v);
		}
	}
}
int lca(int u,int v)
{
	if(deg[u]>deg[v])
		swap(u,v);
	int hu=deg[u],hv=deg[v];
	int tu=u,tv=v;
	for(int det=hv-hu,i=0;det;det>>=1,i++)
		if(det&1)
			tv=fa[tv][i];
	if(tu==tv)
		return tu;
	for(int i=DEG-1;i>=0;i--)
	{
		if(fa[tu][i]==fa[tv][i])
			continue;
		tu=fa[tu][i];
		tv=fa[tv][i];
	}
	return fa[tu][0];
}
int flag[maxn];
int main()
{
	int t,v,u,n,k;
	ll val;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%d%d",&n,&k);
		init();
		memset(flag,0,sizeof flag);
		for(int i=1;i<n;i++)
		{
			scanf("%d%d%lld",&u,&v,&val);
			add(u,v,val);
			add(v,u,val);
			flag[v]=true;
		}
		int root;
		for(int i=1;i<=n;i++)
			if(!flag[i])
			{
				root=i;
				break;
			}
		BFS(root);
		for(int i=0;i<k;i++)
		{
			scanf("%d%d",&u,&v);
			int temp=lca(u,v);
			printf("%lld\n",dis[u]+dis[v]-2ll*dis[temp]);
		}
	}
	return 0;
}