#include<iostream>
#include<string>
#include<cstdio>
#include<vector>
#include<cmath>
#include<algorithm>
#include<functional>
#include<iomanip>
#include<queue>
#include<ciso646>
#include<random>
#include<map>
#include<set>
#include<complex>
#include<bitset>
#include<stack>
#include<unordered_map>
#include<utility>
#include<tuple>
#include<cassert>
using namespace std;
typedef long long ll;
typedef unsigned int ui;
const ll mod = 1000000007;
const ll INF = (ll)1000000007 * 1000000007;
typedef pair<int, int> P;
#define stop char nyaa;cin>>nyaa;
#define rep(i,n) for(int i=0;i<n;i++)
#define per(i,n) for(int i=n-1;i>=0;i--)
#define Rep(i,sta,n) for(int i=sta;i<n;i++)
#define Per(i,sta,n) for(int i=n-1;i>=sta;i--)
#define rep1(i,n) for(int i=1;i<=n;i++)
#define per1(i,n) for(int i=n;i>=1;i--)
#define Rep1(i,sta,n) for(int i=sta;i<=n;i++)
typedef long double ld;
const ld eps = 1e-8;
const ld pi = acos(-1.0);
typedef pair<ll, ll> LP;
int dx[4]={1,-1,0,0};
int dy[4]={0,0,1,-1};

struct SCC{
  vector<vector<int> > G,R,T,C;//T:強連結分解後のグラフ
  vector<int> vs,used,blg;
  SCC(int n):G(n),R(n),used(n),blg(n){}

  void add_edge(int u,int v){
    G[u].emplace_back(v);
    R[v].emplace_back(u);
  }

  int size(){
    return T.size();
  }

  void dfs(int v){
    used[v]=1;
    for(int u:G[v])
      if(!used[u]) dfs(u);
    vs.emplace_back(v);
  }

  void rdfs(int v,int k){
    used[v]=1;
    blg[v]=k;
    C[k].emplace_back(v);
    for(int u:R[v])
      if(!used[u]) rdfs(u,k);
  }

  int build(){
    int n=G.size();
    for(int v=0;v<n;v++)
      if(!used[v]) dfs(v);

    fill(used.begin(),used.end(),0);
    int k=0;
    for(int i=n-1;i>=0;i--){
      if(!used[vs[i]]){
        T.emplace_back();
        C.emplace_back();
        rdfs(vs[i],k++);
      }
    }
    for(int v=0;v<n;v++)
      for(int u:G[v])
        if(blg[v]!=blg[u])
          T[blg[v]].push_back(blg[u]);

    for(int i=0;i<k;i++){
      sort(T[i].begin(),T[i].end());
      T[i].erase(unique(T[i].begin(),T[i].end()),T[i].end());
    }
    return k;
  }
  int operator[](int k) const{return blg[k];};//頂点kの属する強連結成分
};

ll n;int m;
int num[200010];ll sco[200010];
ll v[200010];
vector<int> zaatu;
P e[200010];
ll ans=0;

void dfs(int s, vector<vector<int>> &G){
    if(v[s]!=0) return;
    v[s]=max(v[s],sco[s]);
    for(int t:G[s]){
        if(v[t]==0) dfs(t,G);
        v[s]=max(v[t],v[s]);
    }
    ans+=v[s]*(ll)num[s];
}

void solve(){
    cin >> n >> m;
    rep(i,m){
        int b,c;cin >> b >> c;
        zaatu.push_back(b);
        zaatu.push_back(c);
        e[i]=P(b,c);
    }
    sort(zaatu.begin(),zaatu.end());
    zaatu.erase(unique(zaatu.begin(),zaatu.end()),zaatu.end());
    int V=zaatu.size();
    SCC scc(V);
    // rep(i,V){
    //     cout << zaatu[i] << endl;
    // }
    rep(i,m){
        int u=lower_bound(zaatu.begin(),zaatu.end(),e[i].first)-zaatu.begin();
        int v=lower_bound(zaatu.begin(),zaatu.end(),e[i].second)-zaatu.begin();
        scc.add_edge(u,v);
    }
    int k=scc.build();
    ll d=0;
    rep(i,V){
        sco[scc[i]]=max(sco[scc[i]],(ll)zaatu[i]);
        num[scc[i]]+=1;
        d+=zaatu[i];
    }
    // rep(i,k){
    //     cout << i << " " << sco[i] << " " << num[i] << endl;
    // }
    vector<vector<int>> G=scc.T;
    rep(i,k){
        dfs(i,G);
    }
    cout << n*(n+1)/2+ans-d << endl;
}

int main(){
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout << fixed << setprecision(50);
    solve();
}