CSU 1810 Reverse(贡献)

$给你一个n长度的数字串，问你对所有区间翻转之后得到的数字的和$
$这种数不清的区间的题一般都是贡献$
$计算每一位对答案的贡献$
$一共有\frac{n(n+1)}{2}个区间，然后有包含i的区间有i(n-i+1)个$
$然后不包含i的区间，贡献就是i*(\frac{n(n+1)}{2}-i(n-i+1))*10^i$
$然后考虑包含的，首先就是考虑[i,n-i+1]这段区间，i往这里面每个位置交换$
$都有min(i,n-i+1)个区间$
$所以可以提前预处理出10的次方的和sum[i]$
$贡献就是min(i,n-i+1)*(sum[i]-sum[pos-1])*i$
$然后考虑这个区间外面的两侧，那么就是随着位置的变化，区间的个数在减少$
$但是可以找到规律$
$所以也可以预处理10的次方的前缀和后缀$
$套路比较明显，慢慢推就行了$

---

代码：

#include <map>
#include <set>
#include <stack>
#include <queue>
#include <cmath>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <sstream>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#pragma comment(linker,"/STACK:102400000,102400000")

using namespace std;
#define   MAX           100005
#define   MAXN          1000005
#define   maxnode       205
#define   sigma_size    2
#define   lson          l,m,rt<<1
#define   rson          m+1,r,rt<<1|1
#define   lrt           rt<<1
#define   rrt           rt<<1|1
#define   middle        int m=(r+l)>>1
#define   LL            long long
#define   ull           unsigned long long
#define   mem(x,v)      memset(x,v,sizeof(x))
#define   lowbit(x)     (x&-x)
#define   pii           pair<int,int>
#define   bits(a)       __builtin_popcount(a)
#define   mk            make_pair
#define   limit         10000

//const int    prime = 999983;
const int    INF   = 0x3f3f3f3f;
const LL     INFF  = 0x3f3f;
//const double pi    = acos(-1.0);
const double inf   = 1e18;
const double eps   = 1e-9;
const LL     mod   = 1e9+7;
const ull    mx    = 133333331;

/*****************************************************/
inline void RI(int &x) {
      char c;
      while((c=getchar())<'0' || c>'9');
      x=c-'0';
      while((c=getchar())>='0' && c<='9') x=(x<<3)+(x<<1)+c-'0';
 }
/*****************************************************/

char s[MAX];
LL fact[MAX];
LL pre[MAX];
LL suf[MAX];
LL sum[MAX];

int main() {
    int n;
    while(cin>>n){
        scanf("%s",s+1);
        suf[0]=0;
        fact[n]=1;
        sum[n+1]=0;
        pre[n+1]=0;
        fact[n+1]=1;
        sum[n]=fact[n];
        pre[n]=1;
        for(int i=n-1;i>0;i--){
            fact[i]=fact[i+1]*10%mod;
            sum[i]=(sum[i+1]+fact[i])%mod;
            pre[i]=(pre[i+1]+fact[i]*(n-i+1))%mod;
        }
        for(int i=1;i<=n;i++) suf[i]=(fact[i]*i%mod+suf[i-1])%mod;
        LL tmp=(LL)n*(n+1)/2;
        LL ans=0;
        for(int i=1;i<=n;i++){
            LL l=i;
            LL r=(n-i+1);
            ans+=((tmp-l*r)%mod)*(s[i]-'0')*fact[i]%mod;
            if(ans>=mod) ans-=mod;
            //cout<<ans<<endl;
            if(l<r){
                int pos=n+1-i;
                ans+=l*((sum[i]-sum[pos+1]+mod)%mod)*(s[i]-'0')%mod;
                //cout<<ans<<endl;
                if(ans>=mod) ans-=mod;
                if(pos<n){
                    ans+=pre[pos+1]*(s[i]-'0')%mod;
                    if(ans>=mod) ans-=mod;
                }
                if(i-1){
                    ans+=suf[i-1]*(s[i]-'0')%mod;
                    if(ans>=mod) ans-=mod;
                }
                //cout<<ans<<endl;
            }
            else{
                int pos=n+1-i;
                ans+=r*((sum[pos]-sum[i+1]+mod)%mod)*(s[i]-'0')%mod;
                //cout<<ans<<endl;
                if(ans>=mod) ans-=mod;
                if(pos-1){
                    ans+=suf[pos-1]*(s[i]-'0')%mod;
                    if(ans>=mod) ans-=mod;
                }
                if(i<n){
                    ans+=pre[i+1]*(s[i]-'0')%mod;
                    if(ans>=mod) ans-=mod;
                }
                //cout<<ans<<endl;
            }
        }
        cout<<(ans%mod+mod)%mod<<endl;
    }
    return 0;
}