CF1174E Expected GCD Problem 数论,DP

https://codeforces.com/contest/1174/problem/E

题意:定义gi是排列p1,p2...pi的的GCD(长度为i的前缀GCD),f(p)是 g1,g2..gn中独特的元素个数

fmax(n) 成为f(p)在所有整数1,2...n的排列中的最大值,给出整数n,给出满足f(p)=fmax(n)的排列个数 mod(1e9+7)

思路:对于所有满足f(p)=fmax(n)的排列,它们的首个元素s必然有最多的质因数。每次gcd改变时,只能从其中拿走一个质因数,这样我们可以保证有尽可能多的独特gcd。 而对于s有两条推论:

推论1:s=2^x*3^y,即s只能被2和3整除,因为如果s有其他的质因数p(p>4),我们可以s/p*4,这样可以得到更多质因数。

推论2:y<=1, 因为如果s=2^x*3^y,y>=2,我们可以让s=s/9*8,变为2^(x+3)*3^(y-2),由此得到更多质因数

建立dp[i][x][y], 表示到下标i为止,满足要求的排列个数,并且gcd为2^x*3^y. 定义一个函数f(x,y)表示2^x*3^y的倍数个数,且这些倍数<=n

对于p(i+1)有3个方程:

加上2^x*3^y的倍数,gcd不改变,可以加f(x,y)个数,但是已经加了i个

dp[i+1][x][y]=dp[i+1][x][y]+dp[i][x][y](f(x,y)i)

x减少1,加上是2^(x-1)*3^y的倍数同时不是2^x*3^y的倍数

dp[i+1][x1][y]=dp[i+1][x1][y]+dp[i][x][y](f(x1,y)f(x,y))

y减少1,加上是2^x*3^(y-1)的倍数同时不是2^x*3^y的倍数

dp[i+1][x][y1]=dp[i+1][x][y1]+dp[i][x][y](f(x,y1)f(x,y))

总是可以以2^x开始,所以dp[1][x][0]=1;

同时如果2^(x-1)*3<=n,也可以以此开始,所以dp[1][x-1][1]=1

ans=dp[n][0][0]

#include <iostream>
using namespace std;
#define mod 1000000007
int n,dp[1000005][21][2];
int f(int x,int y)
{
    int tmp=(1<<x);
    if (y)
    tmp*=3; return n/tmp; } int main() { scanf("%d",&n); int p=0; while ((1<<p)<=n) p++; p--; dp[1][p][0]=1; if ((1<<(p-1))*3<=n) dp[1][p-1][1]=1; for (int i=1;i<n;i++) { for (int x=0;x<=p;x++) { for (int y=0;y<=1;y++) { dp[i+1][x][y]=(dp[i+1][x][y]+1LL*dp[i][x][y]*(f(x,y)-i))%mod; if (x) dp[i+1][x-1][y]=(dp[i+1][x-1][y]+1LL*dp[i][x][y]*(f(x-1,y)-f(x,y)))%mod; if (y) dp[i+1][x][y-1]=(dp[i+1][x][y-1]+1LL*dp[i][x][y]*(f(x,y-1)-f(x,y)))%mod; } } } printf("%d",dp[n][0][0]); }

 

转载于:https://www.cnblogs.com/hanker99/p/10978931.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
IndentationError: expected an indented block是Python语言中常见的错误,提示需要缩进一个代码块。这种错误通常是由于缩进不正确导致的,可能是使用了混合tab和空格的缩进,或者没有正确缩进代码。在编译时会报错,并提示出现的行数。解决方法是在出错的那一行使用空格或Tab键进行缩进,但不能混用。比如,对于if语句、函数定义、for循环等,都需要在冒号后面进行缩进。 关于Flake8(E999)错误,这是指Flake8代码检查工具发现的错误,E999是其错误代码。具体的错误类型可能有多种,需要查看具体的错误提示来确定。Flake8是一个常用的Python代码检查工具,可以帮助开发者发现代码中的潜在问题并提供修复建议。如果出现Flake8(E999)错误,可以根据错误提示来解决问题,遵循Python的编码规范和最佳实践来修改代码。 综上所述,IndentationError: expected an indented block是Python代码缩进错误的提示,而Flake8(E999)是Flake8代码检查工具发现的错误。需要根据具体的错误提示来进行相应的修改和调整。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* *3* [IndentationError:expected an indented block错误解决](https://blog.csdn.net/qq_31457413/article/details/98947128)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] - *2* [vs没报错leetcode报错-C-Practice:C-实践](https://download.csdn.net/download/weixin_38517997/19923434)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] [ .reference_list ]

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值