[Violet]樱花/阶乘分解
又到了一年樱花盛开的时节。Vani 和妹子一起去看樱花的时候,找到了一棵大大的樱 花树,上面开满了粉红色的樱花。Vani 粗略估计了一下,一共有足足 ! n 片花瓣。
Vani 轻柔地对她说:“你知道吗?这里面的一片花瓣代表着你,我从里面随机摘一片,能和你相遇的概率只有1/n!那么小。我该是多么的幸运,才让你今天这么近地站在我面前。
相信我,我一定会把这亿万分之一的缘分变为永远。”
粉红的樱花漫天飞舞,妹子瞬间被 Vani 感动了。她轻轻地牵起了他的手,和他相依而 坐。这时,她突然看到田野的尽头也长着两棵樱花树,于是慢慢地把头靠在 Vani 的肩上,在他耳边低语:“看到夕阳里的那两棵樱花树了吗?其中一棵树上的一片花瓣是你,另一棵树上的一片花瓣是我,如果有人从这棵摘下一片,从那棵采下一瓣,我们相遇的概率会不会正好是1/n!呢?”
Vani 的大脑飞速运作了一下,立即算出了答案。正要告诉妹子,她突然又轻轻地说:“以 前你总是说我数学不好,但是这种简单的题我还是会算的。你看假如左边那棵树上有 x片花瓣,右边那个有 y 片花瓣,那么我们相遇的概率不就是1/x+1/y么,不过有多少种情况能使
它正好可以等于1/n!呢?这个你就帮我算一下吧~”
显然,面对天然呆的可爱妹子,Vani 不但不能吐槽她的渣数学,而且还要老老实实地 帮她算出答案哦。
题目描述
求不定方程
1/x+1/y=1/n!
的正整数解(x,y)的数目。
输入格式
一个整数n。
输出格式
一个整数,表示有多少对 (x,y) 满足题意。答案对10^9+7取模。
秒速五厘米,那是樱花飘落的速度,那么怎样的速度,才能走完我与你的距离
1.化简可得(x+y)*(n!)=xy; [式子1]
根据原式我们知道x,y>(n!)
设x=(n!)+c;y=(n!)+d;
代入式子1得
c*d=(n!)*(n!)
也就是说我们要求(n!)*(n!)的约数个数
2.我们先用线性筛一遍质数
根据算术基本定理:每个大于1的自然数均可写为质数的积
那么他的约数个数就是指选取其中任意个质因数组合,论有多少种方法
#include<bits/stdc++.h> #define re return #define ll long long #define inc(i,l,r) for(ll i=l;i<=r;++i) using namespace std; template<typename T>inline void rd(T&x) { char c;bool f=0; while((c=getchar())<'0'||c>'9')if(c=='-')f=1; x=c^48; while((c=getchar())>='0'&&c<='9')x=x*10+(c^48); if(f)x=-x; } const int maxn=1000001; ll n,ans=1,num; int cnt[maxn],notprime[1000005],prime[1000005]; inline void Prime() { inc(i,2,n) { if(!notprime[i]) { prime[++num]=i; } for(int j=1;j<=num;++j) { if(prime[j]*i>n)break; notprime[prime[j]*i]=1; if(i%prime[j]==0)break; } } } int main() { rd(n); Prime(); inc(i,1,num) { for(int j=prime[i];j<=n;j+=prime[i]) { int x=j; while(x%prime[i]==0) { ++cnt[i]; x/=prime[i]; } } cnt[i]<<=1; } inc(i,1,num) ans=(ans*(cnt[i]+1))%1000000007;
//加一的原因是因为可以选1~cnt[i]个,也可以不选 printf("%d",ans); re 0; }
正确的数学姿态是:我们发现N!中质数因子p的个数,就是1~N中每个数含有的质因数p个数.既然如此的话,那么我们发现,至少有一个质因子p的显然有n/p个,而至少有两个质因子p数的显然是有(n/p)/p个