Morass

Limak

Mar 1st, 2016
113
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 1.32 KB | None | 0 0
  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. #define PB push_back
  4. #define ZERO (1e-10)
  5. #define INF (1<<29)
  6. #define CL(A,I) (memset(A,I,sizeof(A)))
  7. #define DEB printf("DEB!\n");
  8. #define D(X) cout<<"  "<<#X": "<<X<<endl;
  9. #define EQ(A,B) (A+ZERO>B&&A-ZERO<B)
  10. typedef long long ll;
  11. typedef long double ld;
  12. typedef pair<ll,ll> pll;
  13. typedef vector<int> vi;
  14. typedef pair<int,int> ii;
  15. typedef vector<ii> vii;
  16. #define IN(n) int n;scanf("%d",&n);
  17. #define FOR(i, m, n) for (int i(m); i < n; i++)
  18. #define REP(i, n) FOR(i, 0, n)
  19. #define F(n) REP(i, n)
  20. #define FF(n) REP(j, n)
  21. #define FT(m, n) FOR(k, m, n)
  22. #define aa first
  23. #define bb second
  24. void ga(int N,int *A){F(N)scanf("%d",A+i);}
  25. #define MX (1<<19)
  26. #define IX(C) (C=='A'?0:C=='C'?1:C=='G'?2:3)
  27. char s[MX];
  28. int N,A[MX][4],q[4],O;
  29. void gt(int b,int e,int*q){
  30.     F(4)q[i]=A[N][i];
  31.     F(4)q[i]-=A[e][i];
  32.     F(4)q[i]+=A[b-1][i];
  33. }
  34. bool ok(int d){
  35.     F(N-d+1){
  36.         int S(0);
  37.         gt(i+1,i+d,q);
  38.         F(4)S+=abs(O-q[i]);
  39.         if(S<=d)return 1;
  40.     }
  41.     return 0;
  42. }
  43. int main(void){
  44.     scanf("%d%s",&N,s+1);
  45.     FT(1,N+1)memcpy(A[k],A[k-1],sizeof(A[k])),++A[k][IX(s[k])];
  46.     O=N>>2;
  47.     int B(0),E(N),M;
  48.     gt(1,5,q);
  49.     while(B+1<E){
  50.         if(!ok(M=(B+E+1)/2))B=M+1;
  51.         else E=M;
  52.     }
  53.     printf("%d\n",B+!ok(B));
  54.     return 0;
  55. }
Advertisement
Add Comment
Please, Sign In to add comment