結果
| 問題 |
No.1171 Runs in Subsequences
|
| コンテスト | |
| ユーザー |
👑 |
| 提出日時 | 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;
}