Submission #572394

# Submission time Handle Problem Language Result Execution time Memory
572394 2022-06-04T10:33:10 Z piOOE Cubeword (CEOI19_cubeword) C++17
84 / 100
218 ms 13368 KB
#include <bits/stdc++.h>

using namespace std;

#define sz(x) ((int)size(x))
#define all(x) begin(x), end(x)
#define trace(x) cout << #x << ": " << (x) << endl;

namespace Ment {
    template<typename T>
    T inverse(T a, T m) {
        T u = 0, v = 1;
        while (a != 0) {
            T t = m / a;
            m -= t * a;
            swap(a, m);
            u -= t * v;
            swap(u, v);
        }
        assert(m == 1);
        return u;
    }

    template<typename T>
    class Modular {
    public:
        using Type = typename decay<decltype(T::value)>::type;

        constexpr Modular() : value() {}

        template<typename U>
        Modular(const U &x) {
            value = normalize(x);
        }

        template<typename U>
        static Type normalize(const U &x) {
            Type v;
            if (-mod() <= x && x < mod()) v = static_cast<Type>(x);
            else v = static_cast<Type>(x % mod());
            if (v < 0) v += mod();
            return v;
        }

        const Type &operator()() const { return value; }

        template<typename U>
        explicit operator U() const { return static_cast<U>(value); }

        constexpr static Type mod() { return T::value; }

        Modular &operator+=(const Modular &other) {
            value += other.value - mod();
            value += (value >> 31) & mod();
            return *this;
        }

        Modular &operator-=(const Modular &other) {
            value -= other.value;
            value += (value >> 31) & mod();
            return *this;
        }

        template<typename U>
        Modular &operator+=(const U &other) { return *this += Modular(other); }

        template<typename U>
        Modular &operator-=(const U &other) { return *this -= Modular(other); }

        Modular &operator++() { return *this += 1; }

        Modular &operator--() { return *this -= 1; }

        Modular operator++(int) {
            Modular result(*this);
            *this += 1;
            return result;
        }

        Modular operator--(int) {
            Modular result(*this);
            *this -= 1;
            return result;
        }

        Modular operator-() const { return Modular(-value); }

        template<typename U = T>
        typename enable_if<is_same<typename Modular<U>::Type, int>::value, Modular>::type &
        operator*=(const Modular &rhs) {
#ifdef _WIN32
            uint64_t x = static_cast<int64_t>(value) * static_cast<int64_t>(rhs.value);
            uint32_t xh = static_cast<uint32_t>(x >> 32), xl = static_cast<uint32_t>(x), d, m;
            asm(
                    "divl %4; \n\t"
                    : "=a" (d), "=d" (m)
                    : "d" (xh), "a" (xl), "r" (mod())
                    );
            value = m;
#else
            value = normalize(static_cast<int64_t>(value) * static_cast<int64_t>(rhs.value));
#endif
            return *this;
        }

        template<typename U = T>
        typename enable_if<is_same<typename Modular<U>::Type, int64_t>::value, Modular>::type &
        operator*=(const Modular &rhs) {
            int64_t q = static_cast<int64_t>(static_cast<long double>(value) * rhs.value / mod());
            value = normalize(value * rhs.value - q * mod());
            return *this;
        }

        template<typename U = T>
        typename enable_if<!is_integral<typename Modular<U>::Type>::value, Modular>::type &
        operator*=(const Modular &rhs) {
            value = normalize(value * rhs.value);
            return *this;
        }

        Modular &operator/=(const Modular &other) { return *this *= Modular(inverse(other.value, mod())); }

        template<typename U>
        friend bool operator==(const Modular<U> &lhs, const Modular<U> &rhs);

        template<typename U>
        friend bool operator<(const Modular<U> &lhs, const Modular<U> &rhs);

        template<typename U>
        friend std::istream &operator>>(std::istream &stream, Modular<U> &number);

    private:
        Type value;
    };

    template<typename T>
    bool operator==(const Modular<T> &lhs, const Modular<T> &rhs) { return lhs.value == rhs.value; }

    template<typename T, typename U>
    bool operator==(const Modular<T> &lhs, U rhs) { return lhs == Modular<T>(rhs); }

