#include <cstdio>
#include <cassert>
#include <algorithm>
const int MAXN = 3000;
int a[MAXN+5];
int sorted[MAXN+5];
int cel[MAXN+5];
int mapa[MAXN+5];
bool wolny[MAXN+5];
int zamiany[MAXN+5][2];
int ppol;
int cykl[MAXN+5];
int ncykl;
#define NDEBUG
int main() {
int n;
scanf("%d", &n);
for (int i = 0; i < n; ++i) {
scanf("%d", a+i);
sorted[i] = a[i];
assert(a[i] <= 3000 && a[i] >= 1);
}
std::sort(sorted, sorted+n);
for (int i = 0; i < n; ++i) {
mapa[sorted[i]] = i;
}
for (int i = 0; i < n; ++i) {
cel[i] = mapa[a[i]];
}
for (int i = 0; i < n; ++i) {
wolny[i] = true;
}
int kolejny, dlmax = 0;
for (int i = 0; i < n; ++i) {
ncykl = 0;
if (wolny[i]) {
kolejny = i;
while(wolny[kolejny]){
cykl[ncykl++] = kolejny;
wolny[kolejny] = false;
kolejny = cel[kolejny];
}
if (ncykl > 1) {
int pary = ncykl / 2;
for (int i = 0; i < pary; ++i) {
zamiany[ppol][0] = cykl[i];
zamiany[ppol][1] = cykl[ncykl-i-1];
++ppol;
}
}
dlmax = std::max(ncykl, dlmax);
}
}
printf("%d\n", dlmax > 2 ? 2 : (dlmax == 2 ? 1 : 0));
int cnt = 0;
while(true) {
if (ppol == 0) {
break;
} else {
++cnt;
printf("%d\n", 2*ppol);
for (int i = 0; i < ppol; ++i) {
printf("%d ", zamiany[i][0]+1);
}
for (int i = ppol-1; i >= 0; --i) {
printf("%d ", zamiany[i][1]+1);
}
putchar('\n');
for (int i = 0; i < ppol; ++i) {
std::swap(a[zamiany[i][0]], a[zamiany[i][1]]);
std::swap(cel[zamiany[i][0]], cel[zamiany[i][1]]);
}
}
for (int i = 0; i < n; ++i) {
wolny[i] = true;
}
ppol = 0;
for (int i = 0; i < n; ++i) {
if (cel[i] == i)
continue;
if (wolny[i] && wolny[cel[i]]) {
zamiany[ppol][0] = i;
zamiany[ppol][1] = cel[i];
++ppol;
wolny[i] = wolny[cel[i]] = false;
}
}
}
if (dlmax == 2) {
assert (cnt == 1);
} else if (dlmax == 1) {
assert (cnt == 0);
} else {
assert (cnt == 2);
}
}
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 | #include <cstdio> #include <cassert> #include <algorithm> const int MAXN = 3000; int a[MAXN+5]; int sorted[MAXN+5]; int cel[MAXN+5]; int mapa[MAXN+5]; bool wolny[MAXN+5]; int zamiany[MAXN+5][2]; int ppol; int cykl[MAXN+5]; int ncykl; #define NDEBUG int main() { int n; scanf("%d", &n); for (int i = 0; i < n; ++i) { scanf("%d", a+i); sorted[i] = a[i]; assert(a[i] <= 3000 && a[i] >= 1); } std::sort(sorted, sorted+n); for (int i = 0; i < n; ++i) { mapa[sorted[i]] = i; } for (int i = 0; i < n; ++i) { cel[i] = mapa[a[i]]; } for (int i = 0; i < n; ++i) { wolny[i] = true; } int kolejny, dlmax = 0; for (int i = 0; i < n; ++i) { ncykl = 0; if (wolny[i]) { kolejny = i; while(wolny[kolejny]){ cykl[ncykl++] = kolejny; wolny[kolejny] = false; kolejny = cel[kolejny]; } if (ncykl > 1) { int pary = ncykl / 2; for (int i = 0; i < pary; ++i) { zamiany[ppol][0] = cykl[i]; zamiany[ppol][1] = cykl[ncykl-i-1]; ++ppol; } } dlmax = std::max(ncykl, dlmax); } } printf("%d\n", dlmax > 2 ? 2 : (dlmax == 2 ? 1 : 0)); int cnt = 0; while(true) { if (ppol == 0) { break; } else { ++cnt; printf("%d\n", 2*ppol); for (int i = 0; i < ppol; ++i) { printf("%d ", zamiany[i][0]+1); } for (int i = ppol-1; i >= 0; --i) { printf("%d ", zamiany[i][1]+1); } putchar('\n'); for (int i = 0; i < ppol; ++i) { std::swap(a[zamiany[i][0]], a[zamiany[i][1]]); std::swap(cel[zamiany[i][0]], cel[zamiany[i][1]]); } } for (int i = 0; i < n; ++i) { wolny[i] = true; } ppol = 0; for (int i = 0; i < n; ++i) { if (cel[i] == i) continue; if (wolny[i] && wolny[cel[i]]) { zamiany[ppol][0] = i; zamiany[ppol][1] = cel[i]; ++ppol; wolny[i] = wolny[cel[i]] = false; } } } if (dlmax == 2) { assert (cnt == 1); } else if (dlmax == 1) { assert (cnt == 0); } else { assert (cnt == 2); } } |
English