題意

給你一個森林(若干棵樹),向其中加邊，使得最後形成一棵樹，要求最後形成的樹的直徑最小，輸出這棵樹的直徑和所加的邊

分析

每加一條邊，都可以將兩棵樹合併成一棵樹，總共有t = n - m棵樹

也就是說我們需要在這t棵樹裡面，分別選一個點，將他們連線起來，使得最後形成一棵直徑最小的樹

這個點的性質是，以這個點為根節點轉化為有根樹的畫，樹的高度一定是最小的（與重心類似，最長鏈最短可以保證合併後直徑最短）

稍加分析可以知道，這個點一定是直徑的中點（直徑兩端到該點的距離相差不超過1）

所以我們維護一棵樹，可以維護以下幾點性質：

1. 樹的直徑dig

2. 樹的直徑的中點序號center

3. 樹的直徑兩端離center較遠的點(相同則任意)site

這樣合併A,B兩顆樹後，新的樹

dig = max ( A.dig , B.dig , ceil(A.dig/2) + ceil(B.dig/2) )

遍歷這棵樹，仍然求得新的center和site

那麼以怎樣的順序合併呢？

（以下未嚴格證明，只是YY）

我們想要直徑儘量小，那麼就是儘量向直徑最大的子樹連邊

這樣，如果另一棵樹直徑小於這棵樹，則直徑不變，若等於，形成一棵新的直徑大1的子樹

所以我們希望最小化這個影響，所以考慮每次將最大的兩顆樹合併起來，直到合成一棵樹

這樣可以用優先佇列來維護這個過程

程式碼

  1
 
 #include<bits/stdc++.h>
  2 using namespace std;
  3 
  4 int n,m;
  5 vector<int>edge[1005];
  6 int color[1005];
  7 void rs(int pre,int now,int co) {
  8     color[now] = co;
  9     for(auto it = edge[now].begin();it!=edge[now].end();it++) {
 10     if(*it == pre) continue;
 11     rs(now,*it,co);
 
 12     }
 13     return;
 14 }
 15 struct Tree{
 16     int dig,center,site;
 17     bool operator < (const Tree& o) const {
 18     return dig < o.dig;
 19     }
 20 }tree;
 21 priority_queue<Tree>pq;
 22 int far,X,Y,id;
 23 void dfs(int pre,int now,int dis){
 24     if(dis>far) {
 25     far = dis;id =now;
 26     }
 27     for(auto it = edge[now].begin();it!=edge[now].end();it++) {
 28     if(*it == pre) continue;
 29     dfs(now,*it,dis+1);
 30     }
 31 }
 32 stack<int>st;
 33 int flag;
 34 vector<pair<int,int> >ans;
 35 void dfs2(int pre,int now,int aim){
 36     if(now == aim) {
 37     tree.site = now;
 38     for(int i=1;i<(far+1)/2;i++){ st.pop();}
 39     tree.center = st.top();
 40     tree.dig = far;
 41     pq.push(tree);
 42     flag = 1;
 43     while(!st.empty())st.pop();
 44     return;
 45     }
 46     st.push(now);
 47     for(auto it = edge[now].begin();it!=edge[now].end();it++) {
 48     if(*it == pre) continue;
 49     dfs2(now,*it,aim);
 50     if(flag) return;
 51     }
 52     st.pop();
 53 }
 54 int main() {
 55     ios::sync_with_stdio(false);
 56     cin>>n>>m;
 57     int x,y;
 58     for(int i=1;i<=m;i++) {
 59     cin>>x>>y;
 60     edge[x].push_back(y);
 61     edge[y].push_back(x);
 62     }
 63     int co = 0;
 64     for(int i=1;i<=n;i++) {
 65     if(color[i]==0) {
 66         rs(0,i,++co);
 67         far=0;id = i;
 68         dfs(0,i,0);
 69         X = id;
 70         far=0;id = X;
 71         dfs(0,X,0);
 72         Y = id;
 73         if(X == Y) {
 74         tree.site = X;
 75         tree.center = X;
 76         tree.dig = 0;
 77         pq.push(tree);
 78         }
 79         else {
 80         flag =0;
 81         dfs2(0,X,Y);
 82         }
 83     }
 84     }
 85     while(pq.size()>1) {
 86     Tree a = pq.top();pq.pop();
 87     Tree b = pq.top();pq.pop();
 88     far = max(a.dig,((a.dig+1)/2) + ((b.dig+1)/2) +1);
 89     far = max(b.dig,far);
 90     edge[a.center].push_back(b.center);
 91     edge[b.center].push_back(a.center);
 92     flag =  0;
 93     dfs2(0,a.site,b.site);
 94     ans.push_back(make_pair(a.center,b.center));
 95     }
 96     tree = pq.top();
 97     cout<<tree.dig<<endl;
 98     
 99     for(auto it = ans.begin();it!=ans.end();it++) {
100     cout<<it->first<<" "<<it->second<<endl;
101     }
102     
103 
104 }