結果

問題 No.1951 消えたAGCT(2)
ユーザー 蜜蜂
提出日時 2022-05-19 11:53:59
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 615 ms / 3,000 ms
コード長 1,797 bytes
コンパイル時間 3,422 ms
コンパイル使用メモリ 232,468 KB
実行使用メモリ 10,000 KB
最終ジャッジ日時 2024-09-17 17:37:02
合計ジャッジ時間 9,847 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 28
権限があれば一括ダウンロードができます

ソースコード

diff # 

//g++ 2.cpp -std=c++14 -O2 -I .
#include <bits/stdc++.h>
using namespace std;

#include <atcoder/all>
using namespace atcoder;

using ll = long long;
using ld = long double;
 
using vi = vector<int>;
using vvi = vector<vi>;
using vll = vector<ll>;
using vvll = vector<vll>;
using vld = vector<ld>;
using vvld = vector<vld>;
using vst = vector<string>;
using vvst = vector<vst>;
 
#define fi first
#define se second
#define pb push_back
#define eb emplace_back
#define pq_big(T) priority_queue<T,vector<T>,less<T>>
#define pq_small(T) priority_queue<T,vector<T>,greater<T>>
#define all(a) a.begin(),a.end()
#define rep(i,start,end) for(ll i=start;i<(ll)(end);i++)
#define per(i,start,end) for(ll i=start;i>=(ll)(end);i--)
#define uniq(a) sort(all(a));a.erase(unique(all(a)),a.end())

int op(int a,int b){
  return (int)(a+b);
}
int e(){
  return 0;
}

string agct = "AGCT";

int add=0;
vi cnt(26,0); // [i] := 元々'A'+i だったのがいくつあるか = 今は 'A'+i+add

int cntagct(){
  int res=0;
  for(char c:agct){
    int d=c-'A';
    res+=cnt[(d-add%26+26*100)%26];
  }
  return res;
}

int main(){
  ios::sync_with_stdio(false);
  cin.tie(nullptr);

  int n;
  string s;
  cin>>n>>s;

  segtree<int,op,e> seg(vi(n,1));

  for(char c:s){
    cnt[c-'A']++;
  }

  int ans=0;

  while(true){
    int c1=cntagct();
    if(!c1){
      cout<<ans<<endl;
      return 0;
    }
    ans++;

    int ng=-1,ok=n;
    while(ok-ng>1){
      int ch=(ng+ok)/2;
      if(seg.prod(0,ch+1)<c1){
        ng=ch;
      }
      else{
        ok=ch;
      }
    }

    assert(seg.get(ok)==1);
    assert(seg.prod(0,ok)==c1-1);
    assert(seg.prod(0,ok+1)==c1);

    // 元々ok番目にあったのを消す
    int c2=cnt[s[ok]-'A']-1;
    seg.set(ok,0);
    cnt[s[ok]-'A']--;

    add+=c2;
    add%=26;
  }
}
0