C. Multiplicity(DP + 分解因数)

Problem - C - Codeforces

给定一个整数数组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 
#include 
#include 
#include
#include
#include
#include