bzoj3209:3209: 花神的数论题

觉得还是数位dp的那种解题形式但是没有认真的想,一下子就看题解。其实还是设置状态转移。一定要多思考啊
f[i][j]=f[i-1][j]+g[i-1][j] g[i][j]=f[i-1][j-1]+g[i-1][j];

然后我就开始gang。然后先是for for j没有从0开始。然后是cnt增加的时候忘了*,接着是1<<tmp没有用ll,还有是读入优化没有用longlong,最后成功的过了若干较小的数据。WAWAWAWAWA5发。最后发现指数不能直接取模!。终于AC了喜极而泣TAT

=>关键是前两个错误思考了好久输出调试了好久渐渐的才知道自己哪里写错了。妈呀以后一遇到数位dp就把过程什么的全都输出一下。。。有关的东西打表出来确认每一步都没有错。。。

=>不过用了两个小时调试收获还是很大的。

#include<cstdio>
#include<cstring>
#include<cctype>
#include<algorithm>
using namespace std;
#define rep(i,s,t) for(int i=s;i<=t;i++)
#define dwn(i,s,t) for(int i=s;i>=t;i--)
#define clr(x,c) memset(x,c,sizeof(x))
#define ll long long
ll read(){ //注意读入不要忘了longlong
	ll x=0;char c=getchar();
	while(!isdigit(c)) c=getchar();
	while(isdigit(c)) x=x*10+c-‘0‘,c=getchar();
	return x;
}
const int mod=1e7+7;
const int nmax=65;
ll f[nmax][nmax],g[nmax][nmax];
ll pw(ll a,ll b){
	ll ans=1;
	while(b){
		if(b&1) ans=ans*a%mod;
		b>>=1;a=a*a%mod;
	}
	return ans;
}
ll cal(ll x){//10011
	ll tmp,temp,cnt=0,ans=1;
	for(tmp=0;(1ll<<tmp)<=x;++tmp);//如果直接1<<tmp会爆
	dwn(i,tmp,1){
		temp=1ll<<(i-1);
		if(x&temp) {
			rep(j,1,i-1) ans=ans*pw(j+cnt,f[i][j])%mod;
			if(cnt) ans=ans*cnt%mod;++cnt;
		}
	}
	return ans*cnt%mod;
}
int main(){
	ll n=read();
	f[1][0]=g[1][1]=1;
	rep(i,2,60) rep(j,0,i){//指数不能直接取模!!!
		 f[i][j]=f[i-1][j]+g[i-1][j];
		 if(j) g[i][j]=f[i-1][j-1]+g[i-1][j-1];
	}
	printf("%lld\n",cal(n));
	return 0;
}

  

3209: 花神的数论题

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 1802  Solved: 822
[Submit][Status][Discuss]

Description

背景
众所周知,花神多年来凭借无边的神力狂虐各大 OJ、OI、CF、TC …… 当然也包括 CH 啦。
描述
话说花神这天又来讲课了。课后照例有超级难的神题啦…… 我等蒟蒻又遭殃了。
花神的题目是这样的
设 sum(i) 表示 i 的二进制表示中 1 的个数。给出一个正整数 N ,花神要问你
派(Sum(i)),也就是 sum(1)—sum(N) 的乘积。

Input

一个正整数 N。

Output

一个数,答案模 10000007 的值。

Sample Input

样例输入一

3

Sample Output

样例输出一

2

HINT

对于样例一,1*1*2=2;

数据范围与约定

对于 100% 的数据,N≤10^15

Source

原创 Memphis

时间: 2024-10-05 14:41:14

bzoj3209:3209: 花神的数论题的相关文章

HYSBZ 3209: 花神的数论题

题目连接:http://www.lydsy.com/JudgeOnline/problem.php?id=3209 花神的数论题 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1148  Solved: 535[Submit][Status][Discuss] Description 背景众所周知,花神多年来凭借无边的神力狂虐各大 OJ.OI.CF.TC …… 当然也包括 CH 啦.描述话说花神这天又来讲课了.课后照例有超级难的神题啦…… 我等蒟

HYSBZ 3209 花神的数论题(数位dp)

