結果

問題 No.2911 位相の公理
ユーザー momoyuumomoyuu
提出日時 2024-10-06 15:34:14
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 67 ms / 2,000 ms
コード長 3,031 bytes
コンパイル時間 1,262 ms
コンパイル使用メモリ 108,900 KB
実行使用メモリ 6,120 KB
最終ジャッジ日時 2024-10-06 15:34:17
合計ジャッジ時間 2,634 ms
ジャッジサーバーID
(参考情報)
judge1 / judge3
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,248 KB
testcase_02 AC 2 ms
5,248 KB
testcase_03 AC 2 ms
5,248 KB
testcase_04 AC 1 ms
5,248 KB
testcase_05 AC 2 ms
5,248 KB
testcase_06 AC 2 ms
5,248 KB
testcase_07 AC 2 ms
5,248 KB
testcase_08 AC 2 ms
5,248 KB
testcase_09 AC 2 ms
5,248 KB
testcase_10 AC 2 ms
5,248 KB
testcase_11 AC 2 ms
5,248 KB
testcase_12 AC 2 ms
5,248 KB
testcase_13 AC 2 ms
5,248 KB
testcase_14 AC 2 ms
5,248 KB
testcase_15 AC 2 ms
5,248 KB
testcase_16 AC 2 ms
5,248 KB
testcase_17 AC 2 ms
5,248 KB
testcase_18 AC 2 ms
5,248 KB
testcase_19 AC 6 ms
5,248 KB
testcase_20 AC 67 ms
6,120 KB
testcase_21 AC 66 ms
6,028 KB
testcase_22 AC 5 ms
5,248 KB
testcase_23 AC 67 ms
6,112 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include<iostream>
#include<vector>
#include<algorithm>

using namespace std;
using ll = long long;

#include<atcoder/modint>
using m1 = atcoder::modint998244353;
using m2 = atcoder::modint1000000007;
#include<array>
const int B = 4;
using dat = array<ll,B>;
ll mods[] = {998244353,1000000007,1000000009,999999999};
dat add(dat a,dat b){
    dat e{};
    for(int i = 0;i<B;i++) e[i] = (a[i]+b[i]) % mods[i];
    return e;
}

dat minuss(dat a,dat b){
    for(int i = 0;i<B;i++){
        b[i] *= -1;
        if(b[i]<0) b[i] += mods[i];
    }
    return add(a,b);
}

dat Set(ll a){
    dat e{};
    for(int i = 0;i<B;i++) e[i] = a;
    return e;
}

ll mpow(ll a,ll k,ll mod){
    ll e = 1;
    while(k){
        if(k&1) e = (e*a) % mod;
        a = (a*a) % mod;
        k >>= 1;
    }
    return e;
}
dat make_dat(ll c){
    dat e{};
    for(int i = 0;i<B;i++){
        e[i] = mpow(2,c,mods[i]);
        e[i] -= 1;
        if(e[i]<0) e[i] += mods[i];
    }
    return e;
}
int main(){
    cin.tie(nullptr);
    ios::sync_with_stdio(false);

    int n,m;
    cin>>n>>m;
    vector<int> cnt(1<<n,0);
    for(int i = 0;i<m;i++){
        string s;
        cin>>s;
        int now = 0;
        for(int j = 0;j<n;j++) if(s[j]=='1') now |= 1 << j;
        cnt[now] = 1;
    }

    int mask = (1<<n) - 1;
    if(cnt[mask]==0){
        cout<<"No\n";
        return 0;
    }
    if(cnt[0]==0){
        cout<<"No\n";
        return 0;
    }

    {
        vector<ll> c(1<<n,0);
        for(int i = 0;i<1<<n;i++) c[i] = cnt[i];

        for(int i = 0;i<n;i++){
        for(int k = 0;k<1<<n;k++) if(k>>i&1){
                int nk = k ^ (1<<i);
                c[nk] += c[k];
            }
        }
        for(int i = 0;i<1<<n;i++) c[i] = ((c[i])*(c[i]-1)) / 2;

        for(int i = 0;i<n;i++){
        for(int k = 0;k<1<<n;k++)
            if(k>>i&1){
                int nk = k ^ (1<<i);
                c[nk] -= c[k];
            }
        }
        for(int i = 0;i<1<<n;i++) if(c[i]!=0){
            if(cnt[i]==0){
                cout<<"No\n";
                return 0;
            }
        }
    }

    {
        vector<ll> c(1<<n,0);
        for(int i = 0;i<1<<n;i++) c[i] = cnt[i];
        for(int i = 0;i<n;i++){
            for(int k = 0;k<1<<n;k++)
                if(k>>i&1){
                    int nk = k ^ (1<<i);
                    c[k] += c[nk];
                }
        }
        vector<dat> now(1<<n);
        for(int i = 0;i<1<<n;i++){
            now[i] = make_dat(c[i]);
        }
        for(int i = 0;i<n;i++) {
            for(int k = 0;k<1<<n;k++)
                if(k>>i&1){
                    int nk = k ^ (1<<i);
                    now[k] = minuss(now[k],now[nk]);
                }
        }

        //for(int i = 0;i<1<<n;i++) cout<<i<<" "<<now[i][0]<<endl;
        for(int i = 0;i<1<<n;i++){
            if(now[i]==Set(0)) continue;
            if(cnt[i]==0) {
                //cout<<i<<endl;
                cout<<"No\n";
                return 0;
            }
        }
    }
    cout<<"Yes\n";

}   


0