#include<stdio.h>
#include<algorithm>
#include<vector>
using namespace std;
vector<pair<int,int> >g[110000];
bool th=true;
int v[110000];
int num[110000];
int cnt[3];
int tc;
void dfs(int a,int b){
 
	num[a]=b;
	cnt[b]++;
	for(int i=0;i<g[a].size();i++){
		if(g[a][i].second==1)tc++;
		if(!v[g[a][i].first]){
			v[g[a][i].first]=1;
			dfs(g[a][i].first,(b+g[a][i].second)%3);
		}else{
			if(num[g[a][i].first]!=(b+g[a][i].second)%3)th=false;
		}
	}
}
 
int main(){
	int a,b;scanf("%d%d",&a,&b);
	for(int i=0;i<b;i++){
		int p,q;scanf("%d%d",&p,&q);
		p--;q--;
		g[p].push_back(make_pair(q,1));
		g[q].push_back(make_pair(p,2));
	}
	long long ret=0;
	for(int i=0;i<a;i++){
		if(!v[i]){
			for(int j=0;j<3;j++)cnt[j]=0;
			th=true;
			v[i]=1;
			tc=0;
			dfs(i,0);
			if(th){
				if(cnt[0]&&cnt[1]&&cnt[2])ret+=(long long)cnt[0]*cnt[1]+(long long)cnt[1]*cnt[2]+(long long)cnt[2]*cnt[0];
				else ret+=tc;
			}else{
				int tot=cnt[0]+cnt[1]+cnt[2];
				ret+=(long long)tot*tot;
			}
		}
	}
	printf("%lld\n",ret);
}