CF1061C Multiplicity 题解

CF1061C Multiplicity

题意：

从序列$ {a_1, a_2, … , a_n}$ 中选出非空子序列 $b_1,b_2,..,b_k$，一个子序列合法需要满足 $\forall i\in [1,k]$, $i\mid b_i$。求有多少互不相等的合法子序列，答案对 $10^9+7$ 取模。

做法

：先考虑暴力，设 $f_{i,j}$ 为枚举到 $i$ ，且该数作为 $b_j$时的方案数。这样转移是 $O(n^2)$ 的。观察到只会在因数处产生贡献，所以直接枚举 $a_i$ 的因数，在因数处转移。滚动数组优化。

Code:

#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
const int N=1e5+5;
const int M=1e6+5;
int n;
int cnt[M];
int a[M],b[M],tot;
bool cmp(int x,int y){
	return x>y;
}
int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++){
		scanf("%d",&a[i]);
	}
	cnt[0]=1;
	for(int i=1;i<=n;i++){
		int t=sqrt(a[i]);
		tot=0;
		for(int j=1;j<=t;j++){
			if(a[i]%j==0){
				b[++tot]=j;
				if(j*j!=a[i]) b[++tot]=a[i]/j;
			}
		}
		sort(b+1,b+1+tot,cmp);
		for(int j=1;j<=tot;j++){
			cnt[b[j]]=(cnt[b[j]-1]+cnt[b[j]])%mod;
		}
	}	
	int ans=0;
	for(int i=1;i<=n;i++){
		ans=(ans+cnt[i])%mod;
	}
	printf("%d\n",ans);
	return 0;
}

启示：

1. 对于方案数较为稀疏的，我们可以设"到k且选k"而非"前k"来表示状态。
2. 注意复杂度的计算。$n$ 以内的因数个数有小于根号n个。