-
Notifications
You must be signed in to change notification settings - Fork 4
/
Copy pathE.cpp
40 lines (39 loc) · 1.23 KB
/
E.cpp
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
#include <bits/stdc++.h>
#include <tr2/dynamic_bitset>
using namespace std;
using namespace tr2;
int main() {
ios_base::sync_with_stdio(false);
cin.tie(nullptr);
int n; cin >> n;
string s, t; cin >> s >> t;
vector<int64_t> a(1 << n), b(1 << n);
for (int msk = 0; msk < (1 << n); msk++) {
a[msk] = (1LL * (s[msk] - '0')) << (2 * __builtin_popcount(msk));
b[msk] = (1LL * (t[msk] - '0')) << (2 * __builtin_popcount(msk));
}
for (int bit = 0; bit < n; bit++) {
for (int msk = 0; msk < (1 << n); msk++) {
if (msk >> bit & 1) {
a[msk] += a[msk ^ (1 << bit)];
b[msk] += b[msk ^ (1 << bit)];
}
}
}
vector<int64_t> ans(1 << n);
for (int msk = 0; msk < (1 << n); msk++) {
ans[msk] += a[msk] * b[msk];
}
for (int bit = n - 1; bit >= 0; bit--) {
for (int msk = 0; msk < (1 << n); msk++) {
if (msk >> bit & 1) {
ans[msk] -= ans[msk ^ (1 << bit)];
}
}
}
string ret;
for (int msk = 0; msk < (1 << n); msk++) {
ret += (char)((ans[msk] >> ((2 * __builtin_popcount(msk))) & 3) + '0');
}
cout << ret;
}