3209: 花神的数论题 Time Limit: 10 Sec  Memory Limit: 128 MB Submit: 1148  Solved: 535 [Submit][Status][Discuss] Description 背景 众所周知,花神多年来凭借无边的神力狂虐各大 OJ.OI.CF.TC -- 当然也包括 CH 啦. 描述 话说花神这天又来讲课了.课后照例有超级难的神题啦-- 我等蒟蒻又遭殃了. 花神的题目是这样的 设 sum(i) 表示 i 的二进制表示中 1 的个数.给出

BZOJ:3209: 花神的数论题

今天居然没参考任何资料解决了这道数位DP,事先只是搞一道数论题练练: 思路:求SUM[1]-SUM[N]的二进制的乘积mod1000000007: N<=10^15;BZOJ的题不会是几个简单的FOR就完结的, 假如N的二进制是:1001001: 位数:1234567 那么当第一位是0的情况:那么数就是000000-111111 序列号:  123456-123456 发现我们可以用排列组合算出二进制有1-6个1的数比如:1个1就是C[6][1],2个1就是C[6][2],...... 当第一位

bzoj 3209: 花神的数论题 数位dp

3209: 花神的数论题 Time Limit: 10 Sec  Memory Limit: 128 MB[Submit][Status][Discuss] Description 背景众所周知,花神多年来凭借无边的神力狂虐各大 OJ.OI.CF.TC …… 当然也包括 CH 啦.描述话说花神这天又来讲课了.课后照例有超级难的神题啦…… 我等蒟蒻又遭殃了.花神的题目是这样的设 sum(i) 表示 i 的二进制表示中 1 的个数.给出一个正整数 N ,花神要问你派(Sum(i)),也就是 sum(

【BZOJ3209】花神的数论题 数位DP

[BZOJ3209]花神的数论题 Description 背景众所周知,花神多年来凭借无边的神力狂虐各大 OJ.OI.CF.TC …… 当然也包括 CH 啦.描述话说花神这天又来讲课了.课后照例有超级难的神题啦…… 我等蒟蒻又遭殃了.花神的题目是这样的设 sum(i) 表示 i 的二进制表示中 1 的个数.给出一个正整数 N ,花神要问你派(Sum(i)),也就是 sum(1)—sum(N) 的乘积. Input 一个正整数 N. Output 一个数,答案模 10000007 的值. Samp

【BZOJ】3209: 花神的数论题

题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=3209 显然是按照二进制位进行DP. 考虑预处理$F[i][j]$表示到了二进制的第$i$位,有$j$个$1$的数字有多少个. 显然:${F[i][j]=F[i-1][j-1]+F[i-1][j]}$ 组合数... 接下来只需补充不漏的计算比$n+1$小的每一个数字对应的1的多少.数位统计即可. 1 llg work(llg x) 2 { 3 llg tot=0; 4 for (llg i

BZOJ 3209: 花神的数论题【数位dp】

Description 背景众所周知,花神多年来凭借无边的神力狂虐各大 OJ.OI.CF.TC …… 当然也包括 CH 啦.描述话说花神这天又来讲课了.课后照例有超级难的神题啦…… 我等蒟蒻又遭殃了.花神的题目是这样的设 sum(i) 表示 i 的二进制表示中 1 的个数.给出一个正整数 N ,花神要问你派(Sum(i)),也就是 sum(1)—sum(N) 的乘积. Input 一个正整数 N. Output 一个数,答案模 10000007 的值. Sample Input 样例输入一3 S

【bzoj3209】 花神的数论题

http://www.lydsy.com/JudgeOnline/problem.php?id=3209 (题目链接) 题意 ${sum(i)}$表示${i}$的二进制表示中${1}$的个数.求${\prod^n sum(i)}$ Solution ${f_{i,s}}$表示dp到第${i}$位,已经有${s}$个${1}$时的乘积.然后一路dfs就可以了. 细节 LL,返回值要与1取个max 代码 // bzoj3598 #include<algorithm> #include<ios

BZOJ 3209 花神的数论题 数位DP+数论

题目大意:令Sum(i)为i在二进制下1的个数 求∏(1<=i<=n)Sum(i) 一道很简单的数位DP 首先我们打表打出组合数 然后利用数位DP统计出二进制下1的个数为x的数的数量 最后输出∏(1<=x<=logn)x^ans[x]即可 此题的坑在于这题的组合数和数位DP的结果都是指数 对指数取模不能直接取 要取Phi(p) 于是我们对10000006取模 然后这题就WA了 因为10000007不是个质数! 10000007=941*10627 于是我们得到Phi(p)=940*