結果
| 問題 |
No.1848 Long Prefixes
|
| コンテスト | |
| ユーザー |
沙耶花
|
| 提出日時 | 2022-02-18 22:32:06 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 136 ms / 2,000 ms |
| コード長 | 2,992 bytes |
| コンパイル時間 | 5,257 ms |
| コンパイル使用メモリ | 258,088 KB |
| 最終ジャッジ日時 | 2025-01-28 00:22:16 |
|
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 40 |
ソースコード
#include <stdio.h>
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace atcoder;
using mint = modint1000000007;
using namespace std;
#define rep(i,n) for(int i=0;i<(n);i++)
#define Inf 1000000001
struct rolling_hash{
long long t_hash;
static vector<long long> power;
static const long long MOD = (1LL<<61)-1;
static const long long b = 123456;
int sz;
rolling_hash(){
sz = 0;
t_hash = 0;
}
rolling_hash(long long c){
sz = 1;
t_hash = b*c;
}
long long mul(__int128 x,__int128 y){
__int128 t = x*y;
t = (t>>61) + (t&MOD);
if(t>=MOD)t -= MOD;
return t;
}
long long get_pow(int sz){
if(power.size()>sz)return power[sz];
while(power.size()<=sz){
if(power.size()==0)power.push_back(1);
else power.push_back(mul(power.back(),b));
}
return power.back();
}
rolling_hash &operator+=(const rolling_hash &another){
(*this).t_hash = mul((*this).t_hash,get_pow(another.sz));
(*this).t_hash += another.t_hash;
if((*this).t_hash>=MOD)(*this).t_hash -= MOD;
(*this).sz += another.sz;
return (*this);
}
rolling_hash operator+(const rolling_hash &another)const{
return (rolling_hash(*this)+=another);
}
rolling_hash &operator-=(const rolling_hash &another){
(*this).t_hash += MOD - mul(another.t_hash,get_pow((*this).sz-another.sz));
if((*this).t_hash>=MOD)(*this).t_hash -= MOD;
(*this).sz -= another.sz;
return (*this);
}
rolling_hash operator-(const rolling_hash &another)const{
return (rolling_hash(*this)-=another);
}
bool operator<(const rolling_hash &another)const{
if((*this).t_hash!=another.t_hash)return (*this).t_hash<another.t_hash;
return (*this).sz<another.sz;
}
bool operator==(const rolling_hash &another)const{
return ((*this).t_hash==another.t_hash && (*this).sz==another.sz);
}
};
vector<long long> rolling_hash::power;
mint sum(long long n){
mint ret = n;
ret *= n+1;
ret /= 2;
return ret;
}
int main(){
int n;
cin>>n;
vector<long long> a(n);
rep(i,n)cin>>a[i];
string s;
cin>>s;
{
vector<long long> b;
string t = "";
rep(i,n){
if(t.size()==0 || t.back()!=s[i]){
b.push_back(a[i]);
t += s[i];
}
else{
b.back() += a[i];
}
}
a = b;
s = t;
n = a.size();
}
vector<rolling_hash> R(n+1);
rep(i,n){
long long t = s[i];
t *= a[i];
R[i+1] = rolling_hash(t);
R[i+1] = R[i] + R[i+1];
}
vector<long long> SS(n+1,0);
rep(i,n){
SS[i+1] = SS[i] + a[i];
}
mint ans = 0;
rep(i,n){
if(s[i] != s[0])continue;
if(a[i] < a[0]){
ans += sum(a[i]);
}
else{
ans += sum(a[0]);
ans -= a[0];
ans += a[0] * (a[i] - a[0]);
int ok = 1,ng = n+1;
while(ng-ok>1){
int mid = (ok+ng)/2;
if(i+mid>n)ng = mid;
else if(R[mid] - R[1] == R[i+mid] - R[i+1])ok = mid;
else ng= mid;
}
mint temp = a[0] + SS[ok] - SS[1];
if(i+ok<n){
if(s[i+ok]==s[ok])ans += min(a[i+ok],a[ok]);
}
ans += temp;
}
}
cout<<ans.val()<<endl;
return 0;
}
沙耶花