//include
//------------------------------------------
#include <vector>
#include <list>
#include <map>
#include <set>
#include <deque>
#include <queue>
#include <stack>
#include <bitset>
#include <algorithm>
#include <functional>
#include <numeric>
#include <utility>

#include <sstream>
#include <iostream>
#include <iomanip>
#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <cctype>
#include <string>
#include <cstring>
#include <ctime>

#include <stdexcept>

using namespace std;

//conversion
//------------------------------------------
inline int toInt(string s) { int v; istringstream sin(s); sin >> v; return v; }
template<class T> inline string toString(T x) { ostringstream sout; sout << x; return sout.str(); }

//debug
#define DUMP(x)  cerr << #x << " = " << (x) << endl
#define DEBUG(x) cerr << #x << " = " << (x) << " (L" << __LINE__ << ")" << " " << __FILE__ << endl

// type alias
using ll = long long;
using ull = unsigned long long;

// constant value
const ll MOD = 1000000007LL;

template<class ll>
inline ll getCeilExponent(ll a)
{
	if (a <= 0)
	{
		throw invalid_argument("Vaule equal to or less than 0 should not be valid");
	}

	ll index = 0;
	--a;
	for (; a != 0; a >>= 1)
	{
		++index;
	}
	return index;
}

template<class ll>
inline ll getFloorExponent(ll a)
{
	if (a <= 0)
	{
		throw invalid_argument("Vaule equal to or less than 0 should not be valid");
	}

	ll index = 0;
	a >>= 1;
	for (; a != 0; a >>= 1)
	{
		++index;
	}
	return index;
}


map<ll, ull> getPrimeFactorization(ll num)
{
	map<ll, ull> primeToCount;

	int rootNum = sqrt(num);
	vector<bool> isCompositeNumber(rootNum + 1, false);
	for (int p = 2; num != 1 && p <= rootNum; ++p)
	{
		if (!isCompositeNumber[p])
		{
			for (int i = p * 2; i <= rootNum; i += p)
			{
				isCompositeNumber[i] = true;
			}

			int count = 0;
			while (num % p == 0)
			{
				++count;
				num /= p;
			}
			if (count >= 1)
			{
				primeToCount.insert(make_pair(p, count));
			}
		}
	}
	if (num != 1)
	{
		primeToCount.insert(make_pair(num, 1));
	}

	return primeToCount;
}

// Computational Complexity: O(log(max(a, b))) 
template<class ll>
inline ll getGCD(ll a, ll b)
{
	if (b == 0)
	{
		return a;
	}
	return getGCD(b, a % b);
}

// Computational Complexity: O(log(max(a, b))) 
template<class ll>
inline pair<ll, ll> getBezoutsIdentitySoulution(ll a, ll b)
{
	if (b == 0)
	{
		return { 1, 0 };
	}
	auto sol = getBezoutsIdentitySoulution(b, a % b);
	return { sol.second, sol.first - (a / b) * sol.second };
}

inline ll getModValue(ll a, ll m)
{
	return (a % m + m) % m;
}

// Computational Complexity: O(log(max(a, n))) 
inline ll getModInverse(ll a, ll m)
{
	auto sol = getBezoutsIdentitySoulution(a, m);
	if (a * sol.first + m * sol.second != 1)
	{
		return -1;
	}
	return getModValue(sol.first, m);
}

ll getModCombinatorialNumber(ll n, ll r, ll m)
{
	ll res = 1;
	for (; n > 0; --n, --r)
	{
		res = (((res * n) % m) * getModInverse(n, m)) % m;
	}
	return getModValue(res, m);
}

ll getModPower(ll base, unsigned t, ll m)
{
	if (t == 0)
	{
		return 1LL;
	}

	auto res = getModPower(base, t >> 1, m);
	res = getModValue(res * res, m);
	if (t % 2 == 1)
	{
		res = getModValue(res * base, m);
	}
	return res;
}

int main()
{
	cin.tie(0);
	ios::sync_with_stdio(false);

	ll n, k;
	cin >> n >> k;

	vector<pair<ll, ll>> sushis(n);
	ll t, d;
	for (ll i = 0; i < n; ++i)
	{
		cin >> t >> d;
		sushis[i] = {t, d};
	}

	sort(sushis.begin(), sushis.end(), [](pair<ll, ll> a, pair<ll, ll> b) { return a.second > b.second; });

	priority_queue<ll, vector<ll>, greater<ll>> qs;
	set<ll> in_sum;
	set<ll> type;
	vector<ll> max_ds(n + 1);
	ll sum_d = 0;
	for (ll i = 0; i < n; ++i)
	{
		t = sushis[i].first;
		d = sushis[i].second;

		max_ds[t] = max(max_ds[t], d);
		type.insert(t);
		if (i >= k)
		{
			continue;
		}
		sum_d += sushis[i].second;

		if (in_sum.count(t) == 1)
		{
			qs.push(d);
		}
		in_sum.insert(t);
	}
	sort(max_ds.begin(), max_ds.end(), greater<ll>());

	ll res = sum_d + in_sum.size() * in_sum.size();
	for (ll i = in_sum.size() + 1; i <= type.size() && !qs.empty(); ++i)
	{
		sum_d -= qs.top();
		qs.pop();
		sum_d += max_ds[i - 1];
		res = max(res, sum_d + i * i);
	}

	cout << res;

	return 0;
}




Measured memory_bytes: 9600