P1445-[Violet]樱花【数学】
生活随笔
收集整理的這篇文章主要介紹了
P1445-[Violet]樱花【数学】
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
正題
題目鏈接:https://www.luogu.org/problemnew/show/P1445
題目大意
求有多少個x,yx,yx,y滿足
1x+1y=1n!\frac{1}{x}+\frac{1}{y}=\frac{1}{n!}x1?+y1?=n!1?
解題思路
x+yxy=1n!\frac{x+y}{xy}=\frac{1}{n!}xyx+y?=n!1?
n!(x+y)=xyn!(x+y)=xyn!(x+y)=xy
xy?n!(x+y)=0xy-n!(x+y)=0xy?n!(x+y)=0
xy?n!(x+y)+(n!)2=(n!)2xy-n!(x+y)+(n!)^2=(n!)^2xy?n!(x+y)+(n!)2=(n!)2
(x?n!)(y?n!)=(n!)2(x-n!)(y-n!)=(n!)^2(x?n!)(y?n!)=(n!)2
然后設a=x?n!,b=y?n!a=x-n!,b=y-n!a=x?n!,b=y?n!
所以求有多少個aaa和bbb就對了
然后n!=∏picin!=\prod p_i^{c_i}n!=∏pici??
(n!)2=∏pi2?ci(n!)^2=\prod p_i^{2*c_i}(n!)2=∏pi2?ci??
然后計算aaa個數就可以推出bbb了
答案就是∏(2?ci+1)\prod (2*c_i+1)∏(2?ci?+1)
codecodecode
#include<cstdio> #include<algorithm> #include<cstring> #define ll long long using namespace std; const ll XJQ=1e9+7,N=1e6+10; ll n,ans=1,pri[N],cnt,c[N]; bool v[N]; void Prime(ll x) {for(ll i=2;i<=x;i++){if(v[i]) continue;pri[++cnt]=i;for(ll j=i;j<=x;j+=i)v[j]=1;} } void Primes(ll n) {ll l=0;for(ll i=1;i<=cnt;i++){if(pri[i]*pri[i]>n) break;if(n%pri[i]) continue;while(n%pri[i]==0)n/=pri[i],c[i]++;}if(n==1) return;int mark=lower_bound(pri+1,pri+1+cnt,n)-pri;c[mark]++; } int main() {scanf("%lld",&n);Prime(n);for(ll i=2;i<=n;i++)Primes(i);for(ll i=1;i<=cnt;i++)(ans*=c[i]*2+1)%=XJQ;printf("%lld",ans); }總結
以上是生活随笔為你收集整理的P1445-[Violet]樱花【数学】的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 如何在网上寻找异地律师协作
- 下一篇: 人们从哪种动物身上得到启发发明了防毒面具