QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#132397#6329. Colorful GraphUrgantTeam#TL 451ms16808kbC++233.2kb2023-07-29 19:59:022023-07-29 19:59:05

Judging History

你现在查看的是最新测评结果

  • [2023-08-10 23:21:45]
  • System Update: QOJ starts to keep a history of the judgings of all the submissions.
  • [2023-07-29 19:59:05]
  • 评测
  • 测评结果:TL
  • 用时:451ms
  • 内存:16808kb
  • [2023-07-29 19:59:02]
  • 提交

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...

output:


result: