牛客练习赛53 B.美味果冻分块

最新推荐文章于 2021-01-22 22:14:49 发布

原创最新推荐文章于 2021-01-22 22:14:49 发布 · 360 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

数论同时被 2 个专栏收录

51 篇文章

订阅专栏

牛客网

26 篇文章

订阅专栏

本文深入探讨了一种将中文题意转化为数学表达式的算法优化方法，通过枚举和统计技巧，将复杂的问题简化为易于计算的形式。文章详细解释了为何选择反向枚举，并介绍了如何利用数组模拟快速幂运算，实现O(nlogn)的时间复杂度。

传送门

题意：中文题意很好理解

思路：

我们将转换成

如何证明这步转化正确呢？
对于i=1，有j=1
对于i=2，有j=1,2
。。。
对于i=n，有j=1,2,...,n
因此如果我们反过来枚举j的话可以发现，
某个i出现的前提就是i≥j，①得证。

还有一个问题，为什么要反过来枚举j呢？
因为原来的角度入手的话， $\small \left \lfloor\frac{i}{j} \right \rfloor$ 是考虑 $\small j$ 是 $\small i$ 的因子，
但变成枚举j之后， $\small \left \lfloor\frac{i}{j} \right \rfloor$ 就变成考虑i是j的倍数。在计算上会很方便。
同时， $\small \left \lfloor\frac{i}{j} \right \rfloor ^{j}$ 的指数 j 在 i 变化时是相对不变的，便于统计。
有了①后，我们设 $\small \left \lfloor\frac{i}{j} \right \rfloor=k$ ,可以得到此时的 i∈[ kj , min((k+1)*j−1,n) ]

那对于每一段 $\small 那对于每一段⌊⌊ijji⌋=k⌋=k,我们有$ ,我们有 $\small \left \lfloor\frac{i}{j} \right \rfloor =k$ ，我们有 (k*j+min((k+1)*j−1,n))∗(min((k+1)*j−1,n)−k*j+1)/2∗k^j （统计kj的总和）

因此我们可以得到

注意到上述式子中的 $\small k^j$ ，它的取值其实是 $\small 1^1,2^1,3^1...,1^2,2^2,3^2,...$
所以其实这个部分可以利用数组去模拟(n≤3e6)，一开始数组里每个值都是1，每次更新n/j个（因为只用到n/j个），让a[i]∗=i，这样就可以等效替代快速幂。
总复杂度大约为O(nlogn)。

原博主

#include<cstdio>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<ctime>
#include<cmath>
#define ll long long
using namespace std;
const ll mod=1e9+7;
const int maxn=3e6+10;
ll power(ll a,ll b)
{
    ll ans=1;
    while(b)
    {
        if(b&1)
            ans=ans*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return ans;
}
ll per[maxn];
void init(ll limit)
{
    for(int i=1;i<=limit;i++)
        per[i]=per[i]*i%mod;
}
int main()
{
    ll inv2=power(2,mod-2);
    ll n;
    cin>>n;
    for(int i=1;i<=n;i++)
    {
        per[i]=1;
    }
    ll sum=0;
    for(int j=1;j<=n;j++)
    {
        init(n/j);
        for(int i=1;i<=n/j;i++)
        {
            ll l=min(n,1ll*(i+1)*j-1);
            sum=(sum+(l-i*j+1)*(i*j+l)%mod*inv2%mod*per[i]%mod)%mod;
        }
    }
    cout<<sum<<endl;
    return 0;
}