初中OJ1599. 【GDKOI2004】香樟树

44 篇文章 0 订阅
32 篇文章 0 订阅

这道题有些难度,分享一下解题方法。
题目:https://jzoj.net/junior/#main/show/1599

60分:

简单的DP。

设F[i]表示到一直第i棵树时选出的树数量的最大值。
F[i]=max(F[j])   (gcd(a[i],a[j])>1)

100分:

用DP会超时。

两个数不互质,就说明这两个数有相同的质因子。
而每棵树的叶子个数<=100000,而2*3*5*7*11*13=30030,所以一个数不同的质因子最多只有6个。

设g[i]表示结尾包含质因数i的序列的最长长度。
第i个数能构成最长的序列,就是其含有的质因子能构成最长序列加上第i个数本身。

则到第i个数能构成的最长序列的长度就是:
max(g[p[j]])+1 (P为其质因子)

算出长度以后,用这个长度去更新所有的g[p[j]]。

分解质因数:

其实只需要把1~sqrt(n)中的质数去除就可以了。
因为有可能会剩下一个很大的质数,所以如果剩下的数>1,就把它加入质因子中。

提示:

sqrt(100000)≈317
1~317中共有66个质数

核心代码:

                while (j<m) and (k<>1) do //分解质因数
                begin
                        inc(j);
                        if k mod ss[j]=0 then
                        begin
                                inc(l);
                                p[l]:=ss[j];
                                while k mod ss[j]=0 do
                                k:=k div ss[j];
                        end;
                end;
                if k>1 then //剩余的大质数
                begin
                        inc(l);
                        p[l]:=k;
                end;
                s:=0;
                for j:=1 to l do
                s:=max(s,g[p[j]]); //寻找最大长度
                inc(s);
                for j:=1 to l do
                g[p[j]]:=s; //更新长度
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值