给定一个序列 A,求有多少非空序列 B 满足 B 是 A 的子序列 并且 ∀ k ∈ [1,lenb], k∣bk∀ k ∈ [1,lenb], k∣bk,
1 | // f[i][j]= f[i-1][j] + (f[i-1][j-1] ,a[i]%j==0 ) |
优化: 枚举 a[i] 的因数作为 j
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 | #include <iostream> #include<vector> #include <algorithm> #define IOS std::ios::sync_with_stdio(0) using namespace std; const int N =1e5+3,M=1e6+3, mod=1e9+7; int n,a[N],f[M] ; vector< int > fac; // f[i][j]= f[i-1][j] + (f[i-1][j-1] ,a[i]%j==0 ) void sov(){ int i,j; f[0]=1; for (i=1;i<=n;i++){ fac.clear(); for (j=1;j*j<=a[i];j++){ if (a[i]%j==0){ fac.push_back(j); if (j*j!=a[i]) fac.push_back(a[i]/j); } } sort(fac.begin(),fac.end()); for (j=fac.size()-1;j>=0;j--) f[fac[j]]=f[fac[j]]+f[fac[j]-1],f[fac[j]]%=mod; } int ans=0; for (i=1;i<=n;i++) ans+=f[i],ans%=mod; cout<<ans<<endl; } signed main(){ cin>>n; for ( int i=1;i<=n;i++) cin>>a[i]; sov(); } |
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 阿里最新开源QwQ-32B,效果媲美deepseek-r1满血版,部署成本又又又降低了!
· 单线程的Redis速度为什么快?
· SQL Server 2025 AI相关能力初探
· AI编程工具终极对决:字节Trae VS Cursor,谁才是开发者新宠?
· 展开说说关于C#中ORM框架的用法!