// Slower C++ solution:
#include<cstdio>

#include<cstring>

#include<vector>

#include<algorithm>

using namespace std;

const long long MAXN=100005;

long long n,m;

vector<int>w[MAXN];

bool vis[MAXN];

long long cl[MAXN];

bool dfs(long long u)

{

	vis[u]=1;

	long long re=1;

	for(long long i=0;i<int(w[u].size());i++)

	{

		long long v=w[u][i];

		if(vis[v])

			re&=(cl[u]!=cl[v]);

		else

		{

			cl[v]=cl[u]^1;

			re&=dfs(v);

		}

	}

	return re;

}

int main()

{

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

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

	scanf("%lld %lld",&n,&m);

	for(long long i=1;i<=m;i++)

	{

		long long u,v;

		scanf("%lld %lld",&u,&v);

		w[u].push_back(v);

		w[v].push_back(u);

	}

	long long I=0,O=0,J=0;

	for(long long i=1;i<=n;i++)

	{

		if(!vis[i])

		{

			if(int(w[i].size())==0)

				I++;

			else if(dfs(i))

				O++;

			else

				J++;

		}

	}

	printf("%lld\n",I*I+I*(n-I)*2+O*O*2+J*J+O*J*2);

	fclose(stdin);

	fclose(stdout);

}

// Faster C++ solution for the same problem (measured speedup: 3.380958493448055x):
#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;

}