結果
| 問題 | No.3410 Happiest Art |
| コンテスト | |
| ユーザー |
 C
|
| 提出日時 | 2026-11-04 00:49:46 |
| 言語 | C++23(gcc16) (gcc 16.2.0 + boost 1.92.0 + ACL) |
| 結果 |
AC
不安定
|
| 実行時間 | 57 ms / 4,000 ms |
| + 379µs | |
| コード長 | 13,809 bytes |
| 記録 | |
| コンパイル時間 | 5,842 ms |
| コンパイル使用メモリ | 404,044 KB |
| 実行使用メモリ | 54,548 KB |
| 最終ジャッジ日時 | 2026-11-04 00:52:13 |
| 合計ジャッジ時間 | 10,436 ms |
|
ジャッジサーバーID （参考情報） |
judge1_0 / judge7_0 |
（要ログイン）
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 48 |
コンパイルメッセージ
hashmap/hashmap-base.hpp:15:7: warning: 'template<class _Category, class _Tp, class _Distance, class _Pointer, class _Reference> struct std::iterator' is deprecated [-Wdeprecated-declarations]
In file included from /home/linuxbrew/.linuxbrew/Cellar/gcc/16.2.0/include/c++/16/bits/stl_algobase.h:64,
from /home/linuxbrew/.linuxbrew/Cellar/gcc/16.2.0/include/c++/16/algorithm:62,
from /home/linuxbrew/.linuxbrew/Cellar/gcc/16.2.0/include/c++/16/x86_64-pc-linux-gnu/bits/stdc++.h:51,
from main.cpp:1:
/home/linuxbrew/.linuxbrew/Cellar/gcc/16.2.0/include/c++/16/bits/stl_iterator_base_types.h:129:34: note: declared here
129 | struct _GLIBCXX17_DEPRECATED iterator
| ^~~~~~~~
hashmap/hashmap-base.hpp:17:7: warning: 'template<class _Category, class _Tp, class _Distance, class _Pointer, class _Reference> struct std::iterator' is deprecated [-Wdeprecated-declarations]
/home/linuxbrew/.linuxbrew/Cellar/gcc/16.2.0/include/c++/16/bits/stl_iterator_base_types.h:129:34: note: declared here
129 | struct _GLIBCXX17_DEPRECATED iterator
| ^~~~~~~~
ソースコード
#include <bits/stdc++.h>
#include <atcoder/all>
#include <boost/polygon/point_data.hpp>
#include <boost/polygon/voronoi.hpp>
#define rep(i,a,b) for(int i=(a);i<(b);i++)
#define rrep(i,a,b) for(int i=(b)-1;i>=(a);i--)
#define all(x) begin(x),end(x)
#define sz(x) (int)(x).size()
using namespace std;
using namespace atcoder;
using ll=long long;
using ld=long double;
using vll=vector<ll>;
using vvll=vector<vll>;
using vi=vector<int>;
using vvi=vector<vi>;
using pll=pair<ll,ll>;
using pi=pair<int,int>;
// using mint=modint998244353;
using mint=modint1000000007;
// using mint2=static_modint<2>;
// using BoostPoint=boost::polygon::point_data<int>;
// using Voronoi=boost::polygon::voronoi_diagram<double>;
// using mint=LazyMontgomeryModInt<998244353>;
using vm=vector<mint>;
// using vvm=vector<vm>;
#line 2 "hashmap/hashmap.hpp"
#line 2 "hashmap/hashmap-base.hpp"
#include <cstdint>
using namespace std;
namespace HashMapImpl {
using u32 = uint32_t;
using u64 = uint64_t;
template <typename Key, typename Data>
struct HashMapBase;
template <typename Key, typename Data>
struct itrB
: iterator<bidirectional_iterator_tag, Data, ptrdiff_t, Data*, Data&> {
using base =
iterator<bidirectional_iterator_tag, Data, ptrdiff_t, Data*, Data&>;
using ptr = typename base::pointer;
using ref = typename base::reference;
u32 i;
HashMapBase<Key, Data>* p;
explicit constexpr itrB() : i(0), p(nullptr) {}
explicit constexpr itrB(u32 _i, HashMapBase<Key, Data>* _p) : i(_i), p(_p) {}
explicit constexpr itrB(u32 _i, const HashMapBase<Key, Data>* _p)
: i(_i), p(const_cast<HashMapBase<Key, Data>*>(_p)) {}
friend void swap(itrB& l, itrB& r) { swap(l.i, r.i), swap(l.p, r.p); }
friend bool operator==(const itrB& l, const itrB& r) { return l.i == r.i; }
friend bool operator!=(const itrB& l, const itrB& r) { return l.i != r.i; }
const ref operator*() const {
return const_cast<const HashMapBase<Key, Data>*>(p)->data[i];
}
ref operator*() { return p->data[i]; }
ptr operator->() const { return &(p->data[i]); }
itrB& operator++() {
assert(i != p->cap && "itr::operator++()");
do {
i++;
if (i == p->cap) break;
if (p->occupied_flag[i] && !p->deleted_flag[i]) break;
} while (true);
return (*this);
}
itrB operator++(int) {
itrB it(*this);
++(*this);
return it;
}
itrB& operator--() {
do {
i--;
if (p->occupied_flag[i] && !p->deleted_flag[i]) break;
assert(i != 0 && "itr::operator--()");
} while (true);
return (*this);
}
itrB operator--(int) {
itrB it(*this);
--(*this);
return it;
}
};
template <typename Key, typename Data>
struct HashMapBase {
using u32 = uint32_t;
using u64 = uint64_t;
using iterator = itrB<Key, Data>;
using itr = iterator;
protected:
template <typename K>
inline u64 randomized(const K& key) const {
return u64(key) ^ r;
}
template <typename K,
enable_if_t<is_same<K, Key>::value, nullptr_t> = nullptr,
enable_if_t<is_integral<K>::value, nullptr_t> = nullptr>
inline u32 inner_hash(const K& key) const {
return (randomized(key) * 11995408973635179863ULL) >> shift;
}
template <
typename K, enable_if_t<is_same<K, Key>::value, nullptr_t> = nullptr,
enable_if_t<is_integral<decltype(K::first)>::value, nullptr_t> = nullptr,
enable_if_t<is_integral<decltype(K::second)>::value, nullptr_t> = nullptr>
inline u32 inner_hash(const K& key) const {
u64 a = randomized(key.first), b = randomized(key.second);
a *= 11995408973635179863ULL;
b *= 10150724397891781847ULL;
return (a + b) >> shift;
}
template <typename K,
enable_if_t<is_same<K, Key>::value, nullptr_t> = nullptr,
enable_if_t<is_integral<typename K::value_type>::value, nullptr_t> =
nullptr>
inline u32 inner_hash(const K& key) const {
static constexpr u64 mod = (1LL << 61) - 1;
static constexpr u64 base = 950699498548472943ULL;
u64 res = 0;
for (auto& elem : key) {
__uint128_t x = __uint128_t(res) * base + (randomized(elem) & mod);
res = (x & mod) + (x >> 61);
}
__uint128_t x = __uint128_t(res) * base;
res = (x & mod) + (x >> 61);
if (res >= mod) res -= mod;
return res >> (shift - 3);
}
template <typename D = Data,
enable_if_t<is_same<D, Key>::value, nullptr_t> = nullptr>
inline u32 hash(const D& dat) const {
return inner_hash(dat);
}
template <
typename D = Data,
enable_if_t<is_same<decltype(D::first), Key>::value, nullptr_t> = nullptr>
inline u32 hash(const D& dat) const {
return inner_hash(dat.first);
}
template <typename D = Data,
enable_if_t<is_same<D, Key>::value, nullptr_t> = nullptr>
inline Key data_to_key(const D& dat) const {
return dat;
}
template <
typename D = Data,
enable_if_t<is_same<decltype(D::first), Key>::value, nullptr_t> = nullptr>
inline Key data_to_key(const D& dat) const {
return dat.first;
}
void reallocate(u32 ncap) {
vector<Data> ndata(ncap);
vector<bool> nf(ncap);
shift = 64 - __lg(ncap);
for (u32 i = 0; i < cap; i++) {
if (occupied_flag[i] && !deleted_flag[i]) {
u32 h = hash(data[i]);
while (nf[h]) h = (h + 1) & (ncap - 1);
ndata[h] = move(data[i]);
nf[h] = true;
}
}
data.swap(ndata);
occupied_flag.swap(nf);
cap = ncap;
occupied = s;
deleted_flag.resize(cap);
fill(std::begin(deleted_flag), std::end(deleted_flag), false);
}
inline bool extend_rate(u32 x) const { return x * 2 >= cap; }
inline bool shrink_rate(u32 x) const {
return HASHMAP_DEFAULT_SIZE < cap && x * 10 <= cap;
}
inline void extend() { reallocate(cap << 1); }
inline void shrink() { reallocate(cap >> 1); }
public:
u32 cap, s, occupied;
vector<Data> data;
vector<bool> occupied_flag, deleted_flag;
u32 shift;
static u64 r;
static constexpr uint32_t HASHMAP_DEFAULT_SIZE = 4;
explicit HashMapBase()
: cap(HASHMAP_DEFAULT_SIZE),
s(0),
occupied(0),
data(cap),
occupied_flag(cap),
deleted_flag(cap),
shift(64 - __lg(cap)) {}
itr begin() const {
u32 h = 0;
while (h != cap) {
if (occupied_flag[h] && !deleted_flag[h]) break;
h++;
}
return itr(h, this);
}
itr end() const { return itr(this->cap, this); }
friend itr begin(const HashMapBase& h) { return h.begin(); }
friend itr end(const HashMapBase& h) { return h.end(); }
itr find(const Key& key) const {
u32 h = inner_hash(key);
while (true) {
if (occupied_flag[h] == false) return this->end();
if (data_to_key(data[h]) == key) {
if (deleted_flag[h] == true) return this->end();
return itr(h, this);
}
h = (h + 1) & (cap - 1);
}
}
bool contain(const Key& key) const { return find(key) != this->end(); }
itr insert(const Data& d) {
u32 h = hash(d);
while (true) {
if (occupied_flag[h] == false) {
if (extend_rate(occupied + 1)) {
extend();
h = hash(d);
continue;
}
data[h] = d;
occupied_flag[h] = true;
++occupied, ++s;
return itr(h, this);
}
if (data_to_key(data[h]) == data_to_key(d)) {
if (deleted_flag[h] == true) {
data[h] = d;
deleted_flag[h] = false;
++s;
}
return itr(h, this);
}
h = (h + 1) & (cap - 1);
}
}
// tips for speed up :
// if return value is unnecessary, make argument_2 false.
itr erase(itr it, bool get_next = true) {
if (it == this->end()) return this->end();
s--;
if (!get_next) {
this->deleted_flag[it.i] = true;
if (shrink_rate(s)) shrink();
return this->end();
}
itr nxt = it;
nxt++;
this->deleted_flag[it.i] = true;
if (shrink_rate(s)) {
Data d = data[nxt.i];
shrink();
it = find(data_to_key(d));
}
return nxt;
}
itr erase(const Key& key) { return erase(find(key)); }
int count(const Key& key) { return find(key) == end() ? 0 : 1; }
bool empty() const { return s == 0; }
int size() const { return s; }
void clear() {
fill(std::begin(occupied_flag), std::end(occupied_flag), false);
fill(std::begin(deleted_flag), std::end(deleted_flag), false);
s = occupied = 0;
}
void reserve(int n) {
if (n <= 0) return;
n = 1 << min(23, __lg(n) + 2);
if (cap < u32(n)) reallocate(n);
}
};
template <typename Key, typename Data>
uint64_t HashMapBase<Key, Data>::r =
chrono::duration_cast<chrono::nanoseconds>(
chrono::high_resolution_clock::now().time_since_epoch())
.count();
} // namespace HashMapImpl
/**
* @brief Hash Map(base) (ハッシュマップ・基底クラス)
*/
#line 4 "hashmap/hashmap.hpp"
template <typename Key, typename Val>
struct HashMap : HashMapImpl::HashMapBase<Key, pair<Key, Val>> {
using base = typename HashMapImpl::HashMapBase<Key, pair<Key, Val>>;
using HashMapImpl::HashMapBase<Key, pair<Key, Val>>::HashMapBase;
using Data = pair<Key, Val>;
Val& operator[](const Key& k) {
typename base::u32 h = base::inner_hash(k);
while (true) {
if (base::occupied_flag[h] == false) {
if (base::extend_rate(base::occupied + 1)) {
base::extend();
h = base::hash(k);
continue;
}
base::data[h].first = k;
base::data[h].second = Val();
base::occupied_flag[h] = true;
++base::occupied, ++base::s;
return base::data[h].second;
}
if (base::data[h].first == k) {
if (base::deleted_flag[h] == true) {
base::data[h].second = Val();
base::deleted_flag[h] = false;
++base::s;
}
return base::data[h].second;
}
h = (h + 1) & (base::cap - 1);
}
}
typename base::itr emplace(const Key& key, const Val& val) {
return base::insert(Data(key, val));
}
};
/*
* @brief ハッシュマップ(連想配列)
* @docs docs/hashmap/hashmap.md
**/
struct node{
int score=0,i=-1,k=-1;
};
int main()
{
ios::sync_with_stdio(false);
cin.tie(nullptr);
/**
* O(NM) ,M=H*W
*
* M<=10 の場合：
* 考えられるグリッドは高々 2^10 個なので全探索する
* ここはO(2^10*NM)
*
* M>=11、正の幸福度の場合：
* ハミング距離を2つのグリッド間で異なるマスの数とし
* S[i] から距離 F[i] 以下のグリッドの集合をハミング球と考える
*
* 幸福度が正なら、少なくとも1人の好きな人を満足させる必要があるため、
* 好きな人のハミング球の和集合だけを候補にすれば良い
*
* まず好きな人のハミング球を列挙し、候補をハッシュ表に登録する
* 次に全員のハミング球を列挙し、表に存在する候補に重みを加算する
* このように両側から列挙してハッシュで照合することで、
* 球同士を比較せずに、共通するグリッドを直接見つけられる
*
* 各球は高々 M+1 個のグリッドしか含まず、
* Zobrist hash により1マス反転後のハッシュも O(1) で求まる
* ここは O(NM)
*
* M>=11、幸福度が０の場合：
* 全員と2マス以上異なるグリッドが必ず存在し、幸福度0を達成できる
* 最初の11マスだけ考えると、各 S[i] が禁止するパターンは
* 自身とその1ビット反転の12通りだけ
* 12N<=1200<2048=2^11 ので、禁止されないパターンが必ず存在する
* 禁止パターンをすべて記録し、残ったものを先頭11マスに採用する
* それ以降のマスは任意で良い
* ここはO(2^11+N+M)
*/
int N,U,H,W;
cin>>N>>U>>H>>W;
int M=H*W;
vi F(N),w(N);
vector<string>S(N);
rep(i,0,N){
cin>>F[i];
rep(j,0,H){
string row;
cin>>row;
S[i]+=row;
}
w[i]=(i<U)?1:-1;
}
//M<=10 の場合：
if(M<=10){
int mx=INT_MIN,argmx=0;
rep(mask,0,1<<M){
int score=0;
rep(i,0,N){
int d=0;
rep(k,0,M){
char c=((mask>>k)&1)?'#':'.';
if(c!=S[i][k])d++;
if(d>F[i])break;
}
if(d<=F[i])score+=w[i];
}
if(score>mx){
mx=score;
argmx=mask;
}
}
cout<<mx<<'\n';
rep(i,0,H){
rep(j,0,W){
int k=i*W+j;
cout<<(((argmx>>k)&1)?'#':'.');
}
cout<<'\n';
}
return 0;
}
// M>=11、正の幸福度の場合：
mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count());
#define ull unsigned long long
vector<ull>z(M),h(N);
rep(i,0,M)z[i]=rng();
rep(i,0,N)rep(j,0,M)if(S[i][j]=='#')h[i]^=z[j];
HashMap<ull,node>mp;
int lim=U;
rep(i,0,U)lim+=F[i]*M;
mp.reserve(lim);
rep(i,0,U){
mp.emplace(h[i],{0,i,-1});
if(F[i]){
rep(k,0,M)mp.emplace(h[i]^z[k],{0,i,k});
}
}
rep(i,0,N){
auto it=mp.find(h[i]);
if(it!=mp.end())it->second.score+=w[i];
if(F[i]){
rep(k,0,M){
it=mp.find(h[i]^z[k]);
if(it!=mp.end())it->second.score+=w[i];
}
}
}
int mx=INT_MIN;
pi argmx;
for(auto& [key,v]:mp){
if(v.score>mx){
mx=v.score;
argmx={v.i,v.k};
}
}
if(mx>=0){
cout<<mx<<'\n';
auto& ans=S[argmx.first];
int k=argmx.second;
if(k!=-1)ans[k]=(ans[k]=='.')?'#':'.';
rep(j,0,H)cout<<ans.substr(j*W,W)<<'\n';
return 0;
}
//M>=11、幸福度が０の場合：
vector<bool> bad(1<<11);
rep(i,0,N){
int mask=0;
rep(k,0,11)if(S[i][k]=='#')mask|=(1<<k);
bad[mask]=true;
rep(k,0,11)bad[mask^(1<<k)]=true;
}
rep(mask,0,1<<11){
if(bad[mask])continue;
cout<<0<<'\n';
string ans(M,'.');
rep(k,0,11)if((mask>>k)&1)ans[k]='#';
rep(j,0,H)cout<<ans.substr(j*W,W)<<'\n';
return 0;
}
}