// Slower C++ solution:
#define _CRT_SECURE_NO_WARNINGS

#include <iostream>

#include <sstream>

#include <iomanip>

#include <cstdio>

#include <cstdlib>

#include <cmath>

#include <cctype>

#include <cstring>

#include <vector>

#include <list>

#include <queue>

#include <deque>

#include <stack>

#include <map>

#include <set>

#include <algorithm>

#include <iterator>

#include <bitset>

#include <ctime>

#include<complex>

using namespace std;



#define FOR(i,a,b) for (int i = (a); i < (b); i++)

#define RFOR(i,b,a) for (int i = (b)-1; i >= (a); i--)

#define ITER(it,a) for (__typeof(a.begin()) it = a.begin(); it != a.end(); it++)

#define FILL(a,value) memset(a, value, sizeof(a))



#define SZ(a) (int)a.size()

#define ALL(a) a.begin(), a.end()

#define PB push_back

#define MP make_pair



typedef long long LL;

typedef vector<int> VI;

typedef pair<int, int> PII;



const double PI = acos(-1.0);

const int INF = 1000 * 1000 * 1000 + 7;

const LL LINF = INF * (LL)INF;



const int MAX = 100000 + 47;

const int MOD = 1000 * 1000 * 1000 + 7;



vector<int> g[MAX];

int C[MAX];

int numb;

bool is;

LL o, e, iz;



void dfs(int v, int c)

{

	if (C[v] != -1)

	{

		if (C[v] != c) is = false;

		return;

	}



	numb++;

	C[v] = c;

	FOR(i, 0, SZ(g[v]))

	{

		int to = g[v][i];

		dfs(to, c ^ 1);

	}

}



int main()

{

	//freopen("in.txt", "r", stdin);

	//freopen("out.txt", "w", stdout);

	ios::sync_with_stdio(false); cin.tie(0);

	int n, m;

	cin >> n >> m;

	FOR(i, 0, m)

	{

		int a, b;

		cin >> a >> b;

		a--;

		b--;

		g[a].push_back(b);

		g[b].push_back(a);

	}

	

	FILL(C, -1);

	FOR(i, 0, n)

	{

		if (C[i] == -1)

		{

			numb = 0;

			is = true;

			dfs(i, 0);



			if (numb == 1)

			{

				iz++;

				continue;

			}



			if (is) e++; else o++;

		}

	}



	//cout << res << endl;

	LL res = 2 * e * e + o * o + 2 * o * e + n * iz + iz * (n - iz);

	cout << res << endl;

}

// Faster C++ solution for the same problem (measured speedup: 3.8818319401847963x):
#include "bits/stdc++.h"

using namespace std;

#ifndef LOCAL

#define endl '\n'

#endif



#define fr(i, a, b) for(int i = a; i <= b; i++)

#define rf(i, a, b) for(int i = a; i >= b; i--)

#define pf push_front

#define pb push_back

#define eb emplace_back

#define fi first

#define se second

#define all(x) x.begin(), x.end()

#define rall(x) x.rbegin(), x.rend()

#define sz(x) (int)x.size()

#define lb lower_bound

#define ub upper_bound



typedef long long ll;

typedef long double f80;

typedef pair<int,int> pii;

typedef pair<ll,ll> pll;



int pct(int x) { return __builtin_popcount(x); }

int pct(ll x) { return __builtin_popcountll(x); }

int bt(int x) { return 31 - __builtin_clz(x); } // floor(log2(x))

int bt(ll x) { return 63 - __builtin_clzll(x); } // floor(log2(x))

int cdiv(int a, int b) { return a / b + !(a < 0 || a % b == 0); }

ll cdiv(ll a, ll b) { return a / b + !(a < 0 || a % b == 0); }

int nxt_C(int x) { int c = x & -x, r = x + c; return (((r ^ x) >> 2) / c) | r; }

ll nxt_C(ll x) { ll c = x & -x, r = x + c; return (((r ^ x) >> 2) / c) | r; }



