#include <iostream>
#include <cstdio>
#include <cstring>
#define N 100005
using namespace std;
typedef long long ll;
struct edge{
	int k,next;
}e[N<<2];
int n,m,home[N],cnt=-1,H[N];
void add(int *home,int x,int y){
	cnt++;
	e[cnt].k=y;
	e[cnt].next=home[x];
	home[x]=cnt;
}
int col[N],D,sum[5],P,E;
bool vis[N];
void dfs(int k){
	sum[col[k]]++;
	int w=col[k]%3+1;
	vis[k]=1;P++;
	for(int i=home[k];~i;i=e[i].next){
		E++;
		if(col[e[i].k]&&col[e[i].k]!=w){
			D=0;
		}else{
			col[e[i].k]=w;
		}
	}
	int W=col[k]-1;
	if(W==0) W=3;
	for(int i=H[k];~i;i=e[i].next) if(!vis[e[i].k]){
		if(!col[e[i].k]) col[e[i].k]=W;
		dfs(e[i].k);
	}
}
int main(){
	memset(home,-1,sizeof(home));
	memset(H,-1,sizeof(H));
	scanf("%d%d",&n,&m);
	for(int i=1,x,y;i<=m;i++){
		scanf("%d%d",&x,&y);
		add(home,x,y);
		add(H,x,y);add(H,y,x);
	}
	ll ans=0;
	for(int i=1;i<=n;i++){
		if(!col[i]){
			sum[1]=sum[2]=sum[3]=E=P=0;
			D=1,col[i]=1,dfs(i);
			if(!D) ans+=1LL*P*(P+1)-P;
			else if(!sum[3]||!sum[2]||!sum[1]) ans+=E;
			else ans+=1LL*sum[1]*sum[2]+1LL*sum[2]*sum[3]+1LL*sum[3]*sum[1];
		}
	}
	printf("%lld",ans);
	return 0;
}