// Copyright (C) 2017 __debug.

// This program is free software; you can redistribute it and/or
// modify it under the terms of the GNU General Public License as
// published by the Free Software Foundation; version 3

// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU General Public License for more details.

// You should have received a copy of the GNU General Public License
// along with this program; If not, see <http://www.gnu.org/licenses/>.


#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/priority_queue.hpp>

#define x first
#define y second
#define MP std::make_pair
#define SZ(x) ((int)(x).size())
#define ALL(x) (x).begin(), (x).end()
#define DEBUG(...) fprintf(stderr, __VA_ARGS__)
#ifdef __linux__
#define getchar getchar_unlocked
#define putchar putchar_unlocked
#endif

using std::pair;
using std::vector;
using std::string;

typedef long long LL;
typedef pair<int, int> Pii;

const int oo = 0x3f3f3f3f;

template<typename T> inline bool chkmax(T &a, T b) { return a < b ? a = b, true : false; }
template<typename T> inline bool chkmin(T &a, T b) { return b < a ? a = b, true : false; }
string procStatus()
{
    std::ifstream t("/proc/self/status");
    return string(std::istreambuf_iterator<char>(t), std::istreambuf_iterator<char>());
}
template<typename T> T read(T &x)
{
    int f = 1;
    char ch = getchar();
    for (; !isdigit(ch); ch = getchar())
        f = (ch == '-' ? -1 : 1);
    for (x = 0; isdigit(ch); ch = getchar())
        x = 10 * x + ch - '0';
    return x *= f;
}
template<typename T> void write(T x)
{
    if (x == 0) {
        putchar('0');
        return;
    }
    if (x < 0) {
        putchar('-');
        x = -x;
    }
    static char s[20];
    int top = 0;
    for (; x; x /= 10)
        s[++top] = x % 10 + '0';
    while (top)
        putchar(s[top--]);
}
// EOT

const int MAXN = 1e5 + 5;

int N, M;
vector<int> G[MAXN], Gr[MAXN];

void input()
{
    read(N); read(M);
    for (int i = 1; i <= M; ++i) {
        int u, v;
        read(u); read(v);
        G[u].push_back(v);
        Gr[v].push_back(u);
    }
}

int color[MAXN];
int cnt[3], cntv, cnte;
int flag;

void dfs(int u, int c)
{
    ++cnt[color[u] = c];
    ++cntv;
    cnte += SZ(G[u]);
    int nc = (color[u] + 1) % 3, pc = (color[u] + 2) % 3;
    for (int i = 0; i < SZ(G[u]); ++i) {
        int v = G[u][i];
        if (!~color[v])
            dfs(v, nc);
        else if (color[v] != nc)
            flag = false;
    }
    for (int i = 0; i < SZ(Gr[u]); ++i) {
        int v = Gr[u][i];
        if (!~color[v])
            dfs(v, pc);
        else if (color[v] != pc)
            flag = false;
    }
}

void solve()
{
    memset(color, -1, sizeof(color));
    LL ans = 0;
    for (int i = 1; i <= N; ++i) {
        if (~color[i])
            continue;
        cnt[0] = cnt[1] = cnt[2] = 0;
        cntv = cnte = 0;
        flag = true;
        dfs(i, 0);
        if (!flag)
            ans += (LL)cntv * cntv;
        else if (cnt[0] && cnt[1] && cnt[2])
            ans += (LL)cnt[0] * cnt[1] + (LL)cnt[1] * cnt[2] + (LL)cnt[2] * cnt[0];
        else
            ans += cnte;
    }

    printf("%lld\n", ans);
}

int main()
{
#ifdef __DEBUG
    freopen("F.in", "r", stdin);
    freopen("F.out", "w", stdout);
#endif

    input();
    solve();

    return 0;
}

// 章台柳，章台柳，昔日青青今在否？
//     -- 韩翃《章台柳·寄柳氏》


Measured memory_bytes: 19456