The construction of subway in Bertown is almost finished! The President of Berland will visit this city soon to look at the new subway himself.
There are n stations in the subway. It was built according to the Bertown Transport Law:
- For each station i there exists exactly one train that goes from this station. Its destination station is pi, possibly pi = i;
- For each station i there exists exactly one station j such that pj = i.
The President will consider the convenience of subway after visiting it. The convenience is the number of ordered pairs (x, y) such that person can start at station x and, after taking some subway trains (possibly zero), arrive at station y (1 ≤ x, y ≤ n).
The mayor of Bertown thinks that if the subway is not convenient enough, then the President might consider installing a new mayor (and, of course, the current mayor doesn't want it to happen). Before President visits the city mayor has enough time to rebuild some paths of subway, thus changing the values of pi for not more than two subway stations. Of course, breaking the Bertown Transport Law is really bad, so the subway must be built according to the Law even after changes.
The mayor wants to do these changes in such a way that the convenience of the subway is maximized. Help him to calculate the maximum possible convenience he can get!
The first line contains one integer number n (1 ≤ n ≤ 100000) — the number of stations.
The second line contains n integer numbers p1, p2, ..., pn (1 ≤ pi ≤ n) — the current structure of the subway. All these numbers are distinct.
Print one number — the maximum possible value of convenience.
3 2 1 3
9
5 1 5 4 3 2
17
In the first example the mayor can change p2 to 3 and p3 to 1, so there will be 9 pairs: (1, 1), (1, 2), (1, 3), (2, 1), (2, 2), (2, 3), (3, 1), (3, 2), (3, 3).
In the second example the mayor can change p2 to 4 and p3 to 5.
代码:
#include<iostream>
#include<vector>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn = 1e5 + 5;
ll vis[maxn];
vector<ll>num;
vector<ll>g[maxn];
int main()
{
ll n; cin >> n;
for (int i = 1; i <= n; i++)
{
ll v; cin >> v;
g[i].push_back(v);
}
for (int i = 1; i <= n; i++)
{
if (vis[i] == 0)
{
vis[i] = 1;
ll j = i, count = 1;
while ((g[j][0] != i))
vis[g[j][0]] = 1, j = g[j][0], count++;
num.push_back(count);
}
}
ll len = num.size(), sum = 0;
sort(num.begin(), num.end());
if (len == 1)
{
cout << num[0]*num[0] << endl;
}
else
{
for (int i = 0; i < len - 2; i++)
sum += num[i] * num[i];
sum += (num[len - 2] + num[len - 1])*(num[len - 2] + num[len - 1]);
cout << sum << endl;
}
return 0;
}