CSU1810(2016年湖南省第十二届大学生计算机程序设计竞赛 - H)

最新推荐文章于 2021-07-17 16:11:40 发布

FlushHip

最新推荐文章于 2021-07-17 16:11:40 发布

阅读量1.2k

点赞数 1

分类专栏： OJ ACM算法竞赛

本文链接：https://blog.csdn.net/FlushHip/article/details/52429704

版权

OJ 同时被 2 个专栏收录

138 篇文章 2 订阅

订阅专栏

ACM算法竞赛

130 篇文章 37 订阅

订阅专栏

Problem : Reverse
Description :

给你一个大整数。按照图中给的公式反转字符串然后加起来。输出最后答案。

Solution :

列举找规律+乘法逆元。举个列子吧。 $12345$ ，我们按照公式展开。得到
$12345 - 1111 - 111 - 11 - 1 21234 - 2345 - 222 - 22 - 2 33123 - 3234 - 345 - 33 - 3 44412 - 4423 - 434 - 45 - 4 55551 - 5552 - 553 - 54 - 5$ $\begin{matrix} 1 & 2 & 3 & 4 & 5 \\ 2 & 1 & 3 & 4 & 5 \\3 & 2 & 1 & 4 & 5 \\4 & 3 & 2 & 1 & 5 \\5 & 4 & 3 & 2 & 1\\ - & - & - & - & -\\1 & 2 & 3 & 4 & 5 \\1 & 3 & 2 & 4 & 5 \\1 & 4 & 3 & 2 & 5 \\1 & 5 & 4 & 3 & 2 \\ - & - & - & - & -\\1 & 2 & 3 & 4 & 5 \\1 & 2 & 4 & 3 & 5 \\1 & 2 & 5 & 4 & 3 \\- & - & - & - & -\\1 & 2 & 3 & 4 & 5 \\1 & 2 & 3 & 5 & 4 \\- & - & - & - & -\\1 & 2 & 3 & 4 & 5 \\\end{matrix}\quad$
我们观察3这个位置，即 $k=3, num=3$ 。在 $i\geq4$ 时， $num$ 一直待在自己的位置上。而在 $i=1$ 有两个待在自己的位置， $i=2$ 有一个待在自己的位置。同时发现 $i=1$ 时 $num$ 从 $1-3$ 都待过， $i=2$ 时 $num$ 从 $2-4$ 都待过， $i=3$ 时 $num$ 从 $3-5$ 都待过。这样，就可以写出一个公式了。 $⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ f (k) = {[( n - k + 1 ) * ( n - k ) 2 + ( k - 1 ) * k 2] * 10 n - k + a * ( 1 - 1 10 k ) 1 - 1 10} * n u m a = 10 k - 1 * ( 10 n - k + 1 - 1 ) 9$ $\left\{ \begin{aligned}f(k) = \{[ \frac{(n-k+1)*(n-k)}{2}+\frac{(k-1)*k}{2}]*10^{n-k}+\frac{a*(1-\frac {1}{10^k})}{1-\frac {1}{10}}\}*num \\ a=\frac {10^{k-1}*(10^{n-k+1}-1)}{9} \end{aligned} \right.$ 把这个公式化简下。得出了每一位的公式，这样，我们枚举每一位，然后把结果加起来就可以了，但是现在有个问题，结果很大，超过了 $long long$ 的范围，因此要对 $(10^9+7)$ 取模。而公式中又用到了除法，我们知道，如果分子取模了再除分母，答案是会不对的。那么这个时候就要用到乘法逆元了。我用的是费马小定理求乘法逆元。因此这个题目还是考察了多个方面。

Code(C++) :

#include <stdio.h>

typedef long long LL;

const int M=100000+5;
const LL MOD=1000000000LL+7;

char num[M];
LL n;

LL mul(LL x,LL n,LL mod)
{
    LL tmp=x;
    LL ans=1;
    while(n){
        if(n&1){
            ans*=tmp;
            ans%=mod;
        }
        tmp*=tmp;
        tmp%=mod;
        n>>=1;
    }
    return ans;
}

LL cal(LL digit,LL k)
{
    LL tmp=(n-k+1)*(n-k)/2+k*(k-1)/2;
    tmp%=MOD;
    tmp=tmp*mul(10,n-k,MOD)%MOD;
    LL a=(mul(10,k-1,MOD)*(mul(10,n-k+1,MOD)-1+MOD))%MOD*mul(9,MOD-2,MOD)%MOD;
    LL t=(a*(mul(10,k,MOD)-1+MOD)%MOD)*(mul(9*mul(10,k-1,MOD)%MOD,MOD-2,MOD))%MOD;
    tmp=(tmp+t)%MOD;
    tmp=tmp*digit%MOD;
    return tmp;
}

int main()
{
    while(~scanf("%lld",&n)){
        scanf("%s",num);
        LL ans=0;
        for(LL i=0;i<n;i++)
            ans=(ans+cal((LL)(num[i]-'0'),i+1))%MOD;
        printf("%lld\n",ans);
    }
    return 0;
}