結果
問題 | No.1171 Runs in Subsequences |
ユーザー | 👑 ygussany |
提出日時 | 2020-08-14 22:23:06 |
言語 | C (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 19 ms / 2,000 ms |
コード長 | 893 bytes |
コンパイル時間 | 1,742 ms |
コンパイル使用メモリ | 30,336 KB |
実行使用メモリ | 25,536 KB |
最終ジャッジ日時 | 2024-10-10 15:44:45 |
合計ジャッジ時間 | 2,175 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 18 |
ソースコード
#include <stdio.h> #define Mod 1000000007 long long div_mod(long long x, long long y, long long z) { if (x % y == 0) return x / y; else return (div_mod((1 + x / y) * y - x, (z % y), y) * z + x) / y; } int main() { char S[200001]; scanf("%s", S); int i, j, l, pos[26][200001] = {}, k[26] = {}; long long pow[200001], pow_inv[200001]; for (l = 0; S[l] != 0; l++) pos[S[l] - 'a'][k[S[l] - 'a']++] = l; for (i = 1, pow[0] = 1; i <= l; i++) pow[i] = pow[i-1] * 2 % Mod; for (i = l - 1, pow_inv[l] = div_mod(1, pow[l], Mod); i >= 0; i--) pow_inv[i] = pow_inv[i+1] * 2 % Mod; long long ans = l * pow[l-1] % Mod, tmp; for (i = 0; i < 26; i++) { for (j = 1, tmp = 0; j < k[i]; j++) { tmp = (tmp * pow_inv[pos[i][j] - pos[i][j-1]] % Mod + pow[l - (pos[i][j] - pos[i][j-1] + 1)]) % Mod; ans = (ans - tmp + Mod) % Mod; } } printf("%lld\n", ans); fflush(stdout); return 0; }