結果

問題 No.387 ハンコ
ユーザー yosupotyosupot
提出日時 2016-06-23 16:02:30
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 2,242 bytes
コンパイル時間 1,076 ms
コンパイル使用メモリ 82,124 KB
実行使用メモリ 9,216 KB
最終ジャッジ日時 2024-10-11 20:10:12
合計ジャッジ時間 39,589 ms
ジャッジサーバーID
(参考情報)
judge5 / judge4
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 WA -
testcase_02 WA -
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <iostream>
#include <vector>
#include <cstring>
#include <algorithm>
#include <queue>
#include <cmath>
#include <bitset>
#include <cassert>
#include <array>

using namespace std;
using ll = long long;
using ull = unsigned long long;
using uint = unsigned int;

const int MN = 201024;

using B = array<ull, MN/64>;
const int BS = MN/64;

void orassign(B &l, B r) {
    for (int i = 0; i < BS; i++) {
        l[i] |= r[i];
    }
}
void xorassign(B &l, B r) {
    for (int i = 0; i < BS; i++) {
        l[i] ^= r[i];
    }
}

B shiftL(const B &x, int k) { // <<
    B ans; ans.fill(0);
    int shift = k%64;
    if (shift) {
        ans[k/64] |= x[0]<<shift;
        for (int i = 1; i < BS-k/64; i++) {
            ans[i+k/64] |= x[i-1]>>(64-shift);
            ans[i+k/64] |= x[i]<<shift;
        }
    } else {
        for (int i = 0; i < BS-k/64; i++) {
            ans[i+k/64] |= x[i];
        }
    }
    return ans;
}
int main() {
    constexpr int MD = 100128;
    vector<int> g[MD];
    int n;
    scanf("%d", &n);
    for (int i = 0; i < n; i++) {
        int a;
        scanf("%d", &a);
        g[a].push_back(i);
    }
    B mp; mp.fill(0);
    for (int i = 0; i < n; i++) {
        int x;
        scanf("%d", &x);
        if (x) {
            mp[i/64] |= 1ULL<<(i%64);
        }
    }
    int ans[MN] = {};
//    B ans; ans.fill(0);
    constexpr int U = 4;
    constexpr int BU = 15;
    constexpr ull bmap = 0x0101010101010101;
    for (int fe = 0; fe < MD/BU; fe++) {
        B res[BU];
        for (int i = 0; i < BU; i++) {
            res[i].fill(0);
            for (int d: g[fe*BU+i]) {
                orassign(res[i], shiftL(mp, d));
            }
        }
        for (int i = 0; i < BS; i++) {
            for (int b = 0; b < U; b++) {
                ull sm = 0;
                for (int j = 0; j < BU; j++) {
                    sm += ((res[j][i]>>b) & bmap);
                }
                for (int j = 0; j < 64/U; j++) {
                    ans[i*64+j*U+b] += (sm>>(j*U))&((1<<U)-1);
                }
            }
        }
    }
    for (int i = 0; i < 2*n-1; i++) {
        if (ans[i] % 2) {
            printf("ODD\n");
        } else {
            printf("EVEN\n");
        }
    }
    return 0;
}
0