#include <bits/stdc++.h>
using namespace std;

#define _p(...) (void)printf(__VA_ARGS__)
#define forr(x,arr) for(auto&& x:arr)
#define _overload3(_1,_2,_3,name,...) name
#define _rep2(i,n) _rep3(i,0,n)
#define _rep3(i,a,b) for(int i=int(a);i<int(b);++i)
#define rep(...) _overload3(__VA_ARGS__,_rep3,_rep2,)(__VA_ARGS__)
#define _rrep2(i,n) _rrep3(i,0,n)
#define _rrep3(i,a,b) for(int i=int(b)-1;i>=int(a);i--)
#define rrep(...) _overload3(__VA_ARGS__,_rrep3,_rrep2,)(__VA_ARGS__)
#define all(x) (x).begin(),(x).end()
#define bit(n) (1LL<<(n))
#define sz(x) ((int)(x).size())
#define ten(n) ((int)1e##n)
#define fst first
#define snd second
using ll=long long;
using pii=pair<int,int>;using pll=pair<ll,ll>;using pil=pair<int,ll>;using pli=pair<ll,int>;
using vs=vector<string>;using vvs=vector<vs>;using vvvs=vector<vvs>;
using vb=vector<bool>;using vvb=vector<vb>;using vvvb=vector<vvb>;
using vi=vector<int>;using vvi=vector<vi>;using vvvi=vector<vvi>;
using vl=vector<ll>;using vvl=vector<vl>;using vvvl=vector<vvl>;
using vd=vector<double>;using vvd=vector<vd>;using vvvd=vector<vvd>;
using vpii=vector<pii>;using vvpii=vector<vpii>;using vvvpii=vector<vvpii>;
template<class T> bool amax(T &a, const T &b) { if (a < b) { a = b; return 1; } return 0; }
template<class T> bool amin(T &a, const T &b) { if (b < a) { a = b; return 1; } return 0; }
ll ri(){ll l;cin>>l;return l;} string rs(){string s;cin>>s;return s;}
template<class T>T read(){T t;cin>>t;return t;}
template<class T,class U>ostream&operator<<(ostream&o,const pair<T,U>&p){return o<<'('<<p.fst<<", "<<p.snd<<')';}
template<class T>ostream&operator<<(ostream&o,const vector<T>&t){o<<"{";rep(i,sz(t))o<<t[i]<<",";o<<"}"<<endl;return o;}
#ifdef LOCAL
vs s_p_l_i_t(const string&s,char c){vs v;stringstream ss(s);string x;while(getline(ss,x,c))v.emplace_back(x);return move(v);}
void e_r_r(vs::iterator it){cerr<<endl;}
template<class T,class... Args>void e_r_r(vs::iterator it,T a,Args... args){cerr<<it->substr((*it)[0]==' ',it->length())<<" = "<<a<<", ";e_r_r(++it,args...);}
#define out(args...){vs a_r_g_s=s_p_l_i_t(#args,',');e_r_r(a_r_g_s.begin(),args);}
#else
#define out(args...)
#endif

void Main() {
  ll n = ri();
  int m = ri();

  vvi E(n);
  rep(i, m) {
    int u = read<int>() - 1;
    int v = read<int>() - 1;
    E[u].push_back(v);
    E[v].push_back(u);
  }

  vi D(n, -1);
  vb loop(n);

  function<int(int, int, int, int)> dfs = [&](int cur, int pre, int depth, int idx) {
    int ret = 1;
    forr(nex, E[cur]) if (nex != pre) {
      if (D[nex] != -1) {
        if (D[nex] % 2 != (D[cur] + 1) % 2) {
          loop[idx] = 1;
        }
      }
      else {
        D[nex] = D[cur] + 1;
        ret += dfs(nex, cur, depth+1, idx);
      }
    }
    return ret;
  };

  vpii sizes;

  rep(i, n) if (D[i] == -1) {
    D[i] = 0;
    int sz = dfs(i, -1, 1, i);
    sizes.emplace_back(sz, i);
  }

  //out(loop);
  //out(sizes);

  ll x = 0;
  ll o = 0;
  ll s = 0;
  forr(si, sizes) {
    if (si.fst == 1) s++;
    else if (loop[si.snd]) o++;
    else x++;
  }

  //out(o, x);

  ll ans = 1LL * n * n;

  forr(si, sizes) {
    ll size = si.fst;
    if (size == 1) continue;
    ll idx = si.snd;

    //out(si, ans);
    ans -= size * (n-s);
    if (loop[idx]) {
      ans += o + x;
    }
    else {
      ans += x * 2;
      ans += o * 1;
    }
    //out(si, ans);
  }

  cout << ans << endl;
}
int main() { cin.tie(nullptr); ios::sync_with_stdio(false); Main(); return 0; }


Measured memory_bytes: 13696