n个点m条有向边，求在入度为零的点到n号点的所有路 //径中，哪条边被这些路径覆盖的次数最多

 //n个点m条有向边，求在入度为零的点到n号点的所有路
 //径中，哪条边被这些路径覆盖的次数最多
 //有关DAG的知识，先记个模板
 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 using namespace std;
 typedef long long ll;
 const int maxn=5e3+;
 const int maxm=5e4+;
 struct node
 {
     int v;
     int next;
     ll cnt;
 }edge[maxm],Edge[maxm];
 int head[maxn];
 int HEAD[maxn];
 int cnt;
 ll f[maxn],g[maxn];
 int id[maxn],od[maxn];
 void add(int u,int v)
 {
     ++cnt;
     edge[cnt].v=v;
     edge[cnt].next=head[u];
     head[u]=cnt;
     Edge[cnt].v=u;
     Edge[cnt].next=HEAD[v];
     HEAD[v]=cnt;
 }
 int main()
 {
     int i,j,n,m,v;
     scanf("%d%d",&n,&m);
     int a,b;
     for( i=;i<=m;i++){
         scanf("%d%d",&a,&b);
         add(a,b);
         id[b]++; od[a]++;
     }
     for(int i=;i<=n;i++){
        if(!id[i]) f[i]=;
        if(!od[i]) g[i]=;
     }
     for(i=;i<=n;i++){
        for(j=head[i];j;j=edge[j].next){
            v=edge[j].v;
            f[v]+=f[i];
            edge[j].cnt+=f[i];
        }
     }
     for(j=n;j>=;j--){
        for(i=HEAD[j];i;i=Edge[i].next){
            v=Edge[i].v;
            g[v]+=g[j];
            Edge[i].cnt+=g[j];
        }
     }
     ll ans=;
     for(i=;i<=m;i++)
         ans=max(ans,edge[i].cnt*Edge[i].cnt);
     printf("%lld\n",ans);
     return ;
 }