結果
問題 | No.1848 Long Prefixes |
ユーザー |
|
提出日時 | 2022-07-03 23:33:03 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 107 ms / 2,000 ms |
コード長 | 3,614 bytes |
コンパイル時間 | 4,158 ms |
コンパイル使用メモリ | 232,580 KB |
実行使用メモリ | 9,088 KB |
最終ジャッジ日時 | 2024-11-29 23:50:28 |
合計ジャッジ時間 | 8,185 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 40 |
ソースコード
#include<bits/stdc++.h>#include<atcoder/all>#define rep(i,b) for(int i=0;i<b;i++)#define rrep(i,b) for(int i=b-1;i>=0;i--)#define rep1(i,b) for(int i=1;i<b;i++)#define repx(i,x,b) for(int i=x;i<b;i++)#define fore(i,a) for(auto i:a)#define fore1(i,a) for(auto &i:a)#define rng(x) (x).begin(), (x).end()#define rrng(x) (x).rbegin(), (x).rend()#define sz(x) ((int)(x).size())#define pb push_back#define fi first#define se secondusing namespace std;using namespace atcoder;using ll = long long;template<typename T> using mpq = priority_queue<T, vector<T>, greater<T>>;template<typename T> bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; }template<typename T> bool chmin(T &a, const T &b) { if (b<a) { a=b; return 1; } return 0; }template<typename T> ll sumv(const vector<T>&a){ll res(0);for(auto&&x:a)res+=x;return res;}bool yn(bool a) { if(a) {cout << "Yes" << endl; return 1;} else {cout << "No" << endl; return 0;}}#define dame { cout << "No" << endl; return 0;}#define dame1 { cout << -1 << endl; return 0;}#define test(x) cout << "test" << x << endl;#define deb(x,y) cout << x << " " << y << endl;#define deb3(x,y,z) cout << x << " " << y << " " << z << endl;#define deb4(x,y,z,x2) cout << x << " " << y << " " << z << " " << x2 << endl;#define out cout << ans << endl;#define outv fore(yans , ans) cout << yans << "\n";#define show(x) cerr<<#x<<" = "<<x<<endl;#define showv(v) {fore(vy , v) {cout << vy << " ";} cout << endl;}#define showv2(v) fore(vy , v) cout << vy << "\n";using pll = pair<ll,ll>;using pil = pair<int,ll>;using pli = pair<ll,int>;using pii = pair<int,int>;using tp = tuple<int ,int ,int>;using vi = vector<int>;using vl = vector<ll>;using vs = vector<string>;using vb = vector<bool>;using vpii = vector<pii>;using vpli = vector<pli>;using vpll = vector<pll>;using vpil = vector<pil>;using vvi = vector<vector<int>>;using vvl = vector<vector<ll>>;using vvs = vector<vector<string>>;using vvb = vector<vector<bool>>;using vvpii = vector<vector<pii>>;using vvpli = vector<vector<pli>>;using vvpll = vector<vpll>;using vvpil = vector<vpil>;using mint = modint1000000007;using vm = vector<mint>;using vvm = vector<vector<mint>>;vector<int> dx={1,0,-1,0,1,1,-1,-1},dy={0,1,0,-1,1,-1,1,-1};ll gcd(ll a, ll b) { return a?gcd(b%a,a):b;}ll lcm(ll a, ll b) { return a/gcd(a,b)*b;}const double eps = 1e-10;const ll LINF = 1001002003004005006ll;const int INF = 1001001001;int main(){ll n; cin>>n;vl a(n-1),b(n);fenwick_tree<mint> fw(n-1);rep(i,n){cin>>b[i];if (i != 0){a[i-1] = b[i];fw.add(i-1 , b[i]);}}string s; cin>>s;vi sv = z_algorithm(s);vi av = z_algorithm<ll>(a);mint ans = 0;rep(i,n){// show(i);// cout << ans.val() << endl;if (sv[i] == 0) continue;if (b[0] > b[i]){ans += mint(b[i]) * mint(b[i]+1) / 2;continue;}ll k = b[i] - b[0] + 1;ans += mint(b[0]) * mint(b[0]-1) / 2;ans += mint(b[0]) * mint(k);// show(i);// cout << ans.val() << endl;if (i == n-1) continue;int r = min(av[i] , sv[i]-1);ans += fw.sum(0 , r);// show(i);// cout << ans.val() << endl;if (i+av[i]+1 >= n) continue;if (sv[i] > av[i] + 1){ans += min(b[i+av[i]+1] , b[av[i]+1]);// show(i);}// show(i);// cout << ans.val() << endl;}cout << ans.val() << endl;return 0;}