给定一个整数数组a1,a2,...,an。
如果可以从a中删除一些元素得到b,则称数组b为a的子序列。
当且仅当对于每个i(1≤i≤k),bi是i的倍数时,数组b1,b2,...,bk被称为好。
在模109+7下找到a中好的子序列的数量。
如果两个子序列的包含数字的索引集合不同,则认为它们是不同的。也就是说,在比较子序列时,元素的值不重要。特别地,数组a恰好有2n-1个不同的子序列(不包括空子序列)。
输入
第一行包含一个整数n(1≤n≤100000)- 数组a的长度。
接下来一行包含整数a1,a2,...,an(1≤ai≤106)。
输出
仅打印一个整数-在模109+7下取的好子序列数量。
Examples
input
Copy
2 1 2
output
Copy
3
input
Copy
5 2 2 1 22 14
output
Copy
13
在第一个例子中,所有三个非空可能的子序列都是好的:{1},{1,2},{2}。在第二个例子中,可能的好子序列为:{2},{2,2},{2,22},{2,14},{2},{2,22},{2,14},{1},{1,22},{1,14},{22},{22,14},{14}。
请注意,某些子序列会列在多次,因为它们在原始数组中出现多次。
题解:
首先分解因数,把每个数的所有因数,存到其vector数组v[i]中,接着dp即可,
dp[i][j]:前i个且第i个的长度为j的序列个数
dp[i][j]+=dp[i-1][j]; 前i个,且当前i为终点长度为j的序列
if(a[i]%j==0) dp[i][j]+=dp[i-1][j-1];
肯定还要接着进行优化,由于j并不是都满足,所以我们利用前面我们分解的因数可以优化(dp[j]代表长度为j的b数组有多少种答案)
dp[v[i][j]] = dp[v[i][j]] + dp[v[i][j] - 1]肯定由比当前因子少一的转移过来
为啥10000*sqrt(1e6)不会t,因为题中保证了,数组a恰好有2n-1个不同的子序列(不包括空子序列)
还要就是dp过程,类似01背包只能拿一次,所以从后往前DP
#include <cstdio>
#include <cstring>
#include <algorithm>
#include<iostream>
#include<vector>
#include<set>
#include<map>
#include<cmath>
#include<queue>
using namespace std;
typedef long long ll;
#define int long long
typedef pair<int,int> PII;
int mod = 1e9 + 7;
vector<int> v[100050];
int dp[1000050];
void solve()
{
int n;
cin >> n;
for(int i = 1;i <= n;i++)
{
int x;
cin >> x;
for(int j = 1;j*j <= x;j++)
{
if(x%j == 0)
{
if(j*j == x)
{
v[i].push_back(j);
}
else
{
v[i].push_back(j);
v[i].push_back(x/j);
}
}
}
sort(v[i].begin(),v[i].end());
}
dp[0] = 1;
for(int i = 1;i <= n;i++)
{
for(int j = v[i].size() - 1;j >= 0;j--)
{
dp[v[i][j]] = (dp[v[i][j]] + dp[v[i][j] - 1])%mod;
}
}
int ans = 0;
for(int i = 1;i <= 100000;i++)
{
ans = (ans + dp[i])%mod;
}
cout << ans;
}
signed main()
{
// ios::sync_with_stdio(0 );
// cin.tie(0);cout.tie(0);
int t = 1;
// cin >> t;
while(t--)
{
solve();
}
}