Submission #546453

#	Time	Username	Problem	Language	Result	Execution time	Memory
546453		xynex	Mutating DNA (IOI21_dna)	C++17	100 / 100	57 ms	7432 KiB

dna

/**
  * Author: Tenjin
  * Created: 06.04.2022 19:13
  * Why am I so stupid? :c
  * Slishkom slab
**/
#include <bits/stdc++.h>

// #pragma GCC optimize("inline")
// #pragma GCC optimize("-fgcse,-fgcse-lm")
// #pragma GCC optimize("-ftree-pre,-ftree-vrp")
// #pragma GCC optimize("-ffast-math")
// #pragma GCC optimize("-fipa-sra")
// #pragma GCC optimize("-fpeephole2")
// #pragma GCC optimize("-fsched-spec")
// #pragma GCC optimize("Ofast,no-stack-protector")
// #pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,avx2,tune=native")
// #pragma GCC optimize("unroll-loops")

using namespace std;

#define ll long long
#define dl double long
#define ull unsigned long long
#define pr pair
#define vt vector
#define ff first
#define ss second
#define mp make_pair
#define sz(a) (int)a.size()
#define pb push_back
#define pf push_front
#define popB pop_back
#define popF pop_front
#define bit_count __builtin_popcount
#define all(x) x.begin(), x.end()
#define rall(x) x.rbegin(), x.rend()
#define sp(x) fixed << setprecision(x)

template<typename T> T get_rand(T l, T r) {
  random_device rd;
  mt19937 gen(rd());
  return uniform_int_distribution<T>(l, r)(gen);
}
template<typename T> T lcm(T a, T b) { return a * (b / __gcd(a, b)); }

template<class A> void read(vt<A>& v);
template<class A, size_t S> void read(array<A, S>& a);
template<class T> void read(T& x) { cin >> x; }
void read(double& d) { string t; read(t); d = stod(t); }
void read(long double& d) { string t; read(t); d = stold(t); }
template<class H, class... T> void read(H& h, T&... t) { read(h); read(t...); }
template<class A> void read(vt<A>& x) { for (auto& a : x) read(a); }
template<class A, size_t S> void read(array<A, S>& x) { for (auto& a : x) read(a); }

string to_string(char c) { return string(1, c); }
string to_string(bool b) { return b ? "true" : "false"; }
string to_string(const char* s) { return string(s); }
string to_string(string s) { return s; }
string to_string(vt<bool> v) { string res; for (int i = 0; i < sz(v); ++i) res += char('0' + v[i]); return res; }
template<size_t S> string to_string(bitset<S> b) { string res; for (int i = 0; i < S; ++i) res += char('0' + b[i]); return res; }
template<class T> string to_string(T v) { bool f = 1; string res; for (auto x : v) { if (!f) res += ' '; f = 0; res += to_string(x); } return res; }

template<class A> void write(A x) { cout << to_string(x); }
template<class H, class... T> void write(const H& h, const T&... t) { write(h); write(t...); }

void print() { write("\n"); }
template<class H, class... T> void print(const H& h, const T&... t) { write(h); if (sizeof...(t)) write(' '); print(t...); }

void freop(string s) {
  freopen((s + ".in").c_str(), "r", stdin);
  freopen((s + ".out").c_str(), "w", stdout);
}

const int MOD = 1e9 + 7;
const int N = 3e6 + 5;
const ll INF = 9e18;
const int M = 3e3 + 5;
const dl pi = acos(-1);
const dl eps = 1e-12;
const int sq = 700;

int dx[] = {-1, 0, 1, 0};
int dy[] = {0, -1, 0, 1};
/* ll vs int*/
/*
ATACAT ACTATA
3
1 3
4 5
3 5
*/
int pref[N][3][3];
int n;
void init(string a, string b) {
  n = sz(a);
  map<char, int> whr;
  whr['A'] = 0, whr['T'] = 1, whr['C'] = 2;
  for(int i = 1; i <= n; ++i) {
    for(int j = 0; j < 3; ++j) {
      for(int k = 0; k < 3; ++k) {
        pref[i][j][k] = pref[i - 1][j][k];
      }
    }
    pref[i][whr[a[i - 1]]][whr[b[i - 1]]]++;
  }
}

int cur[3][3];
int get_distance(int x, int y) {
  for(int i = 0; i < 3; ++i) {
    for(int j = 0; j < 3; ++j) {
      cur[i][j] = pref[y + 1][i][j] - pref[x][i][j];
    //  print(i, j, cur[i][j]);
    }
  }
  int ans = 0;
  for(int i = 0; i < 3; ++i) {
    for(int j = i + 1; j < 3; ++j) {
      int mn = min(cur[i][j], cur[j][i]);
      ans += mn;
      cur[i][j] -= mn;
      cur[j][i] -= mn;
    }
  }
  vt<int> ord(3);
  for(int i = 0; i < 3; ++i) {
    for(int j = 0; j < 3; ++j) {
      if(i == j) continue;
      ord[i] += cur[i][j];
    }
  }
  //print(ord);
  if((ord[0] + ord[1] + ord[2]) % 3) return -1;
  for(int i = 0; i < 3; ++i) {
    for(int j = 0; j < 3; ++j) {
      if(ord[i] != ord[j]) return -1;
    }
  }
  return ans + (ord[0] + ord[1] + ord[2]) / 3 * 2;
}

/*int main() {
  string a, b; read(a, b);
  init(a, b);
  int q; read(q);
  while(q--) {
    int x, y; read(x, y);
    print(get_distance(x, y));
  }
}*/

Compilation message (stderr)

dna.cpp: In function 'void freop(std::string)':
dna.cpp:71:10: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
   71 |   freopen((s + ".in").c_str(), "r", stdin);
      |   ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
dna.cpp:72:10: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
   72 |   freopen((s + ".out").c_str(), "w", stdout);
      |   ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

• Subtask 1: 21 / 21
• Subtask 2: 22 / 22
• Subtask 3: 13 / 13
• Subtask 4: 28 / 28
• Subtask 5: 16 / 16