【EOJ Monthly 2021.2】 B. 杨柳依依

昨天参加了【EOJ Monthly 2021.2】，感受非常的微妙 ，今天回过头再来看这几道题，依然是那么的让人头昏脑涨 [看来还是我太$vegetable$了 ]。
今天重点订了第二道题，来分享一下好了。

参考博客

【EOJ Monthly 2021.2】

题目

赛时感受

最开始做这道题的时候距比赛结束还有不到50分钟（大部分时间都用在了第一题和第三题），读完题后只想着能拿几分是几分，于是看了一眼数据范围：

子任务1 和 子任务3 都只是一条链，很和谐嘛，so就一通操作猛如虎 ，拿了4分 2333 。
再回过头来看，子任务1 和 子任务2 的 k=1；也很和谐嘛，k=1就意味着只要出发的那个点就可以满足题意了，于是又是一波操作猛如虎 ，拿了9分，这时距比赛结束还有20分钟，就去看了看最后一道题，发现并没有什么思路也不能骗到分，so本场考试到此结束，2333 。

算法分析

这道题考查的是最短路问题，但要解最短路中经过结点 $i$ 的所有路径的条数。两遍$BFS$求最短路。第一遍$BFS$找到起点 $s$ 到点 $i$ 的方案数，第二遍$BFS$找到终点 $t$ 到起点 $i$ 的方案数，相乘就是总的可能路线。

具体代码和个人理解见下：

AC Code

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<vector>
using namespace std;
const int MAXN=5e3+100;
const double eps=1e-6;
vector<int> G[MAXN];
int num1[MAXN],num2[MAXN],d[MAXN];
int anspot,n,m,k;
double ans[MAXN];
inline int sread()
{
 int x=0,f=1;char c=getchar();
 while(c>'9'||c<'0') {if(c=='-') f=-1;c=getchar();}
 while(c>='0'&&c<='9') {x=x*10+c-'0';c=getchar();}
 return f*x; 
}
int main()
{
 n=sread();m=sread();
 int x,y;
 for(int i=1;i<=m;++i)
 {
  x=sread();y=sread();
  G[x].push_back(y);
  G[y].push_back(x);
 } k=sread();
 //输入
 for(int i=1;i<=k;++i)
 {
  x=sread(); y=sread(); 
 
  memset(num1,0,sizeof(num1));
  memset(num2,0,sizeof(num2));
  memset(d,0,sizeof(d));
 
  queue<int> q;//建队列 
  d[x]=0;//距离 
  
  num1[x]=1;q.push(x);
  while(!q.empty())//第一遍BFS找到起点s到点i的方案数
  {
   int cur=q.front(); q.pop();//取出队首 
   if(cur==y) break; //到达终点，跳出。 
   for(auto i:G[cur])
//for(auto a:b)中b为一个容器，效果是利用a遍历并获得b容器中的
//每一个值，但是a无法影响到b容器中的元素。 
   {
    if(num1[i]==0)//遍历到i点的方案数为 0 ，即未遍历到过该点，入队。 
    {
     q.push(i);
     d[i]=d[cur]+1;
    }
    if(d[i]==d[cur]+1)//若i点可从cur点一步得到，则 num1[i]+=num1[cur]
    {
     num1[i]+=num1[cur];
    } 
   }
  }
  
  while(!q.empty()) q.pop();//清空队列 
  
  num2[y]=1;q.push(y);
  while(!q.empty())//第二遍BFS找到终点t到起点i的方案数
  {
   int cur=q.front();q.pop();
   if(cur==x) break;
   for(auto i:G[cur]){
    if(d[i]==d[cur]-1){
     if(num2[i]==0) q.push(i);
     num2[i]+=num2[cur];
    }
   }
  }//理解同第一遍BFS 
  
  for(int i=0;i<n;++i) ans[i]+=1.0*num2[i]*num1[i]/num2[x];
 }
 
 for(int i=1;i<n;++i) 
  if(ans[i]-ans[anspot]>=eps) anspot=i;//ans[]最大的点 
 
 printf("%d\n",anspot)；//输出结果

 return 0; //完结撒花
}

反思与总结

1.比赛经验太少，要多做题，多打比赛。
2.本次比赛的四道题都与图论有关，可见在图论方面掌握的还不是很好啊。
3.本题用到了$bfs$，对于搜索要灵活运用，掌握好了会很好用的。
4.定义变量的时候要让自己看得懂，不然就会出现一段时间后再看代码不知道变量的具体含义。
5.在订正题目的过程中看到了好多$dalao$的代码，发现这个世界原来还有好多好多好多需要学习的东西。
革命尚未完成，同志仍需努力啊！