#include <cstdio>

const int N=100005;

int fst[N],a[N],cnt[3],tot=0,num;
bool flag;

struct Edge{
	int to,v,nxt;
}e[N*2];

void addedge(int x,int y){
	e[++tot].to=y;e[tot].v=1;e[tot].nxt=fst[x];fst[x]=tot;
	e[++tot].to=x;e[tot].v=-1;e[tot].nxt=fst[y];fst[y]=tot;
}

void dfs(int x){
	cnt[a[x]]++;
	for (int i=fst[x];i;i=e[i].nxt){
		int y=e[i].to,z=(a[x]+e[i].v+3)%3;
		num++;
		if (a[y]==-1){
			a[y]=z;
			dfs(y);
		}
		else if (a[y]!=z) flag=1;
	}
}

int main(){
	int n,m;
	scanf("%d%d\n",&n,&m);
	for (int i=1;i<=m;i++){
		int x,y;
		scanf("%d%d\n",&x,&y);
		addedge(x,y);
	}
	for (int i=1;i<=n;i++) a[i]=-1;
	long long ans=0;
	for (int i=1;i<=n;i++)
	if (a[i]==-1){
		for (int j=0;j<3;j++) cnt[j]=0;
		flag=0;a[i]=num=0;dfs(i);
		if (flag) ans+=1ll*(cnt[0]+cnt[1]+cnt[2])*(cnt[0]+cnt[1]+cnt[2]);
		else if (!cnt[0] || !cnt[1] || !cnt[2]) ans+=num/2;
		else for (int i=0;i<3;i++) ans+=1ll*cnt[i]*cnt[(i+1)%3];
	}
	printf("%lld\n",ans);
}