    template<typename T, typename U>
    bool operator==(U lhs, const Modular<T> &rhs) { return Modular<T>(lhs) == rhs; }

    template<typename T>
    bool operator!=(const Modular<T> &lhs, const Modular<T> &rhs) { return !(lhs == rhs); }

    template<typename T, typename U>
    bool operator!=(const Modular<T> &lhs, U rhs) { return !(lhs == rhs); }

    template<typename T, typename U>
    bool operator!=(U lhs, const Modular<T> &rhs) { return !(lhs == rhs); }

    template<typename T>
    bool operator<(const Modular<T> &lhs, const Modular<T> &rhs) { return lhs.value < rhs.value; }

    template<typename T>
    Modular<T> operator+(const Modular<T> &lhs, const Modular<T> &rhs) { return Modular<T>(lhs) += rhs; }

    template<typename T, typename U>
    Modular<T> operator+(const Modular<T> &lhs, U rhs) { return Modular<T>(lhs) += rhs; }

    template<typename T, typename U>
    Modular<T> operator+(U lhs, const Modular<T> &rhs) { return Modular<T>(lhs) += rhs; }

    template<typename T>
    Modular<T> operator-(const Modular<T> &lhs, const Modular<T> &rhs) { return Modular<T>(lhs) -= rhs; }

    template<typename T, typename U>
    Modular<T> operator-(const Modular<T> &lhs, U rhs) { return Modular<T>(lhs) -= rhs; }

    template<typename T, typename U>
    Modular<T> operator-(U lhs, const Modular<T> &rhs) { return Modular<T>(lhs) -= rhs; }

    template<typename T>
    Modular<T> operator*(const Modular<T> &lhs, const Modular<T> &rhs) { return Modular<T>(lhs) *= rhs; }

    template<typename T, typename U>
    Modular<T> operator*(const Modular<T> &lhs, U rhs) { return Modular<T>(lhs) *= rhs; }

    template<typename T, typename U>
    Modular<T> operator*(U lhs, const Modular<T> &rhs) { return Modular<T>(lhs) *= rhs; }

    template<typename T>
    Modular<T> operator/(const Modular<T> &lhs, const Modular<T> &rhs) { return Modular<T>(lhs) /= rhs; }

    template<typename T, typename U>
    Modular<T> operator/(const Modular<T> &lhs, U rhs) { return Modular<T>(lhs) /= rhs; }

    template<typename T, typename U>
    Modular<T> operator/(U lhs, const Modular<T> &rhs) { return Modular<T>(lhs) /= rhs; }

    template<typename T, typename U>
    Modular<T> power(const Modular<T> &a, const U &b) {
        assert(b >= 0);
        Modular<T> x = a, res = 1;
        U p = b;
        while (p > 0) {
            if (p & 1) res *= x;
            x *= x;
            p >>= 1;
        }
        return res;
    }

    template<typename T>
    string to_string(const Modular<T> &number) {
        return to_string(number());
    }

    template<typename T>
    std::ostream &operator<<(std::ostream &stream, const Modular<T> &number) {
        return stream << number();
    }

    template<typename T>
    std::istream &operator>>(std::istream &stream, Modular<T> &number) {
        typename common_type<typename Modular<T>::Type, int64_t>::type x;
        stream >> x;
        number.value = Modular<T>::normalize(x);
        return stream;
    }

    constexpr int md = 998244353;
    using Mint = Modular<std::integral_constant<decay<decltype(md)>::type, md>>;
}
using Ment::Mint;


typedef long long ll;

mt19937 rnd(chrono::steady_clock::now().time_since_epoch().count());

int rand(int l, int r) { return (int) ((ll) rnd() % (r - l + 1)) + l; }

const int N = 8, A = 32, L = 11;
const ll infL = 3e18;

int g[N][3] = {{1, 2, 4},
               {0, 3, 5},
               {0, 3, 6},
               {1, 2, 7},
               {0, 5, 6},
               {1, 4, 7},
               {2, 4, 7},
               {3, 5, 6}};

