ps:可能组合数一不小心打错了,请发现的大佬提出,谢谢。
我们来讨论每一位数$a_{i}$被算了多少次。
总共有$n-1$个空位可以放$'+'$所以,$a_{i}$左边有$i-1$个空位,右边$n-1-(i-1)$个。
举个例子来说~~(手动模拟一下)~~,如果数$a_{i}$右边有一个加号,则剩下$n-2$个空位放$k-1$个加号的方案里,每种方案$a_{i}$都当做各位记录到答案里,所以对答案影响:$C^{k-1}_{n-2}*$ $a_{i}$。
假设右边第一个不放,第二个放加号,则还有$n-3$个空位,$k-2$个加号。(因为$a_{i}$右边的以为不能放加号)对答案影响: $C_{n-3}^{k-2}*10*$ $a_{i}$。
$emm……$规律好像找出来了。
这样一直递推下去~~(找规律)~~的话。到$a_{n-1}$ 和 $a_{n}$ 位填的话。对答案的影响: $C^{1}_{n-k-2}*10^{k-1}*a_{i}$
但是,我们发现,最后一位数的后面好像并不能放加号。也就是说,最后一位需要特判。。。
所以,当右边全空着的时候对答案的贡献 $10^{i}*C_{n-i-1}^{i}*a_{i}$
全部统计起来,输出就好。
代码奉上。
#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
const long long p=1e9+7;
int n,k;
long long ans[100010];
char a[100010];
long long pen[100010];//阶乘
long long rpen[100010];//阶乘逆元
long long num[100010];//系数
long long pow2(long long a,long long b)
{
long long res=1;
for(;b;b>>=1,a=a*a%p) if(b&1) res=res*a%p;
return res%p;
}
void rebiut()
{
num[0]=pen[0]=rpen[0]=1;
for(int i=1;i<=n;i++)
pen[i]=pen[i-1]*i%p,
rpen[i]=pow2(pen[i],p-2),
num[i]=10*num[i-1]%p;
return ;
}
long long C(long long n,long long k)//计算组合数,n下面,k上面
{
if(k<0||k>n||n<0) return 0;
else return ((pen[n]*rpen[k])%p)*rpen[n-k]%p;
}
long long now;
int main()
{
scanf("%d%d",&n,&k);
scanf("%s",a);
rebiut();
for(int i=0;i<n;i++) now=(now+num[i]*(a[n-i-1]-'0'))%p;
if(k==0) {printf("%lld ",now);return 0;}
int d=n-k-1;
ans[0]=C(n-2,k-1);
for(int i=0;i<=d;i++) ans[i]=(ans[i-1]+num[i]*C(n-i-2,k-1))%p;
for(int i=d+1;i<n;i++) ans[i]=ans[i-1];
for(int i=0;i<=d;i++) ans[i]=(ans[i]+num[i]*C(n-i-2,k))%p;
now=0;
for(int i=0;i<n;i++) now=(now+ans[n-i-1]*(a[i]-'0'))%p;
printf("%lld ",now);
return 0;
}