【USACO 3.2.1】階乘

標籤：style   class   blog   code   color   使用   

【描述】

N的階乘寫作N!表示小於等於N的所有正整數的乘積。階乘會很快的變大，如13!就必須用32位整數類型來儲存，70！即使用浮點數也存不下了。你的任務是找到階乘最後面的非零位。舉個例子,5!=1*2*3*4*5=120所以5!的最後面的非零位是2，7！=1*2*3*4*5*6*7=5040，所以最後面的非零位是4。

  【格式】PROGRAM NAME: rect4INPUT FORMAT:(file rect4.in) 

共一行，一個整數不大於4，220的整數N。

 OUTPUT FORMAT:(file rect4.out)

 

共一行，輸出N!最後面的非零位。

【分析】

先打素數表，然後分解質因數，再去掉可以產生末尾0的2和5(不是所有二和五，而是他們的最小值)，最後在把剩下的乘起來MOD10。

 1 #include <iostream> 2 #include <cstdio> 3 #include <cmath> 4 #include <cstring> 5 #include <algorithm> 6 const int maxn=4220; 7 using namespace std;  8 int cnt[maxn],vis[maxn]; 9 int prime[maxn];10 void built_prime();//打表 11 void check(int n);12 int pow_mod(int num,int t);13 int main()14 {15     int n,i,temp,ans=1;16     //檔案操作17     freopen("fact4.in","r",stdin);18     freopen("fact4.out","w",stdout);19     built_prime();20     //printf("1");21     scanf("%d",&n);22     for (i=1;i<=n;i++) check(i);//分解質因數 23     temp=min(cnt[2],cnt[5]);24     cnt[2]-=temp;cnt[5]-=temp;25     for (i=1;i<=prime[0];i++)26     {27         if (cnt[prime[i]]!=0)28         ans=(ans*pow_mod(prime[i],cnt[prime[i]]))%10;//快速冪模數 29     }30     printf("%d",ans%10);31     return 0;32 }33 int pow_mod(int num,int t)34 {35     if (t==1) return num%10;36     if (t==2) return (num*num)%10;37     int temp=pow_mod(num,t/2);38     if (t%2==0) return (temp*temp)%10;39     else return (temp*temp*num)%10;40 }41 void built_prime()42 {43      int i,j;44      memset(cnt,0,sizeof(cnt));45      prime[0]=0;46      for (i=2;i<=maxn;i++)47      {48          if (vis[i]==0)49          {50              prime[++prime[0]]=i;51              for (j=i;j<=maxn;j+=i) vis[j]=1;52          }53      }54 }55 void check(int n)56 {57      int i;58      while (n!=1)59      {60            for (i=1;i<=prime[0];i++)61            if (n%prime[i]==0) {cnt[prime[i]]++;n=n/prime[i];break;}62      }63 }