QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#132397 | #6329. Colorful Graph | UrgantTeam# | TL | 451ms | 16808kb | C++23 | 3.2kb | 2023-07-29 19:59:02 | 2023-07-29 19:59:05 |
Judging History
answer
#pragma GCC optimize("O3,unroll-loops")
#pragma GCC target("avx2,bmi,bmi2,lzcnt,popcnt")
#include<bits/stdc++.h>
using namespace std;
mt19937 rnd(time(0));
int const maxn = 7005;
vector < int > g[maxn], G[maxn], ts, go[maxn], res[2 * maxn];
int used[2 * maxn], cmp, numb[maxn], bad;
vector < int > in[maxn];
bitset < maxn > Gz[maxn];
int one[maxn], two[maxn], out[maxn];
bitset < maxn > free_righ, USED;
void answer(int v, int mx, int cmp) {
used[v] = bad;
if (v <= cmp) {
for (auto key : in[v]) out[key] = mx;
}
for (auto key : res[v]) {
if (used[key] != bad) answer(key, mx, cmp);
}
}
void dfs1(int v) {
used[v] = 2, numb[v] = cmp, in[cmp].push_back(v);
for (auto u : g[v]) {
if (used[u] == 1) {
dfs1(u);
} else if (numb[u] != cmp) {
go[cmp].push_back(numb[u]);
}
}
}
void dfs(int v) {
used[v] = 1;
for (auto u : G[v]) if (!used[u]) dfs(u);
ts.push_back(v);
}
void dfs2(int v, int start) {
used[v] = bad;
if (v != start) Gz[start][v] = 1;
for (auto u : go[v]) {
if (used[u] != bad) dfs2(u, start);
}
}
bool calc(int v) {
used[v] = bad;
int pos = (Gz[v]&free_righ)._Find_first();
if (pos <= cmp) {
one[v] = pos, two[pos] = v, free_righ[pos] = 0;
return true;
}
for (int u = Gz[v]._Find_first(); u <= cmp; u = Gz[v]._Find_next(u)) {
if (used[two[u]] != bad && calc(two[u])) {
one[v] = u, two[u] = v;
return true;
}
}
return false;
}
main() {
#ifdef HOME
freopen("input.txt", "r", stdin);
#endif // HOME
ios_base::sync_with_stdio(0);
cin.tie(0);
cout.tie(0);
int n, m, u, v;
cin >> n >> m;
for (int i = 1; i <= m; i++) {
cin >> u >> v;
g[u].push_back(v);
G[v].push_back(u);
}
for (int i = 1; i <= n; i++) if (!used[i]) dfs(i);
reverse(ts.begin(), ts.end());
for (auto key : ts) {
if (used[key] != 2) {
cmp++;
dfs1(key);
}
}
bad = 2;
for (int i = 1; i <= cmp; i++) {
bad++;
dfs2(i, i);
}
int find = 1;
for (int i = 1; i <= cmp; i++) free_righ[i] = 1;
vector < int > Q;
for (int i = 1; i <= cmp; i++) {
Q.push_back(i);
}
shuffle(Q.begin(), Q.end(), rnd);
for (auto i : Q) {
for (auto j : Q) {
if (Gz[i][j] && !one[i] && !two[j]) {
one[i] = j, two[j] = i, free_righ[j] = 0;
}
}
}
while (find) {
find = 0;
bad++;
for (auto i : Q) {
if (used[i] != bad && !one[i] && calc(i)) {
find = 1;
}
}
}
int mx = 0;
bad++;
for (int i = 1; i <= cmp; i++) {
if (one[i]) res[i].push_back(one[i] + cmp), res[one[i] + cmp].push_back(i);
res[i].push_back(i + cmp), res[i + cmp].push_back(i);
}
for (int i = 1; i <= cmp; i++) {
if (used[i] != bad) {
mx++;
answer(i, mx, cmp);
}
}
for (int i = 1; i <= n; i++) cout << out[i] << " ";
cout << '\n';
return 0;
}
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 2ms
memory: 6600kb
input:
5 5 1 4 2 3 1 3 2 5 5 1
output:
1 1 2 1 1
result:
ok AC
Test #2:
score: 0
Accepted
time: 0ms
memory: 6084kb
input:
5 7 1 2 2 1 4 3 5 1 5 4 4 1 4 5
output:
2 2 1 1 1
result:
ok AC
Test #3:
score: 0
Accepted
time: 2ms
memory: 4576kb
input:
8 6 6 1 3 4 3 6 2 3 4 1 6 4
output:
4 4 4 4 3 4 2 1
result:
ok AC
Test #4:
score: 0
Accepted
time: 451ms
memory: 16808kb
input:
7000 6999 4365 4296 2980 3141 6820 4995 4781 24 2416 5844 2940 2675 3293 2163 3853 5356 262 6706 1985 1497 5241 3803 353 1624 5838 4708 5452 3019 2029 6161 3849 4219 1095 1453 4268 4567 1184 1857 2911 3977 1662 2751 6353 6496 2002 6628 1407 4623 425 1331 4445 4277 1259 3165 4994 1044 2756 5788 5496 ...
output:
752 1629 733 564 1682 1208 1278 1052 231 1750 1369 1716 803 1708 1434 607 1749 203 513 1748 626 1747 427 1282 466 99 1112 110 1174 1257 1746 953 1673 140 907 766 1745 942 157 397 1744 1305 1438 1743 569 1742 464 165 338 1571 1064 919 1741 564 1742 1740 173 768 1739 613 1049 269 141 1092 1738 50 474 ...
result:
ok AC
Test #5:
score: 0
Accepted
time: 332ms
memory: 15548kb
input:
7000 6999 4832 1603 5984 6985 5355 3687 6007 2170 5984 3486 3267 2189 538 2123 4343 4553 5855 6168 5984 257 4239 2304 5984 2063 3298 1869 5984 6353 5984 2018 5984 5387 5984 3382 3164 3978 2690 2816 4810 2638 5984 3773 5984 1634 5984 2786 5984 3671 5984 5140 2943 5721 5984 414 1105 4060 3093 796 5984...
output:
1719 52 2333 2332 2331 1602 2330 2329 1441 77 1715 1991 1526 2328 1514 1868 1092 2100 804 2327 1242 2326 2325 2324 1831 2221 2323 2322 2321 2320 2319 1122 340 1857 1743 517 2318 524 2317 2316 1659 2315 2314 1972 2313 1618 2312 25 930 2311 2310 2309 2308 2307 2306 217 1704 1925 2237 284 1029 2299 118...
result:
ok AC
Test #6:
score: 0
Accepted
time: 168ms
memory: 11780kb
input:
7000 6999 1649 5337 1701 3344 4394 2172 3330 39 5932 1141 5381 5340 5453 3300 125 2172 6810 5263 804 2172 6635 2172 676 4740 3015 1183 1710 5769 611 5915 3419 1581 2094 2172 4508 2172 6604 2433 6113 1466 1604 696 1518 1123 1287 2940 4825 2172 5130 4524 2693 2172 106 2172 5157 2172 3693 2172 5198 217...
output:
778 380 2333 2332 508 217 2331 2330 2329 1623 893 465 1385 1432 16 2328 2157 2327 2326 2325 1926 513 1585 2324 1697 2323 1068 2322 156 1155 2321 276 2320 1281 2319 863 2094 2318 2317 2257 1386 2316 2315 2040 2304 1827 1547 2314 331 573 2313 1031 2312 2311 674 2310 338 2025 2309 561 1216 2133 2016 23...
result:
ok AC
Test #7:
score: 0
Accepted
time: 193ms
memory: 11952kb
input:
7000 6999 2896 6321 881 2623 5058 2623 4833 2623 4669 2623 4781 5007 1447 2623 4781 4768 4781 3834 2758 4792 797 5055 3784 2623 4781 5510 6606 3040 597 3459 4136 2037 1291 3989 4781 837 4781 4379 5637 2053 1642 2665 4781 4664 4781 952 4924 2511 4781 4201 4781 2352 4781 5362 3901 197 137 2623 2706 19...
output:
1454 1413 1079 1516 879 660 1353 314 1750 1424 1120 373 920 1456 479 1749 1748 682 224 623 1747 967 1746 1745 1744 1743 1742 1741 727 1438 1695 365 1542 1264 1318 1036 1740 1739 857 1738 167 898 1737 750 263 869 1726 119 1166 1483 526 793 620 1449 1279 956 442 354 1736 167 1139 1701 1347 1735 1230 1...
result:
ok AC
Test #8:
score: 0
Accepted
time: 65ms
memory: 10232kb
input:
6999 6998 1269 3969 1269 2429 1269 2609 1269 2515 1269 6166 1269 6614 3108 1269 2105 1269 4670 1269 578 1269 4661 1269 1421 1269 2576 1269 6152 1269 1269 6636 3011 1269 305 1269 5189 1269 1683 1269 6861 1269 1269 5798 1499 1269 282 1269 914 1269 80 1269 677 1269 701 1269 1269 359 6521 1269 1269 1754...
output:
565 1414 1891 3273 3499 1756 602 3498 3497 3496 3495 3494 2721 3493 2590 2742 1902 3492 749 3206 3491 124 144 3490 32 431 1112 1325 2123 3489 814 3488 1888 2343 3487 3229 3486 3485 2743 3484 3483 2360 1697 3482 3481 1584 1004 2010 3496 3015 3480 3479 3478 3072 3477 3476 476 863 1232 1796 493 3475 34...
result:
ok AC
Test #9:
score: 0
Accepted
time: 27ms
memory: 6812kb
input:
7000 0
output:
7000 6999 6998 6997 6996 6995 6994 6993 6992 6991 6990 6989 6988 6987 6986 6985 6984 6983 6982 6981 6980 6979 6978 6977 6976 6975 6974 6973 6972 6971 6970 6969 6968 6967 6966 6965 6964 6963 6962 6961 6960 6959 6958 6957 6956 6955 6954 6953 6952 6951 6950 6949 6948 6947 6946 6945 6944 6943 6942 6941 ...
result:
ok AC
Test #10:
score: -100
Time Limit Exceeded
input:
7000 6999 3138 1903 3285 5919 6182 1430 1164 961 1577 6445 1390 3384 935 5723 6614 6387 4799 2877 3915 5128 5366 5455 2287 3941 2053 2326 4022 6993 488 2922 4327 4701 4674 3221 1666 4773 4356 3232 3888 937 4318 6942 577 1299 4491 1938 5154 1254 790 5532 4286 5478 2918 6725 2853 304 2554 5207 5140 77...