#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;
}