数论——[Violet]樱花

Posted 71-111

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了数论——[Violet]樱花相关的知识,希望对你有一定的参考价值。

题意:求方程 1/X+1/Y=1/(N!) 的正整数解的组数,其中N≤10^6。


 

设n!=z,y=z+d

1/x+1/y=1/z

1/x+1/(z+d)=1/z

(x+z+d)/(x*z+dx)=1/z

z(x+z+d)=x*z+dx

z^2+dz=dx

x=z^2/d+z

所以只要求有z^2有多少约数(即d的集合元素个数)即可


 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 #define N 1000003
 4 #define mo 1000000007
 5 #define ll long long
 6 using namespace std;
 7 int n,pri[N],mp[N],num[N],p[N];
 8 
 9 void oula()
10 {
11     for (int i=2;i<=N;i++)
12     {
13         if (!p[i]) pri[++pri[0]]=i,mp[i]=pri[0];
14         for (int j=1;i*pri[j]<=N&&j<=pri[0];j++)
15         {
16             p[i*pri[j]]=1;
17             if (i%pri[j]==0) break;
18         }
19     }
20 }
21 void cal(int x)
22 {
23     for (int i=1;pri[i]*pri[i]<=x;i++)
24     {
25         if (x%pri[i]==0)
26         {
27             while (x%pri[i]==0) num[i]++,x/=pri[i];
28         }
29         
30     }
31     if (x>1) num[mp[x]]++;
32 }
33 int main()
34 {
35     scanf("%d",&n);
36     oula();
37     for (int i=1;i<=n;i++)
38         cal(i);
39     ll ans=1;
40     for (int i=1;i<=pri[0];i++) ans=ans*(ll)(num[i]*2+1)%mo;//计算 
41     printf("%lld",ans);
42     return 0;
43 }

 

以上是关于数论——[Violet]樱花的主要内容,如果未能解决你的问题,请参考以下文章

bzoj2721[Violet 5]樱花 数论

BZOJ_2721_[Violet 5]樱花_数学

Luogu1445 [Violet]樱花

P4167 [Violet]樱花

2721: [Violet 5]樱花|约数个数

P1445 [Violet]樱花(质因数分解)