#include "bits/stdc++.h" #include #define rep(i,n) for(ll i = 0; i < n; i++) typedef long long ll; typedef unsigned long long ull; using namespace std; #define vll vector> #define vl vector #define vi vector #define vii vector> #define pb push_back #define pf push_front #define ld long double #define Sort(a) sort(a.begin(),a.end()) #define cSort(a,cmp) sort(a.begin(),a.end(),cmp) #define reSort(a) sort(a.rbegin(), a.rend()) static const ll llMAX = numeric_limits::max(); static const int intMAX = numeric_limits::max(); static const ll llMIN = numeric_limits::min(); static const int intMIN = numeric_limits::min(); static const ll d_5 = 100000; static const ll d9_7 = 1000000007; static const ll d_9 = 1000000000; static const double PI=3.14159265358979323846; template void Printvector(std::vector a){ int size = a.size(); rep(i,size){ cout< void Printvector(std::vector> a){ int size = a.size(); rep(i,size){ int size2=a[i].size(); rep(j,size2){ cout< prime_factor(int64_t n) { unordered_map ret; for(int64_t i = 2; i * i <= n; i++) { while(n % i == 0) { ret[i]++; n /= i; } } if(n != 1) ret[n] = 1; return ret; } template vector getaccum(vector a){ int size=a.size(); vector ans(size); ans[0]=a[0]; for(int i=0;i accum,int l,int r){//閉区間[l,r]の総和 if(l==0){ return accum[r]; }else{ return accum[r]-accum[l-1]; } } struct datas{ ll path; ll total; };/* bool cmp(const datas &a, const datas &b) { return a.num < b.num; }*/ struct opera{ int x; int y; datas state; }; int LCS(string s,string t) { int n=s.size(); int m=t.size(); vector> dp=vector>(n+1,vector(m+1,0)); for (int i=0; i class delaySegmentTree{ public: int n; vector nodes; vector delaynodes; //constructor delaySegmentTree(int size,T init){ initialize(size,init); } //関数を定義(minなのかmaxなのかとか) T thisoperator(T a, T b){ return max(a,b); } void update(int x,T a){//xがindex x+=n-1; nodes[x]=a; while(x>0){ x=(x-1)/2; nodes[x]=thisoperator(nodes[2*x+1],nodes[2*x+2]); } } void initialize(int inputn,T init){ int k=0; while(inputn>(1< (n*2-1,init); delaynodes=vector (n*2-1,-1); } // k 番目のノードについて遅延評価を行う void eval(int k, int l, int r) { //kが今見てるindex // 遅延配列が空でない場合、自ノード及び子ノードへの // 値の伝播が起こる if(delaynodes[k] != -1) { nodes[k] = max(delaynodes[k],nodes[k]); // 最下段かどうかのチェックをしよう // 子ノードは親ノードの 1/2 の範囲であるため、 // 伝播させるときは半分にする if(r - l > 1) { delaynodes[2*k+1] = max(delaynodes[k],delaynodes[2*k+1]) ; delaynodes[2*k+2] = max(delaynodes[k],delaynodes[2*k+2]) ; } // 伝播が終わったので、自ノードの遅延配列を空にする delaynodes[k] = -1; } } void maximaze(int a, int b, T x, int k=0, int l=0, int r=-1) { if(r < 0) r = n; // k 番目のノードに対して遅延評価を行う eval(k, l, r); // 範囲外なら何もしない if(b <= l || r <= a) return; // 完全に被覆しているならば、遅延配列に値を入れた後に評価 if(a <= l && r <= b) { delaynodes[k] = max(x,delaynodes[k]); eval(k, l, r); } // そうでないならば、子ノードの値を再帰的に計算して、 // 計算済みの値をもらってくる else { maximaze(a, b, x, 2*k+1, l, (l+r)/2); maximaze(a, b, x, 2*k+2, (l+r)/2, r); nodes[k] = max(nodes[2*k+1] , nodes[2*k+2]); } } T sec_get(int reql,int reqr,int nowindex=0,int nowl=0,int nowr=-1){ // 最初に呼び出されたときの対象区間は [0, n) //回区間に注意 if(nowr < 0) nowr = n; // 要求区間と対象区間が交わらない -> 適当に返す if(nowr <= reql || reqr <= nowl) return 0; //関数が呼び出されるので評価 eval(nowindex,nowl,nowr); // 要求区間が対象区間を完全に被覆 -> 対象区間を答えの計算に使う if(reql <= nowl && nowr <= reqr) return nodes[nowindex]; // 要求区間が対象区間の一部を被覆 -> 子について探索を行う // 左側の子を vl ・ 右側の子を vr としている // 新しい対象区間は、現在の対象区間を半分に割ったもの T val1 = sec_get(reql, reqr, 2*nowindex+1, nowl, (nowl+nowr)/2); T val2 = sec_get(reql, reqr, 2*nowindex+2, (nowl+nowr)/2, nowr); return thisoperator(val1, val2); } void Printn(){ cout<<"Printvector"<>n; vl a(n); rep(i,n){ cin>>a[i]; } map,greater> nums; rep(i,n){ if(nums.find(a[i])==nums.end()){ nums.insert({a[i],{(int)i}}); }else{ nums[a[i]].insert(i); } } delaySegmentTree seg(n,(ll)0); rep(i,n){ seg.update(i,a[i]); } /*for(auto itr=nums.begin();itr!=nums.end();){ if(itr->second.size()==1){ nums.erase(itr++); }else{ itr++; } }*/ for(auto i:nums){ int l=*i.second.begin(); int r=*--i.second.end(); //cout<