結果

問題 No.852 連続部分文字列
ユーザー LayCurseLayCurse
提出日時 2019-08-21 05:38:27
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 244 ms / 3,153 ms
コード長 1,782 bytes
コンパイル時間 2,487 ms
コンパイル使用メモリ 210,036 KB
最終ジャッジ日時 2025-01-07 14:26:50
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

diff # 

#pragma GCC optimize ("Ofast")
#include<bits/stdc++.h>
using namespace std;
inline void rd(char &c){
  int i;
  for(;;){
    i = getchar_unlocked();
    if(i!=' '&&i!='\n'&&i!='\r'&&i!='\t'&&i!=EOF){
      break;
    }
  }
  c = i;
}
inline int rd(char c[]){
  int i, sz=0;
  for(;;){
    i = getchar_unlocked();
    if(i!=' '&&i!='\n'&&i!='\r'&&i!='\t'&&i!=EOF){
      break;
    }
  }
  c[sz++] = i;
  for(;;){
    i = getchar_unlocked();
    if(i==' '||i=='\n'||i=='\r'||i=='\t'||i==EOF){
      break;
    }
    c[sz++] = i;
  }
  c[sz]='\0';
  return sz;
}
inline void wt_L(char a){
  putchar_unlocked(a);
}
inline void wt_L(double x){
  printf("%.15f",x);
}
char S[1000002];
int N;
int cnt[26];
int main(){
  double res;
  int i, j, k, m;
  N = rd(S);
  for(i=0;i<(N);i++){
    S[i] -= 'a';
  }
  res = 0;
  for(k=(1);k<(27);k++){
    j = 0;
    m = 0;
    for(i=0;i<(26);i++){
      cnt[i] = 0;
    }
    for(i=0;i<(N);i++){
      while(j<N && m<k){
        if((cnt[S[j]]++)==0){
          m++;
        }
        j++;
      }
      if(j==N && m<k){
        break;
      }
      res += N-j+1;
      if((--cnt[S[i]])==0){
        m--;
      }
    }
  }
  res /= (double)N*(N+1)/2;
  wt_L(res);
  wt_L('\n');
  return 0;
}
// cLay varsion 20190820-1

// --- original code ---
// char S[1000002]; int N;
// int cnt[26];
// {
//   int i, j, k, m;
//   double res;
//   
//   rd(S@N);
//   rep(i,N) S[i] -= 'a';
// 
//   res = 0;
// 
//   rep(k,1,27){
//     j = 0;
//     m = 0;
//     rep(i,26) cnt[i] = 0;
//     rep(i,N){
//       while(j<N && m<k){
//         if((cnt[S[j]]++)==0) m++;
//         j++;
//       }
//       if(j==N && m<k) break;
//       res += N-j+1;
//       if((--cnt[S[i]])==0) m--;
//     }
//   }
// 
//   res /= (double)N*(N+1)/2;
// 
//   wt(res);
// }
0