結果
問題 | No.1171 Runs in Subsequences |
ユーザー | recososo |
提出日時 | 2020-08-14 23:03:41 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 57 ms / 2,000 ms |
コード長 | 1,757 bytes |
コンパイル時間 | 1,632 ms |
コンパイル使用メモリ | 171,732 KB |
実行使用メモリ | 17,276 KB |
最終ジャッジ日時 | 2023-07-31 23:43:11 |
合計ジャッジ時間 | 3,182 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge12 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 13 ms
15,132 KB |
testcase_01 | AC | 13 ms
15,232 KB |
testcase_02 | AC | 13 ms
15,124 KB |
testcase_03 | AC | 13 ms
15,096 KB |
testcase_04 | AC | 13 ms
15,164 KB |
testcase_05 | AC | 13 ms
15,160 KB |
testcase_06 | AC | 13 ms
15,124 KB |
testcase_07 | AC | 13 ms
15,060 KB |
testcase_08 | AC | 13 ms
15,116 KB |
testcase_09 | AC | 14 ms
15,040 KB |
testcase_10 | AC | 14 ms
15,168 KB |
testcase_11 | AC | 13 ms
15,132 KB |
testcase_12 | AC | 14 ms
15,120 KB |
testcase_13 | AC | 13 ms
15,136 KB |
testcase_14 | AC | 13 ms
15,136 KB |
testcase_15 | AC | 53 ms
16,132 KB |
testcase_16 | AC | 53 ms
15,888 KB |
testcase_17 | AC | 56 ms
16,144 KB |
testcase_18 | AC | 55 ms
16,096 KB |
testcase_19 | AC | 57 ms
17,276 KB |
testcase_20 | AC | 55 ms
16,524 KB |
testcase_21 | AC | 54 ms
16,208 KB |
ソースコード
#include<bits/stdc++.h> using namespace std; typedef long long ll; typedef long double ld; #define REP(i, n) for (int i = 0; i < (n); ++i) #define REPR(i, n) for (int i = n - 1; i >= 0; --i) #define FOR(i, m, n) for (int i = m; i < n; ++i) #define FORR(i, m, n) for (int i = m; i >= n; --i) #define ALL(v) (v).begin(),(v).end() template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; } template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return 1; } return 0; } const ll INF=1LL<<60; const int inf=(1<<30)-1; const int mod=1e9+7; int dx[8]={1,0,-1,0,-1,-1,1,1}; int dy[8]={0,1,0,-1,-1,1,-1,1}; const int nmax=500005; ll fac[nmax],finv[nmax],inv[nmax]; void COMinit(){ fac[0]=fac[1]=1; finv[0]=finv[1]=1; inv[1]=1; for(int i=2;i<nmax;i++){ fac[i]=fac[i-1]*i%mod; inv[i]=mod-inv[mod%i]*(mod/i)%mod; finv[i]=finv[i-1]*inv[i]%mod; } } ll com(int n,int k){ if(n<k||n<0||k<0){ return 0; } return fac[n]*(finv[k]*finv[n-k]%mod)%mod; } // 繰り返し二乗法 n*nが64bitにおさまる必要 ll pow_mod(ll n,ll k,ll m){ ll res=1; for(;k>0;k>>=1){ if(k&1){ res=(res*n)%m; } n=(n*n)%m; } return res; } int main(){ cin.tie(0); ios::sync_with_stdio(false); COMinit(); string s;cin >> s; ll n=s.size(); vector<vector<int>> c(26); REP(i,n){ c[s[i]-'a'].push_back(i); } ll ans=n*pow_mod(2,n-1,mod)%mod; REP(i,26){ int t=c[i].size(); vector<ll> a(t+1); REP(j,t){ a[j+1]=(a[j]+pow_mod(2,c[i][j],mod))%mod; } ll res=0; REP(j,t){ (res+=a[j]*pow_mod(2,n-1-c[i][j],mod)%mod)%=mod; } (ans+=mod-res)%=mod; } cout << ans << endl; }