图论专题总结及近年NOIP有关真题

目录

写在前面

1.图的建立和遍历

2.最短路

3.拓扑排序

4.生成树

5.连通块

6.二分图，网络流等

---

写在前面

这篇文章主要总结部分NOIP的图论知识点以及一些NOIP真题。

1.图的建立和遍历

图的存储方式主要就是邻接矩阵和前向星这两种，这里不多赘述。

遍历说白了就是纯粹的暴力。。。。但也是一种很重要的拿部分分的方式。。。

主要就是深度优先搜索和广度优先搜索两种搜索方式

纯考这个的好像就一道信息传递吧。。虽然很水就是了。不过2003年的传染病控制我觉得很迷。。。

信息传递：

直接dfs找最小环即可，连tarjan都不用因为不用存储环

#include<bits/stdc++.h>
using namespace std;
const int MAXN=2e5+10;
const int INF=0x3f3f3f3f;

int n;
int to[MAXN],dep[MAXN],num[MAXN];
int tot,ans=INF;

int Read(){
	int i=0,f=1;
	char c;
	for(c=getchar();(c>'9'||c<'0')&&c!='-';c=getchar());
	if(c=='-')
	  f=-1,c=getchar();
	for(;c>='0'&&c<='9';c=getchar())
	  i=(i<<3)+(i<<1)+c-'0';
	return i*f;
}

void dfs(int u,int d)
{
    dep[u]=d;
    num[u]=tot;
    if(dep[to[u]]){
        if(num[to[u]]<tot){
            return ;
        }
        ans=min(ans,d-dep[to[u]]+1);
    }
    else
      dfs(to[u],d+1);
}

int main()
{
    n=Read();
    for(int i=1;i<=n;++i)
      to[i]=Read();
    for(int i=1;i<=n;++i){
        if(!dep[i]){
            tot++;
            dfs(i,1);
        }
    }
    cout<<ans;
    return 0;
}

2.最短路

求最短路的方法主要就是SPFA和迪杰斯特拉两种，由于前面一种太容易被卡成O(nm)的复杂度，所以一般推荐堆优化的迪杰斯特拉，基本稳定O((n+m)logn)

而需要用Floyd的题不多，需要结合数据范围来看，如果数据范围较小且多次询问两点间最短路，可能需要上Floyd，另外这个也有传递闭包的作用（但CSP选手表示不会）

NOIP近年有关题目：

2017DAY1T3逛公园，2013DAY2T3华容道，2012DAY2T2寻找道路

可以看见一旦考最短路，位置都是在T3位置，但寻找道路这道题真的太水了，让我怀疑这个DAY2T2的水准。

但前面两道都很有意思，第一题侧重于最短路上的DP，第二道则是更侧重于建图，而经典题目如墨墨的等式这种，都是在最短路的基础之上加上数论或者DP，所以刷题时可以侧重这方面。

逛公园：

首先跑出1到所有点的最短路，记作dis[i]。

定义状态方程dp[x][i]表示到点x比最短路大i的路径条数，根据这个数据规模明显直接O(nk)枚举就行。

DP转移方程：，这里因为一个点可以被多次搜到，所以记忆化可以做到剪枝效果。而后面这个式子看着长其实移项后就是，其中表示如果走这条边比最短路长多少，然后用当前枚举长的长度减去这么一坨就是剩下的长度，就是a。

#include<bits/stdc++.h>
using namespace std;
const int MAXN=1e5+10;
const int MAXM=2e5+10;
const int INF=0x3f3f3f3f;
#define pii pair<int,int>
 
int T;
int n,m,p,k,cnt,flag,ans;
int head[MAXN],viss[MAXN],dis[MAXN],dp[MAXN][51],vis[MAXN][51];
int headd[MAXN];
int nxt[MAXM],to[MAXM],w[MAXM];
int nxtt[MAXM],too[MAXM],ww[MAXM];
 
int Read(){
	int i=0,f=1;
	char c=getchar();
	for(;(c>'9'||c<'0')&&c!='-';c=getchar());
	if(c=='-')  f=-1,c=getchar();
	for(;c>='0'&&c<='9';c=getchar())  i=(i<<3)+(i<<1)+c-'0';
	return i*f;
}
 
void add(int x,int y,int z){
	++cnt;
	nxt[cnt]=head[x];
	head[x]=cnt;
	to[cnt]=y;
	w[cnt]=z;
	nxtt[cnt]=headd[y];
	headd[y]=cnt;
	too[cnt]=x;
	ww[cnt]=z;
}
 
void dij(int s){
	memset(dis,INF,sizeof(dis));
	dis[s]=0;
	priority_queue<pii > q;
	q.push(make_pair(0,s));
	while(!q.empty()){
		int u=q.top().second;
		q.pop();
//		viss[u]=0;
		for(int i=head[u];i!=-1;i=nxt[i]){
			int v=to[i];
			if(dis[u]+w[i]<dis[v]){
				dis[v]=dis[u]+w[i];
//				if(!viss[v])
				q.push(make_pair(-dis[v],v));
			}
		}
	}
}
 
int dfs(int u,int k){
//	cout<<u<<' '<<k<<'\n';
	if(~dp[u][k])  return dp[u][k];
	dp[u][k]=0,vis[u][k]=1;
	for(int i=headd[u];i!=-1;i=nxtt[i]){
		int v=too[i];
		int delta=dis[u]+k-dis[v]-ww[i];
		if(delta<0)  continue;
		if(vis[v][delta])  flag=1;
		dp[u][k]=(dp[u][k]+dfs(v,delta))%p;
	}
	vis[u][k]=0;
	return dp[u][k];
}
 
void init(){
	memset(dp,-1,sizeof(dp));
	memset(head,-1,sizeof(head));
	memset(headd,-1,sizeof(headd));
	memset(vis,0,sizeof(vis));
	memset(viss,0,sizeof(viss));
	cnt=0,flag=0;
}