vector<int> get_bits(int mask) {

	vector<int> bb;

	while(mask) { int b = bt(mask); bb.pb(b); mask ^= (1 << b); }

	reverse(all(bb));

	return bb;

}



int get_mask(vector<int> v) {

	int mask = 0;

	for(int x : v) { mask ^= (1 << x); }

	return mask;

}



template<typename T>

void uniq(vector<T> &v) { sort(all(v)); v.resize(unique(all(v)) - v.begin()); }



mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count());



ll rand(ll l, ll r){

	uniform_int_distribution<ll> uid(l, r);

	return uid(rng);

}



void sc() {}



template <typename Head, typename... Tail>

void sc(Head &H, Tail &... T) { cin >> H; sc(T...); }



#ifdef LOCAL

#define debug(...) cerr << "[L:" << __LINE__ << "][" << #__VA_ARGS__ << "]:", debug_out(__VA_ARGS__)

#else

#define debug(...) 42

#endif



// #ifndef LOCAL

// string to_string(__int128 x) {

// 	string s = "";

// 	bool neg = 0;

// 	if(x < 0) { s += "-"; neg = 1; x = -x; }

// 	if(!x) s += '0';

// 	while(x) {

// 		int rem = x % 10;

// 		s += to_string(rem);

// 		x /= 10;

// 	}

// 	reverse(s.begin() + neg, s.end());

// 	return s;

// }

// #endif



const int mod = 1e9 + 7; // 998244353;



int pwr(int a,ll b) {

	int ans = 1;

	while(b) {

		if(b & 1) ans = (ans * 1LL * a) % mod;

		a = (a * 1LL * a) % mod;

		b >>= 1;

	}

	return ans;

}



/*

	Lookout for overflows!!

	Check array sizes!!

	Clear before test cases!!

	Use the correct modulo!!

	Check for corner cases!!

	Are you forgetting something?!

	Read problem statement carefully!!!

*/



const int N = 1e5 + 5;

int dsu[N], sz[N];

bool odd[N], par[N];



int root(int u) {

    while(u != dsu[u]) u = dsu[u];

    return u;

}



int get_par(int u) {

    int p = 0;

    while(u != dsu[u]) {

        p ^= par[u];

        u = dsu[u];

    }

    return p;

}



void uni(int u,int v) {

    int uu = root(u), vv = root(v);

    if(uu == vv) {

        if(get_par(u) == get_par(v))

            odd[uu] = 1;

        return;

    }

    if(sz[uu] > sz[vv]) {

        swap(uu, vv), swap(u, v);

    }

    dsu[uu] = vv;

    sz[vv] += sz[uu];

    odd[vv] |= odd[uu];

    par[uu] = get_par(u) ^ get_par(v) ^ 1;

}



void solve() {

	int n, m;

    sc(n, m);

    fr(i, 1, n) {

        sz[i] = 1;

        dsu[i] = i;

        odd[i] = 0;

        par[i] = 0;

    }

    fr(i, 1, m) {

        int u, v;

        sc(u, v);

        uni(u, v);

    }

    int c1 = 0, c2 = 0, c3 = 0;

    fr(i, 1, n) {

        if(i == dsu[i]) {

            if(sz[i] == 1) {

                c1++;

            }

            else if(odd[i]) {

                c2++;

            }

            else {

                c3++;

            }

        }

    }

    ll ans = 0;

    while(c1) {

        ans += 2 * n - 1;

        n--;

        c1--;

    }

    while(c2) {

        ans += 2 * (c2 + c3) - 1;

        c2--;

    }

    ans += (c3 * 1LL * (c3 - 1)) * 2 + 2LL * c3;

    cout << ans;

}



int main() {

	ios :: sync_with_stdio(0);

	cin.tie(0);

	cout.tie(0);

	int t = 1;

	// cin >> t;

	for(int tt = 1; tt <= t; tt++) {

		solve();

	}

	return 0;

}