#define PROBLEM "https://judge.yosupo.jp/problem/number_of_substrings"
#include"../../template/template.hpp"
#include"../../string/suffix-array.hpp"usingnamespaceNyaan;voidNyaan::solve(){ins(S);SuffixArraysa(S);LCPArraylcp(sa);llans=1LL*sz(S)*(sz(S)+1)/2;each(x,lcp.LCP)ans-=x;out(ans);}
#line 1 "verify/verify-yosupo-string/yosupo-number-of-substrings.test.cpp"
#define PROBLEM "https://judge.yosupo.jp/problem/number_of_substrings"
#line 2 "template/template.hpp"
usingnamespacestd;// intrinstic#include<immintrin.h>#include<algorithm>
#include<array>
#include<bitset>
#include<cassert>
#include<cctype>
#include<cfenv>
#include<cfloat>
#include<chrono>
#include<cinttypes>
#include<climits>
#include<cmath>
#include<complex>
#include<cstdarg>
#include<cstddef>
#include<cstdint>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<deque>
#include<fstream>
#include<functional>
#include<initializer_list>
#include<iomanip>
#include<ios>
#include<iostream>
#include<istream>
#include<iterator>
#include<limits>
#include<list>
#include<map>
#include<memory>
#include<new>
#include<numeric>
#include<ostream>
#include<queue>
#include<random>
#include<set>
#include<sstream>
#include<stack>
#include<streambuf>
#include<string>
#include<tuple>
#include<type_traits>
#include<typeinfo>
#include<unordered_map>
#include<unordered_set>
#include<utility>
#include<vector>// utility#line 3 "template/util.hpp"
namespaceNyaan{usingll=longlong;usingi64=longlong;usingu64=unsignedlonglong;usingi128=__int128_t;usingu128=__uint128_t;template<typenameT>usingV=vector<T>;template<typenameT>usingVV=vector<vector<T>>;usingvi=vector<int>;usingvl=vector<longlong>;usingvd=V<double>;usingvs=V<string>;usingvvi=vector<vector<int>>;usingvvl=vector<vector<longlong>>;template<typenameT>usingminpq=priority_queue<T,vector<T>,greater<T>>;template<typenameT,typenameU>structP:pair<T,U>{template<typename...Args>P(Args...args):pair<T,U>(args...){}usingpair<T,U>::first;usingpair<T,U>::second;P&operator+=(constP&r){first+=r.first;second+=r.second;return*this;}P&operator-=(constP&r){first-=r.first;second-=r.second;return*this;}P&operator*=(constP&r){first*=r.first;second*=r.second;return*this;}template<typenameS>P&operator*=(constS&r){first*=r,second*=r;return*this;}Poperator+(constP&r)const{returnP(*this)+=r;}Poperator-(constP&r)const{returnP(*this)-=r;}Poperator*(constP&r)const{returnP(*this)*=r;}template<typenameS>Poperator*(constS&r)const{returnP(*this)*=r;}Poperator-()const{returnP{-first,-second};}};usingpl=P<ll,ll>;usingpi=P<int,int>;usingvp=V<pl>;constexprintinf=1001001001;constexprlonglonginfLL=4004004004004004004LL;template<typenameT>intsz(constT&t){returnt.size();}template<typenameT,typenameU>inlineboolamin(T&x,Uy){return(y<x)?(x=y,true):false;}template<typenameT,typenameU>inlineboolamax(T&x,Uy){return(x<y)?(x=y,true):false;}template<typenameT>inlineTMax(constvector<T>&v){return*max_element(begin(v),end(v));}template<typenameT>inlineTMin(constvector<T>&v){return*min_element(begin(v),end(v));}template<typenameT>inlinelonglongSum(constvector<T>&v){returnaccumulate(begin(v),end(v),0LL);}template<typenameT>intlb(constvector<T>&v,constT&a){returnlower_bound(begin(v),end(v),a)-begin(v);}template<typenameT>intub(constvector<T>&v,constT&a){returnupper_bound(begin(v),end(v),a)-begin(v);}constexprlonglongTEN(intn){longlongret=1,x=10;for(;n;x*=x,n>>=1)ret*=(n&1?x:1);returnret;}template<typenameT,typenameU>pair<T,U>mkp(constT&t,constU&u){returnmake_pair(t,u);}template<typenameT>vector<T>mkrui(constvector<T>&v,boolrev=false){vector<T>ret(v.size()+1);if(rev){for(inti=int(v.size())-1;i>=0;i--)ret[i]=v[i]+ret[i+1];}else{for(inti=0;i<int(v.size());i++)ret[i+1]=ret[i]+v[i];}returnret;};template<typenameT>vector<T>mkuni(constvector<T>&v){vector<T>ret(v);sort(ret.begin(),ret.end());ret.erase(unique(ret.begin(),ret.end()),ret.end());returnret;}template<typenameF>vector<int>mkord(intN,Ff){vector<int>ord(N);iota(begin(ord),end(ord),0);sort(begin(ord),end(ord),f);returnord;}template<typenameT>vector<int>mkinv(vector<T>&v){intmax_val=*max_element(begin(v),end(v));vector<int>inv(max_val+1,-1);for(inti=0;i<(int)v.size();i++)inv[v[i]]=i;returninv;}vector<int>mkiota(intn){vector<int>ret(n);iota(begin(ret),end(ret),0);returnret;}template<typenameT>Tmkrev(constT&v){Tw{v};reverse(begin(w),end(w));returnw;}template<typenameT>boolnxp(T&v){returnnext_permutation(begin(v),end(v));}// 返り値の型は入力の T に依存// i 要素目 : [0, a[i])template<typenameT>vector<vector<T>>product(constvector<T>&a){vector<vector<T>>ret;vector<T>v;autodfs=[&](autorc,inti)->void{if(i==(int)a.size()){ret.push_back(v);return;}for(intj=0;j<a[i];j++)v.push_back(j),rc(rc,i+1),v.pop_back();};dfs(dfs,0);returnret;}// F : void(T&), mod を取る操作// T : 整数型のときはオーバーフローに注意するtemplate<typenameT,typenameF>TPower(Ta,longlongn,constT&I,F&&f){static_assert(std::is_invocable_r_v<void,F&,T&>,"Power callback must be callable as void(T&)");Tres=I;for(;n;std::invoke(f,a=a*a),n>>=1){if(n&1)std::invoke(f,res=res*a);}returnres;}// T : 整数型のときはオーバーフローに注意するtemplate<typenameT>TPower(Ta,longlongn,constT&I=T{1}){autono_op=[](T&)->void{};returnPower(a,n,I,no_op);}template<typenameT>TRev(constT&v){Tres=v;reverse(begin(res),end(res));returnres;}template<typenameT>vector<T>Transpose(constvector<T>&v){usingU=typenameT::value_type;if(v.empty())return{};intH=v.size(),W=v[0].size();vectorres(W,T(H,U{}));for(inti=0;i<H;i++){for(intj=0;j<W;j++){res[j][i]=v[i][j];}}returnres;}template<typenameT>vector<T>Rotate(constvector<T>&v,intclockwise=true){usingU=typenameT::value_type;intH=v.size(),W=v[0].size();vectorres(W,T(H,U{}));for(inti=0;i<H;i++){for(intj=0;j<W;j++){if(clockwise){res[W-1-j][i]=v[i][j];}else{res[j][H-1-i]=v[i][j];}}}returnres;}}// namespace Nyaan#line 58 "template/template.hpp"
// bit operation#line 1 "template/bitop.hpp"
namespaceNyaan{__attribute__((target("popcnt")))inlineintpopcnt(constu64&a){return__builtin_popcountll(a);}inlineintlsb(constu64&a){returna?__builtin_ctzll(a):64;}inlineintctz(constu64&a){returna?__builtin_ctzll(a):64;}inlineintmsb(constu64&a){returna?63-__builtin_clzll(a):-1;}template<typenameT>inlineintgbit(constT&a,inti){return(a>>i)&1;}template<typenameT>inlinevoidsbit(T&a,inti,boolb){if(gbit(a,i)!=b)a^=T(1)<<i;}constexprlonglongPW(intn){return1LL<<n;}constexprlonglongMSK(intn){return(1LL<<n)-1;}}// namespace Nyaan#line 61 "template/template.hpp"
// inout#line 1 "template/inout.hpp"
namespaceNyaan{template<typenameT,typenameU>ostream&operator<<(ostream&os,constpair<T,U>&p){os<<p.first<<" "<<p.second;returnos;}template<typenameT,typenameU>istream&operator>>(istream&is,pair<T,U>&p){is>>p.first>>p.second;returnis;}template<typenameT>ostream&operator<<(ostream&os,constvector<T>&v){ints=(int)v.size();for(inti=0;i<s;i++)os<<(i?" ":"")<<v[i];returnos;}template<typenameT>istream&operator>>(istream&is,vector<T>&v){for(auto&x:v)is>>x;returnis;}istream&operator>>(istream&is,__int128_t&x){stringS;is>>S;x=0;intflag=0;for(auto&c:S){if(c=='-'){flag=true;continue;}x*=10;x+=c-'0';}if(flag)x=-x;returnis;}istream&operator>>(istream&is,__uint128_t&x){stringS;is>>S;x=0;for(auto&c:S){x*=10;x+=c-'0';}returnis;}ostream&operator<<(ostream&os,__int128_tx){if(x==0)returnos<<0;if(x<0)os<<'-',x=-x;stringS;while(x)S.push_back('0'+x%10),x/=10;reverse(begin(S),end(S));returnos<<S;}ostream&operator<<(ostream&os,__uint128_tx){if(x==0)returnos<<0;stringS;while(x)S.push_back('0'+x%10),x/=10;reverse(begin(S),end(S));returnos<<S;}voidin(){}template<typenameT,class...U>voidin(T&t,U&...u){cin>>t;in(u...);}voidout(){cout<<"\n";}template<typenameT,class...U,charsep=' '>voidout(constT&t,constU&...u){cout<<t;if(sizeof...(u))cout<<sep;out(u...);}structIoSetupNya{IoSetupNya(){cin.tie(nullptr);ios::sync_with_stdio(false);cout<<fixed<<setprecision(15);cerr<<fixed<<setprecision(7);}}iosetupnya;}// namespace Nyaan#line 64 "template/template.hpp"
// debug#line 1 "template/debug.hpp"
namespaceDebugImpl{template<typenameU,typename=void>structis_specialize:false_type{};template<typenameU>structis_specialize<U,typenameconditional<false,typenameU::iterator,void>::type>:true_type{};template<typenameU>structis_specialize<U,typenameconditional<false,decltype(U::first),void>::type>:true_type{};template<typenameU>structis_specialize<U,enable_if_t<is_integral<U>::value,void>>:true_type{};voiddump(constchar&t){cerr<<t;}voiddump(conststring&t){cerr<<t;}voiddump(constbool&t){cerr<<(t?"true":"false");}voiddump(__int128_tt){if(t==0)cerr<<0;if(t<0)cerr<<'-',t=-t;stringS;while(t)S.push_back('0'+t%10),t/=10;reverse(begin(S),end(S));cerr<<S;}voiddump(__uint128_tt){if(t==0)cerr<<0;stringS;while(t)S.push_back('0'+t%10),t/=10;reverse(begin(S),end(S));cerr<<S;}template<typenameU,enable_if_t<!is_specialize<U>::value,nullptr_t>=nullptr>voiddump(constU&t){cerr<<t;}template<typenameT>voiddump(constT&t,enable_if_t<is_integral<T>::value>*=nullptr){stringres;if(t==Nyaan::inf)res="inf";ifconstexpr(is_signed<T>::value){if(t==-Nyaan::inf)res="-inf";}ifconstexpr(sizeof(T)==8){if(t==Nyaan::infLL)res="inf";ifconstexpr(is_signed<T>::value){if(t==-Nyaan::infLL)res="-inf";}}if(res.empty())res=to_string(t);cerr<<res;}template<typenameT,typenameU>voiddump(constpair<T,U>&);template<typenameT>voiddump(constpair<T*,int>&);template<typenameT>voiddump(constT&t,enable_if_t<!is_void<typenameT::iterator>::value>*=nullptr){cerr<<"[ ";for(autoit=t.begin();it!=t.end();){dump(*it);cerr<<(++it==t.end()?"":", ");}cerr<<" ]";}template<typenameT,typenameU>voiddump(constpair<T,U>&t){cerr<<"( ";dump(t.first);cerr<<", ";dump(t.second);cerr<<" )";}template<typenameT>voiddump(constpair<T*,int>&t){cerr<<"[ ";for(inti=0;i<t.second;i++){dump(t.first[i]);cerr<<(i==t.second-1?"":", ");}cerr<<" ]";}voidtrace(){cerr<<endl;}template<typenameHead,typename...Tail>voidtrace(Head&&head,Tail&&...tail){cerr<<" ";dump(head);if(sizeof...(tail)!=0)cerr<<",";trace(std::forward<Tail>(tail)...);}}// namespace DebugImpl#ifdef NyaanDebug
#define trc(...) \
do { \
cerr << "## " << #__VA_ARGS__ << " = "; \
DebugImpl::trace(__VA_ARGS__); \
} while (0)
#else
#define trc(...) (void(0))
#endif
#ifdef NyaanLocal
#define trc2(...) \
do { \
cerr << "## " << #__VA_ARGS__ << " = "; \
DebugImpl::trace(__VA_ARGS__); \
} while (0)
#else
#define trc2(...) (void(0))
#endif
#line 67 "template/template.hpp"
// macro#line 1 "template/macro.hpp"
#define each(x, v) for (auto&& x : v)
#define each2(x, y, v) for (auto&& [x, y] : v)
#define all(v) (v).begin(), (v).end()
#define rep(i, N) for (long long i = 0; i < (long long)(N); i++)
#define repr(i, N) for (long long i = (long long)(N)-1; i >= 0; i--)
#define rep1(i, N) for (long long i = 1; i <= (long long)(N); i++)
#define repr1(i, N) for (long long i = (N); (long long)(i) > 0; i--)
#define reg(i, a, b) for (long long i = (a); i < (b); i++)
#define regr(i, a, b) for (long long i = (b)-1; i >= (a); i--)
#define fi first
#define se second
#define ini(...) \
int __VA_ARGS__; \
in(__VA_ARGS__)
#define inl(...) \
long long __VA_ARGS__; \
in(__VA_ARGS__)
#define ins(...) \
string __VA_ARGS__; \
in(__VA_ARGS__)
#define in2(s, t) \
for (int i = 0; i < (int)s.size(); i++) { \
in(s[i], t[i]); \
}
#define in3(s, t, u) \
for (int i = 0; i < (int)s.size(); i++) { \
in(s[i], t[i], u[i]); \
}
#define in4(s, t, u, v) \
for (int i = 0; i < (int)s.size(); i++) { \
in(s[i], t[i], u[i], v[i]); \
}
#define die(...) \
do { \
Nyaan::out(__VA_ARGS__); \
return; \
} while (0)
#line 70 "template/template.hpp"
namespaceNyaan{voidsolve();}intmain(){Nyaan::solve();}#line 2 "string/suffix-array.hpp"
#line 2 "data-structure/sparse-table.hpp"
#line 6 "data-structure/sparse-table.hpp"
usingnamespacestd;template<typenameT>structSparseTable{inlinestaticconstexprTINF=numeric_limits<T>::max()/2;intN;vector<vector<T>>table;Tf(Ta,Tb){returnmin(a,b);}SparseTable(){}SparseTable(constvector<T>&v):N(v.size()){intb=1;while((1<<b)<=N)++b;table.push_back(v);for(inti=1;i<b;i++){table.push_back(vector<T>(N,INF));for(intj=0;j+(1<<i)<=N;j++){table[i][j]=f(table[i-1][j],table[i-1][j+(1<<(i-1))]);}}}// [l, r)Tquery(intl,intr){assert(0<=landl<=randr<=N);if(l==r)returnINF;intb=31-__builtin_clz(r-l);returnf(table[b][l],table[b][r-(1<<b)]);}};/**
* @brief Sparse Table
*/#line 6 "string/suffix-array.hpp"
// remind: SA including empty string// verify https://judge.yosupo.jp/submission/240structSuffixArray{int_size;vector<int>sa;string&s;SuffixArray(string&str):_size(str.size()),s(str){s.push_back(0);sa.resize(s.size());iota(begin(sa),end(sa),0);sort(begin(sa),end(sa),[&](inta,intb){returns[a]==s[b]?a>b:s[a]<s[b];});vector<int>classes(s.size()),c(s.begin(),s.end()),cnt(s.size());for(intlen=1;len<(int)s.size();len<<=1){for(inti=0;i<(int)s.size();i++){if(i>0&&c[sa[i-1]]==c[sa[i]]&&sa[i-1]+len<(int)s.size()&&c[sa[i-1]+len/2]==c[sa[i]+len/2]){classes[sa[i]]=classes[sa[i-1]];}else{classes[sa[i]]=i;}}iota(begin(cnt),end(cnt),0);copy(begin(sa),end(sa),begin(c));for(inti=0;i<(int)s.size();i++){ints1=c[i]-len;if(s1>=0)sa[cnt[classes[s1]]++]=s1;}classes.swap(c);}s.pop_back();}voidoutput(){cout<<"SA\tidx\tstr"<<endl;for(inti=0;i<size();i++){cout<<i<<": \t"<<sa[i]<<" \t";if(sa[i]!=_size)cout<<s.substr(sa[i],_size-sa[i])<<endl;elsecout<<"$"<<endl;}cout<<endl;}intsize()const{return_size+1;}intoperator[](intk)const{returnsa[k];}};structLCPArray{constSuffixArray&SA;vector<int>LCP,rank;LCPArray(constSuffixArray&sa):SA(sa){LCP.resize(SA.size());rank.resize(SA.size());for(inti=0;i<SA.size();i++){rank[SA[i]]=i;}LCP[0]=0;for(inti=0,h=0;i<SA.size()-1;i++){intj=SA[rank[i]-1];h?h--:h;while((i>j?i:j)+h<SA.size()-1&&SA.s[i+h]==SA.s[j+h]&&++h);LCP[rank[i]-1]=h;}}voidoutput(){cout<<"SA\tidx\tLCP\tstr"<<endl;for(inti=0;i<SA.size();i++){cout<<i<<"\t"<<SA[i]<<" \t"<<LCP[i]<<"\t";if(SA[i]==SA.size()-1)cout<<"$";elsecout<<SA.s.substr(SA[i],SA.size()-1-SA[i]);cout<<endl;}}};// verify// https://onlinejudge.u-aizu.ac.jp/status/users/NyaanNyaan/submissions/1/ALDS1_14_D/judge/3874273/C++14// https://atcoder.jp/contests/abc135/submissions/7574225// https://judge.yosupo.jp/submission/241// https://atcoder.jp/contests/abc141/submissions/7577295structStringSearch{string&s;constSuffixArray&sa;constLCPArray&lcp;SparseTable<int>sparse;StringSearch(LCPArray&lcp):s(lcp.SA.s),sa(lcp.SA),lcp(lcp),sparse(lcp.LCP){}pair<int,int>comp(conststring&t,intlen,intsi,intti=0){intsn=(int)s.size(),tn=(int)t.size();si+=len,ti+=len;while(si<sn&&ti<tn){if(s[si]!=t[ti])returnmake_pair(s[si]<t[ti],ti);si++,ti++;}returnmake_pair((si>=sn&&ti<tn),ti);}pair<int,int>find_range(intleft,intmed,intright,intlen){{intng=left-1,ok=med;while(ng+1<ok){intcur=(ng+ok)/2;if(sparse.query(cur,med)>=len)ok=cur;elseng=cur;}left=ok;}{intok=med,ng=right+1;while(ok+1<ng){intcur=(ng+ok)/2;if(sparse.query(med,cur)>=len)ok=cur;elseng=cur;}right=ok;}returnmake_pair(left,right);}public:// Longest Common Prefix between S[i , N) and S[j , N)intArbitraryLCP(inti,intj){if(i==j)return(int)(s.size())-i;returnsparse.query(min(lcp.rank[i],lcp.rank[j]),max(lcp.rank[i],lcp.rank[j]));}intArbitaryLCP(inti,intj){returnArbitraryLCP(i,j);}// String Search O(|T| + log |S|)// return : [l, r], l and r are indices of Suffix Array// if T doesn't exist, return (-1, -1)pair<int,int>find(string&t){intleft=1,right=sa.size()-1,med=left;intleftlen=0,rightlen=0,tlen=t.size();pair<int,int>ret;while(left+1<right){med=(left+right)/2;intcorres_len=max(min(leftlen,sparse.query(left,med)),min(rightlen,sparse.query(med,right)));if(corres_len<max(leftlen,rightlen)){if(leftlen<rightlen)left=med,leftlen=corres_len;elseright=med,rightlen=corres_len;continue;}ret=comp(t,corres_len,sa[med]);if(ret.second==tlen)returnfind_range(left,med,right,tlen);if(ret.first==0)right=med,rightlen=ret.second;elseleft=med,leftlen=ret.second;}if(sa.size()<=3){if(comp(t,0,sa[left]).second==tlen)returnfind_range(left,left,right,tlen);if(comp(t,0,sa[right]).second==tlen)returnfind_range(left,right,right,tlen);returnmake_pair(-1,-1);}med=left+right-med;ret=comp(t,min(leftlen,rightlen),sa[med]);if(ret.second==tlen)returnfind_range(left,med,right,tlen);returnmake_pair(-1,-1);}};// Usage:// SuffixArray sa(S);// LCPArray lcp(sa);// StringSearch search(lcp);#line 5 "verify/verify-yosupo-string/yosupo-number-of-substrings.test.cpp"
usingnamespaceNyaan;voidNyaan::solve(){ins(S);SuffixArraysa(S);LCPArraylcp(sa);llans=1LL*sz(S)*(sz(S)+1)/2;each(x,lcp.LCP)ans-=x;out(ans);}