Submission #920726

#TimeUsernameProblemLanguageResultExecution timeMemory
920726zeta7532Genetics (BOI18_genetics)C++17
100 / 100
468 ms48468 KiB
#include <bits/stdc++.h>
#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
using namespace std;
using ll = long long;
const ll mod = 998244353;
#define fi first
#define se second
#define rep(i,n) for(ll i=0;i<n;i++)
#define all(x) x.begin(),x.end()
#define faster ios::sync_with_stdio(false);cin.tie(nullptr)
 
int main() {
    random_device seed_gen;
    mt19937 engine(seed_gen());
    uniform_real_distribution<float> dist(0,1);
    ll N,M,K;
    cin >> N >> M >> K;
    vector<string> s(N);
    rep(i,N){
        cin >> s[i];
    }
    vector<ll> P(N);
    ll sum=0;
    rep(i,N) P[i]=mod*dist(engine),sum+=P[i];
    vector<vector<ll>> cum(M,vector<ll>(4,0));
    rep(i,N){
        rep(j,M){
            if(s[i][j]=='A') cum[j][0]+=P[i];
            if(s[i][j]=='C') cum[j][1]+=P[i];
            if(s[i][j]=='G') cum[j][2]+=P[i];
            if(s[i][j]=='T') cum[j][3]+=P[i];
        }
    }
    rep(i,N){
        ll cnt=0;
        rep(j,M){
            if(s[i][j]!='A') cnt+=cum[j][0];
            if(s[i][j]!='C') cnt+=cum[j][1];
            if(s[i][j]!='G') cnt+=cum[j][2];
            if(s[i][j]!='T') cnt+=cum[j][3];
        }
        if(cnt%mod==((sum-P[i])*K)%mod){
            cout << i+1 << endl;
        }
    }
    return 0;
}
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...