int n;
Mint cnt[A][A][L], dp2[A][A][A], dp7[A][A][A], dp1[A][A][A], dp4[A][A][A];
Mint ans = 0;
int len = 0, used[N];

int get(char c) {
    if ('a' <= c && c <= 'z') {
        return c - 'a';
    }
    if ('A' <= c && c <= 'Z') {
        return (c - 'A') + 16;
    }
    return (c - '0') + 52;
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cin >> n;
    vector<string> stt;
    for (int i = 0; i < n; ++i) {
        string s;
        cin >> s;
        string ss = s;
        reverse(all(ss));
        stt.push_back(ss);
        stt.push_back(s);
    }
    sort(all(stt));
    stt.resize(unique(all(stt)) - begin(stt));
    for (string s: stt) {
        ++cnt[get(s[0])][get(s.back())][sz(s)];
    }
    for (len = 1; len < L; ++len) {
        for (int x0 = 0; x0 < A; ++x0) {
            for (int x3 = 0; x3 < A; ++x3) {
                for (int x6 = 0; x6 < A; ++x6) {
                    Mint sum = 0;
                    for (int x = 0; x < A; ++x) {
                        sum += cnt[x][x0][len] * cnt[x][x3][len] * cnt[x][x6][len];
                    }
                    dp2[x0][x3][x6] = dp7[x0][x3][x6] = dp1[x0][x3][x6] = dp4[x0][x3][x6] = sum;
                }
            }
        }
        for (int x0 = 0; x0 < A; ++x0) {
            for (int x3 = 0; x3 < A; ++x3) {
                for (int x5 = 0; x5 < A; ++x5) {
                    for (int x6 = 0; x6 < A; ++x6) {
                        ans = (ans + dp4[x0][x5][x6] * dp7[x3][x5][x6] * dp1[x0][x3][x5] *
                                     dp2[x0][x3][x6]);
                    }
                }
            }
        }
    }
    cout << ans;
    return 0;
}
# Verdict Execution time Memory Grader output
1 Correct 163 ms 8648 KB Output is correct
2 Correct 151 ms 8580 KB Output is correct
3 Correct 149 ms 8648 KB Output is correct
4 Correct 152 ms 8648 KB Output is correct
5 Correct 173 ms 8648 KB Output is correct
6 Correct 151 ms 8676 KB Output is correct
7 Correct 145 ms 8640 KB Output is correct
8 Correct 156 ms 8648 KB Output is correct
9 Correct 162 ms 8648 KB Output is correct
10 Correct 147 ms 8648 KB Output is correct
11 Correct 150 ms 8600 KB Output is correct
12 Correct 171 ms 8648 KB Output is correct
13 Correct 157 ms 8612 KB Output is correct
14 Correct 147 ms 8576 KB Output is correct
15 Correct 194 ms 8648 KB Output is correct
16 Correct 149 ms 8692 KB Output is correct
17 Correct 152 ms 8648 KB Output is correct
18 Correct 149 ms 8648 KB Output is correct
19 Correct 161 ms 8612 KB Output is correct
20 Correct 158 ms 8608 KB Output is correct
21 Correct 150 ms 8648 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 163 ms 8648 KB Output is correct
2 Correct 151 ms 8580 KB Output is correct
3 Correct 149 ms 8648 KB Output is correct
4 Correct 152 ms 8648 KB Output is correct
5 Correct 173 ms 8648 KB Output is correct
6 Correct 151 ms 8676 KB Output is correct
7 Correct 145 ms 8640 KB Output is correct
8 Correct 156 ms 8648 KB Output is correct
9 Correct 162 ms 8648 KB Output is correct
10 Correct 147 ms 8648 KB Output is correct
11 Correct 150 ms 8600 KB Output is correct
12 Correct 171 ms 8648 KB Output is correct
13 Correct 157 ms 8612 KB Output is correct
14 Correct 147 ms 8576 KB Output is correct
15 Correct 194 ms 8648 KB Output is correct
16 Correct 149 ms 8692 KB Output is correct
17 Correct 152 ms 8648 KB Output is correct
18 Correct 149 ms 8648 KB Output is correct
19 Correct 161 ms 8612 KB Output is correct
20 Correct 158 ms 8608 KB Output is correct
21 Correct 150 ms 8648 KB Output is correct
22 Correct 148 ms 8648 KB Output is correct
23 Correct 158 ms 8592 KB Output is correct
24 Correct 147 ms 8604 KB Output is correct
25 Correct 152 ms 8648 KB Output is correct
26 Correct 145 ms 8648 KB Output is correct
27 Correct 165 ms 8636 KB Output is correct
28 Correct 152 ms 8656 KB Output is correct
29 Correct 153 ms 8648 KB Output is correct
30 Correct 160 ms 8648 KB Output is correct
31 Correct 161 ms 8648 KB Output is correct
32 Correct 146 ms 8648 KB Output is correct
33 Correct 147 ms 8648 KB Output is correct
34 Correct 151 ms 8648 KB Output is correct
35 Correct 155 ms 8696 KB Output is correct
36 Correct 147 ms 8648 KB Output is correct
37 Correct 151 ms 8680 KB Output is correct
38 Correct 157 ms 8580 KB Output is correct
39 Correct 145 ms 8648 KB Output is correct
40 Correct 155 ms 8648 KB Output is correct
41 Correct 158 ms 8628 KB Output is correct
42 Correct 155 ms 8648 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 163 ms 8648 KB Output is correct
2 Correct 151 ms 8580 KB Output is correct
3 Correct 149 ms 8648 KB Output is correct
4 Correct 152 ms 8648 KB Output is correct
5 Correct 173 ms 8648 KB Output is correct
6 Correct 151 ms 8676 KB Output is correct
7 Correct 145 ms 8640 KB Output is correct
8 Correct 156 ms 8648 KB Output is correct
9 Correct 162 ms 8648 KB Output is correct
10 Correct 147 ms 8648 KB Output is correct
11 Correct 150 ms 8600 KB Output is correct
12 Correct 171 ms 8648 KB Output is correct
13 Correct 157 ms 8612 KB Output is correct
14 Correct 147 ms 8576 KB Output is correct
15 Correct 194 ms 8648 KB Output is correct
16 Correct 149 ms 8692 KB Output is correct
17 Correct 152 ms 8648 KB Output is correct
18 Correct 149 ms 8648 KB Output is correct
19 Correct 161 ms 8612 KB Output is correct
20 Correct 158 ms 8608 KB Output is correct
21 Correct 150 ms 8648 KB Output is correct
22 Correct 148 ms 8648 KB Output is correct
23 Correct 158 ms 8592 KB Output is correct
24 Correct 147 ms 8604 KB Output is correct
25 Correct 152 ms 8648 KB Output is correct
26 Correct 145 ms 8648 KB Output is correct
27 Correct 165 ms 8636 KB Output is correct
28 Correct 152 ms 8656 KB Output is correct
29 Correct 153 ms 8648 KB Output is correct
30 Correct 160 ms 8648 KB Output is correct
31 Correct 161 ms 8648 KB Output is correct
32 Correct 146 ms 8648 KB Output is correct
33 Correct 147 ms 8648 KB Output is correct
34 Correct 151 ms 8648 KB Output is correct
35 Correct 155 ms 8696 KB Output is correct
36 Correct 147 ms 8648 KB Output is correct
37 Correct 151 ms 8680 KB Output is correct
38 Correct 157 ms 8580 KB Output is correct
39 Correct 145 ms 8648 KB Output is correct
40 Correct 155 ms 8648 KB Output is correct
41 Correct 158 ms 8628 KB Output is correct
42 Correct 155 ms 8648 KB Output is correct
43 Correct 187 ms 8648 KB Output is correct
44 Correct 188 ms 8648 KB Output is correct
45 Correct 197 ms 8648 KB Output is correct
46 Correct 183 ms 8648 KB Output is correct
47 Correct 185 ms 8648 KB Output is correct
48 Correct 190 ms 8608 KB Output is correct
49 Correct 196 ms 8612 KB Output is correct
50 Correct 182 ms 8592 KB Output is correct
51 Correct 185 ms 8648 KB Output is correct
52 Correct 216 ms 8680 KB Output is correct
53 Correct 185 ms 8648 KB Output is correct
54 Correct 193 ms 8648 KB Output is correct
55 Correct 218 ms 8648 KB Output is correct
56 Correct 181 ms 8648 KB Output is correct
57 Correct 187 ms 8648 KB Output is correct
58 Correct 194 ms 8620 KB Output is correct
59 Correct 188 ms 8688 KB Output is correct
60 Correct 187 ms 8648 KB Output is correct
61 Correct 203 ms 8648 KB Output is correct
62 Correct 211 ms 8648 KB Output is correct
63 Correct 185 ms 8608 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 163 ms 8648 KB Output is correct
2 Correct 151 ms 8580 KB Output is correct
3 Correct 149 ms 8648 KB Output is correct
4 Correct 152 ms 8648 KB Output is correct
5 Correct 173 ms 8648 KB Output is correct
6 Correct 151 ms 8676 KB Output is correct
7 Correct 145 ms 8640 KB Output is correct
8 Correct 156 ms 8648 KB Output is correct
9 Correct 162 ms 8648 KB Output is correct
10 Correct 147 ms 8648 KB Output is correct
11 Correct 150 ms 8600 KB Output is correct
12 Correct 171 ms 8648 KB Output is correct
13 Correct 157 ms 8612 KB Output is correct
14 Correct 147 ms 8576 KB Output is correct
15 Correct 194 ms 8648 KB Output is correct
16 Correct 149 ms 8692 KB Output is correct
17 Correct 152 ms 8648 KB Output is correct
18 Correct 149 ms 8648 KB Output is correct
19 Correct 161 ms 8612 KB Output is correct
20 Correct 158 ms 8608 KB Output is correct
21 Correct 150 ms 8648 KB Output is correct
22 Correct 148 ms 8648 KB Output is correct
23 Correct 158 ms 8592 KB Output is correct
24 Correct 147 ms 8604 KB Output is correct
25 Correct 152 ms 8648 KB Output is correct
26 Correct 145 ms 8648 KB Output is correct
27 Correct 165 ms 8636 KB Output is correct
28 Correct 152 ms 8656 KB Output is correct
29 Correct 153 ms 8648 KB Output is correct
30 Correct 160 ms 8648 KB Output is correct
31 Correct 161 ms 8648 KB Output is correct
32 Correct 146 ms 8648 KB Output is correct
33 Correct 147 ms 8648 KB Output is correct
34 Correct 151 ms 8648 KB Output is correct
35 Correct 155 ms 8696 KB Output is correct
36 Correct 147 ms 8648 KB Output is correct
37 Correct 151 ms 8680 KB Output is correct
38 Correct 157 ms 8580 KB Output is correct
39 Correct 145 ms 8648 KB Output is correct
40 Correct 155 ms 8648 KB Output is correct
41 Correct 158 ms 8628 KB Output is correct
42 Correct 155 ms 8648 KB Output is correct
43 Correct 187 ms 8648 KB Output is correct
44 Correct 188 ms 8648 KB Output is correct
45 Correct 197 ms 8648 KB Output is correct
46 Correct 183 ms 8648 KB Output is correct
47 Correct 185 ms 8648 KB Output is correct
48 Correct 190 ms 8608 KB Output is correct
49 Correct 196 ms 8612 KB Output is correct
50 Correct 182 ms 8592 KB Output is correct
51 Correct 185 ms 8648 KB Output is correct
52 Correct 216 ms 8680 KB Output is correct
53 Correct 185 ms 8648 KB Output is correct
54 Correct 193 ms 8648 KB Output is correct
55 Correct 218 ms 8648 KB Output is correct
56 Correct 181 ms 8648 KB Output is correct
57 Correct 187 ms 8648 KB Output is correct
58 Correct 194 ms 8620 KB Output is correct
59 Correct 188 ms 8688 KB Output is correct
60 Correct 187 ms 8648 KB Output is correct
61 Correct 203 ms 8648 KB Output is correct
62 Correct 211 ms 8648 KB Output is correct
63 Correct 185 ms 8608 KB Output is correct
64 Runtime error 74 ms 13368 KB Execution killed with signal 11
65 Halted 0 ms 0 KB -