洛谷-5153 简单的函数

题目描述
HKE有一次发现了一个很有趣的函数。
定义f(2)=1f。对于n≥3,设tt为最小的使得n不能被t整除的正整数,则f(n)=f(t)+1。
举个栗子。比如n=6,此时t=4,f(6)=f(4)+1=f(3)+2=f(2)+3=4。
现在,HKE想知道 f ( 2 ) × f ( 3 ) × ⋯ × f ( n ) f(2)×f(3)×⋯×f(n) f(2)×f(3)××f(n)是多少?答案可能很大,请对 1 0 9 + 7 10^9+7 109+7
+7取模。

输入格式
一行一个正整数n。

输出格式
一行为所求的结果。

输入输出样例
输入 #1
4

输出 #1
6

说明/提示
对于100%的数据, n ≤ 1 0 18 n\leq10^{18} n1018

解释:设 n = 2 p 1 3 p 2 . . . n=2^{p_1}3^{p_2}... n=2p13p2...
如果 p 1 p_1 p1=0,则很明显 f ( n ) = 2 f(n)=2 f(n)=2,即统计奇数个数 n u m 1 num_1 num1
p 1 ≥ 1 p_1\ge1 p11时,t很明显有两种情况,
( 1 ) t = 2 p 1 + 1 (1)t=2^{p_1+1} (1)t=2p1+1
( 2 ) (2) (2)t=其他的素数组合
如果是第一种情况则 f ( n ) = 4 f(n)=4 f(n)=4
因为 t = p 1 + 1 t=^{p_1+1} t=p1+1,所以, [ 1 , t − 1 ] [1,t-1] [1,t1]都被整除,我们取 [ 1 , t − 1 ] [1,t-1] [1,t1]得最小公倍数 g g g,然后统计出来满足
g ∗ x ≤ n g*x\le n gxn得个数 x x x,其中 x x x不能有2因子,因此统计奇数, n u m 2 num_2 num2
第二种情况 f ( n ) = 3 f(n)=3 f(n)=3
用总数减去以上算出得所有个数得到 n u m 3 num_3 num3
最终答案
a n s = 2 n u m 1 ∗ 3 n u m 3 ∗ 4 n u m 2 ans=2^{num_1}*3^{num_3}*4^{num_2} ans=2num13num34num2

#include<iostream>
using namespace std;
const unsigned long long mod=1000000007;
unsigned long long n=0;
unsigned long long num1=0,num2=0,num3=0;
unsigned long long t=0;
unsigned long long gcd(unsigned long long a,unsigned long long b){
	if(b==0) return a;
	return gcd(b,a%b);
}
unsigned long long lcm(unsigned long long a,unsigned long long b){
	return a*b/gcd(a,b);
}
unsigned long long pow(unsigned long long a,unsigned long long b){
	unsigned long long ret=1%mod;
	while(b){
		if(b&1) ret=ret*a%mod;
		a=a*a%mod;b>>=1;
	}
	return ret;
}
int main(){
	ios::sync_with_stdio(false);
	cin>>n;unsigned long long mi=2;
	num1=n-n/2-1;
	for(int i=1;i<=20;i++){
		mi*=2;
		if(mi>n) break;
		unsigned long long temp=1;
		for(unsigned long long j=1;j<mi&&temp<=n;j++) temp=lcm(temp,j);
		if(temp>n) break;
		t=n/temp;
		num2+=t-t/2;
	}
	num3=n-2-num1-num2;
	unsigned long long ret=pow(2,num1)*pow(3,num3)%mod*pow(4,num2)%mod;
	cout<<ret<<endl;
	return 0;
} 
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值