#include<bits/stdc++.h>
#define maxn 100009
#define LL long long
using namespace std;
struct data{
  int nex,to,w;
}e[maxn*2];
int head[maxn],edge=0,tot,t[3],co[maxn],flag;
inline void add(int from,int to,int w){
  e[++edge]=(data){head[from],to,w};
  head[from]=edge;
}
void dfs(int x){
  t[co[x]]++;
  for(int i=head[x];i;i=e[i].nex){
    int v=e[i].to;if(e[i].w==1)tot++;
    if(co[v]==-1) co[v]=(co[x]+e[i].w+3)%3,dfs(v);
    else if(co[v]!=(co[x]+e[i].w+3)%3) flag=1;
  }
}
int main(){
  int n,m;LL ans=0;
  scanf("%d%d",&n,&m);
  for(int i=1,x,y;i<=m;i++)
    scanf("%d%d",&x,&y),add(x,y,1),add(y,x,-1);
  memset(co,-1,sizeof(co));
  for(int i=1;i<=n;i++)
    if(co[i]==-1){
      t[0]=t[1]=t[2]=tot=flag=0;
      co[i]=0;dfs(i);
      if(flag) ans+=(LL)(t[1]+t[2]+t[0])*(t[1]+t[2]+t[0]);
      else if(t[0] && t[1] && t[2]) ans+=(LL)t[0]*t[1]+(LL)t[1]*t[2]+(LL)t[2]*t[0];
      else ans+=tot;
    }
  printf("%lld",ans);
  return 0;
}
