CF14D Two Paths

• 给一棵树，保证联通，求两条没有重叠部分的树链长度乘积的最大值。n<=200
• 两条没有重叠部分的树链，一定可以存在于两颗没有重复的子树，进一步转换一下，段掉一条边，把原来的树分成两个子树，求两棵树树的直径乘积，枚举断掉的这一条边即可，求出最大乘积。

#include <iostream>
#include <cstdio>
using namespace std;
int n,last[4005],cnt,x,y,maxlen,node,ans,maxn;
bool mj[4005][4005],inq[4005];
struct edge{
    int v,next;
}e[4005];
inline void add(int u,int v)
{
    e[++cnt].v=v;
    e[cnt].next=last[u];
    last[u]=cnt;
    e[++cnt].v=u;
    e[cnt].next=last[v];
    last[v]=cnt;
}
void dfs(int u,int len)
{
    inq[u]=1;
    for(int i=last[u];i;i=e[i].next)
    {
        if(mj[u][e[i].v]||inq[e[i].v])
        continue;
        dfs(e[i].v,len+1);
    }
    inq[u]=0;
    if(len>maxlen)
    {
        maxlen=len;
        node=u;
    }
}
int main()
{
    cin>>n;
    for(int i=1;i<=n-1;i++)
    {
        cin>>x>>y;
        add(x,y);
    }
    for(int i=1;i<=2*n-3;i+=2)
    {
        ans=maxlen=-1;
        mj[e[i].v][e[i+1].v]=mj[e[i+1].v][e[i].v]=1;
        dfs(e[i].v,0);
        maxlen=-1;
        dfs(node,0);
        ans=maxlen;
        maxlen=-1;
        dfs(e[i+1].v,0);
        maxlen=-1;
        dfs(node,0);
        ans*=maxlen;
        maxn=max(maxn,ans);
        mj[e[i].v][e[i+1].v]=mj[e[i+1].v][e[i].v]=0;
    }
    cout<<maxn;
    return 0;
}