LightOJ-1341 Aladdin and the Flying Carpet 因子個數
阿新 • • 發佈:2018-12-17
這道題要用到一個結論:
設 N=p1^a1*p2^a2*p3^a3*...*pn^an(其中p1,p2,...,pn為N的因子,a1,a2,... ,an分別為因子的指數)(唯一分解定理);則 N的因子個數 M=(1 + a1)*(1 + a2)*(1 + a3)*...*(1 + an)
然後M/=2,再減去不合條件的個數即可。
附上AC程式碼:
#include<cstdio> #include<cstring> #include<cmath> #include<iostream> #include<algorithm> using namespace std; #define ll long long const ll MAX=1e6+5; bool isp[MAX]; ll pri[MAX]; void init() { memset(isp,true,sizeof(isp)); memset(pri,0,sizeof(pri)); isp[0]=isp[1]=false; for(ll i=2;i<MAX;i++) if(isp[i]) { pri[++pri[0]]=i; for(ll j=2;i*j<MAX;j++) isp[i*j]=false; } } ll fac[MAX],tot;//質因子 ll num[MAX];//每個質因子的冪 void decom(ll a) { tot=0; memset(num,0,sizeof(num)); ll n=a; for(ll j=1;j<=pri[0];j++) { if(pri[j]*pri[j]>n) break; if(n%pri[j]==0) { fac[tot]=pri[j]; while(n%pri[j]==0) { num[tot]++; n/=pri[j]; } tot++; } } if(n>1) { num[tot]++; fac[tot++]=n; } /*for(int i=0;i<tot;i++) cout<<"i="<<i<<" fac="<<fac[i]<<" num="<<num[i]<<endl;*/ } ll a,b; int main() { int t; scanf("%d",&t); init(); for(int tt=1;tt<=t;tt++) { scanf("%lld%lld",&a,&b); if(a<b*b) { printf("Case %d: 0\n",tt); continue; } decom(a); ll sum=1; for(ll i=0;i<tot;i++) sum*=(1+num[i]); sum/=2; //cout<<"sum="<<sum<<endl; for(ll i=1;i<b;i++)//此處不能用fac判斷 eg:36,4 if(a%i==0) sum--; printf("Case %d: %lld\n",tt,sum); } return 0; }