AtCoder Beginner Contest 132 F

F - Small Products

怎么说呢？目前超了自己水平两个档次的题。
第一个档次熟练掌握$dp$
第二个档次熟练掌握分块

这道题：参考大佬

要求求相邻元素乘积$\le N$的序列数。$N\le 1e9$

序列数考虑$dp$
显然：$dp[i][j]={\sum }_{x\ast j\le N}dp[i-1][x]$
但是显然第二维大小不符合要求。

这种$1e9$一般都是分块或者数学(数学一般是$1e18$)
这道题是分块的做法：值得学习。

我们考虑前$\sqrt{N}$个可以直接利用，后面的部分可以划分成$\sqrt{N}$块。
第$i$块代表了使得$i\ast k\le N$成立且$i$唯一的所有$k$。
意思就是说这里的$k$只能和$i$组合，与$i+1$不能组合(在后面的式子可以发现被去除了)，与$i-1$组合的部分在前面已经有了。

第$i$块的数量为$\frac{N}{i}-\frac{N}{i+1}$
数学证明并不会，但举例子比较容易懂。
比如说第一块，只能被$1$唯一整除的。
$N/1=N$是不考虑唯一的情况
$N/2$是不考虑唯一，被$2$整除的。
唯一被$1$整除的显然是$N-N/2$。以此类推。

接下来是正戏。
我们考虑分类讨论：
第一部分是第$i$位放的是$\sqrt{N}$以内的。$S[i][j]$
第二部分是第$i$位放的是以外的，这个时候要利用分块处理。$B[i][j]$

答案显然是${\sum }_{i=1}^{\sqrt{N}}S[k][i]+B[k][i]$
$S[i][j]$怎么用$dp$解决呢，考虑第$j$位放的是小的数，那么第$j-1$位可以放小的数和大的数。
小的数无所谓，对于大的数只能从$[j,\sqrt{N}]$的块用，小于的会导致数过大，超过$N$.
所以$S[i][j]={\sum }_{k=1}^{\sqrt{N}}S[i-1][k]+{\sum }_{k=j}^{\sqrt{N}}B[i-1][k]$
对于$B[i][j]$，前一位只能放小数，不过这里的$j$表示的是块，所以还需要乘上块数里的数量。
所以$B[i][j]=(\frac{N}{i}-\frac{N}{i+1})\ast {\sum }_{k=1}^{j}S[i-1][k]$超过$j$的显然会超。

这样子做$dp$复杂度是$k\ast \sqrt{N}\ast \sqrt{N}$
考虑加的操作，可以简化为前缀和，利用前缀和可以实现$k\ast \sqrt{N}$
有个小细节就是第$\sqrt{N}$块和放$\sqrt{N}$等价。取其一即可。

#include<bits/stdc++.h>
#define FOR(i,l,r) for(int i=l;i<=r;i++)
#define ll long long
using namespace std;

const ll mod = 1e9+7;

ll n;int k;
ll S[110][202000],B[110][202000];

int main(){
    cin>>n>>k;
    ll d=sqrt(n);
    ll ans=0;
    S[0][1]=1;
    FOR(i,1,k){
        FOR(j,1,d)S[i-1][j]=(S[i-1][j]+S[i-1][j-1])%mod;
        FOR(j,1,d)B[i-1][j]=(B[i-1][j]+B[i-1][j-1])%mod;
        FOR(j,1,d){
            S[i][j]=(S[i-1][d]+B[i-1][d]-B[i-1][j-1]+mod)%mod;
            B[i][j]=(n/j-n/(j+1))*S[i-1][j]%mod;
            if(j==n/j)B[i][j]=0;//SжагаСЫ
        }
    }
    FOR(i,1,d)ans=(ans+S[k][i]+B[k][i])%mod;
    cout<<ans<<endl;
}