C. Multiplicity 简单数论+dp(dp[i][j]=dp[i-1][j-1]+dp[i-1][j] 前面序列要满足才能构成后面序列)+sort

题意:给出n 个数 的序列 问 从n个数删去任意个数  删去的数后的序列b1 b2 b3 ......bk  k|bk

思路: 这种题目都有一个特性 就是取到bk 的时候 需要前面有个bk-1的序列前置  这个时候暴力会多一个n 的复杂度

所以只要定义一个状态(j)表示选择了j个数 这个时候就可以转移到j+1 了

定义状态:dp[i][j] 前i个数 选择了j个 

dp[i][j]=dp[i-1][j-1]+dp[i-1][j] ( j|a[i] ) 这个 选+不选

dp[i][j]=dp[i-1][j]    ( j|a[i]不成立 )

这里无法用n^2的复杂度过 而 我们知道 一个数的因子数可以用sqrt(j)的时间求出来 但是j 和a[i]/j 两个因子的大小不确定 所以就会影响dp进程 因为dp要从j到j+1从小到大转移(因为二维开不下 需要滚动 不然可以随便顺序)

( 数的因子是很稀疏的 所以不会超时  )

 

 1 #include<bits/stdc++.h>
 2 #define FOR(i,f_start,f_end) for(int i=f_start;i<=f_end;i++)
 3 #define MS(arr,arr_value) memset(arr,arr_value,sizeof(arr)) 
 4 #define F first 
 5 #define S second
 6 #define pii pair<int ,int >
 7 #define mkp make_pair
 8 #define pb push_back
 9 #define arr(zzz) array<ll,zzz>
10 #define ll long long 
11 using namespace std;
12 const int maxn=1e6+100;
13 const int inf=0x3f3f3f3f;
14 const int mod=1e9+7;
15 int a[maxn];
16 int dp[1000000+5];
17 int main(){
18        int n;
19     scanf("%d",&n);
20     for(int i=0;i<n;i++)scanf("%d",&a[i]);
21     ll ans=0;
22     int p=1;
23     dp[0]=1;
24     for(int i=0;i<n;i++){
25             vector<int>v(sqrt(a[i]));
26         for(int j=1;j*j<=a[i];j++){
27             if(a[i]%j==0){
28                 v.pb(j);
29                 if(a[i]/j!=j)v.pb(a[i]/j);
30             }
31         }
32         sort(v.begin(),v.end(),[](int a,int b){return a>b;});
33         for(auto p:v){
34             dp[p]=(1ll*dp[p-1]+dp[p])%mod;
35         }
36     }
37     for(int i=1;i<=1000000;i++)ans+=dp[i],ans%=mod;
38     cout<<ans<<endl;
39     return 0;
40 }
View Code

 

全部评论

相关推荐

来,说点可能被同行“骂”的大实话。🙊当初接数字马力Offer时,朋友都说:“蚂蚁的“内包”公司?你想清楚啊!”但入职快一年后的今天,我反而对他有了不一样的看法!🔹&nbsp;是偏见?还是信息差!之前没入职之前外面都在说什么岗位低人一等这类。实际上:这种情况不可至否,不能保证每个团队都是其乐融融。但我在的部门以及我了解的周边同事都还是十分好相处的~和蚂蚁师兄师姐之间也经常开一些小玩笑。总之:身份是蚂蚁公司给的,地位是自己挣的(一个傲娇女孩的自述)。🔹&nbsp;待遇?玩的就是真实!试用期工资全额发!六点下班跑得快(早9晚6或者早10晚7,动态打卡),公积金顶格交。别听那些画饼的,到手的钱和下班的时间才是真的(都是牛马何必难为牛马)。🔹&nbsp;能不能学到技术?来了就“后悔”!我们拥有权限直通蚂蚁知识库,技术栈多到学不完。说“学不到东西”的人,来了可能后悔——后悔来晚了(哈哈哈哈,可以不学但是不能没有)!💥&nbsp;内推地址:https://app.mokahr.com/su/ueoyhg❗我的内推码:NTA6Nvs走我的内推,可以直达业务部门,面试流程更快速,进度可查!今天新放HC,之前挂过也能再战!秋招已经正式开始啦~机会就摆在这,敢不敢来试一试呢?(和我一样,做个勇敢的女孩)
下午吃泡馍:数字马力的薪资一般哇,5年经验的java/测试就给人一万出头,而且刚入职第三天就让人出差,而且是出半年
帮你内推|数字马力 校招
点赞 评论 收藏
分享
评论
点赞
收藏
分享

创作者周榜

更多
牛客网
牛客网在线编程
牛客网题解
牛客企业服务