Last non-zero Digit in N!（HDU-1066）

---

1 题意

  求$n!$最右边的非零数字。
  链接：link

---

2 思路

  已知$n!=1\times 2\times 3\times \dots \times (n-2)\times (n-1)\times n$。
  以$5$为长度对$n!$划分，$(1\times 2\times 3\times 4\times 5)\times (6\times 7\times 8\times 9\times 10)\times (11\times 12\times 13\times 14\times 15)\times \dots \times r$。
  每一块去掉一个$2$和一个$5$，对结果没有影响。
$$\frac{(1\times 2\times 3\times 4\times 1)}{2}\times \frac{(6\times 7\times 8\times 9\times 2)}{2}\times \frac{(11\times 12\times 13\times 14\times 3)}{2}\times \dots \times r$$
  取出每一组的最后一个数字。
$$\lfloor \frac{n}{5}\rfloor !\times \frac{(1\times 2\times 3\times 4)}{2}\times \frac{(6\times 7\times 8\times 9)}{2}\times \frac{(11\times 12\times 13\times 14)}{2}\times \dots \times \mathrm{r}$$
  此时，每一组对$10$取余都等于$2$，可将上式化简。
$$\lfloor \frac{n}{5}\rfloor !\times 2^{\lfloor \frac{n}{5}\rfloor }\times r$$
  每一组数中只有最后一个数能被$5$整除，所以$r$一定不能被$10$整除，可以直接对其取余。
  由于$n\%10=5$和$n\%10=0$没有余项，余项值暂定为$1$。$r\%10$以$10$为循环节，$n\%10$与$r\%10$的关系如下表所示。

$n\%10$	$r\%10$
0	1
1	1
2	2
3	6
4	4
5	1
6	6
7	2
8	6
9	4

  $2^n\%10$以$4$为循环节，$n\ge 1$。$n\%4$与$2^n\%10$的关系如下表所示。

$n\%4$	$2^n\%10$
0	6
1	2
2	4
3	6

  因为$r\%10$以$10$位循环节，$2^n\%10$以$4$为循环节，所以$(r\times 2^n)\%10$以$20$，其中$n\ge 1$。$n\%20$与$(r\times 2^n)\%10$的关系如下表所示。

$n\%20$	$(r\times 2^n)\%10$
0	1
1	1
2	2
3	6
4	4
5	2
6	2
7	4
8	2
9	8
10	4
11	4
12	8
13	4
14	6
15	8
16	8
17	6
18	8
19	2

  当$n=0$时$(r\times 2^n)\%10=1$，也满足上述表格。
  事先将$(r\times 2^n)\%10$对$n\%20$的映射值存储在$m[]$数组中，则$\mathrm{L}\mathrm{N}\mathrm{D}(\mathrm{n})=\mathrm{L}\mathrm{N}\mathrm{D}\left(\lfloor \frac{n}{5}\rfloor \right)\times m[\mathrm{n}\%20]\%10$。

2.1 时间复杂度分析

  根据递推式可知，递归或循环求解只需要$\mathcal{O}(\log (n))$的时间，由于除$5$也需要$\mathcal{O}(\log (n))$的时间，所以时间复杂度为$\mathcal{O}({\log }^2(n))$。

2.2 实现

2.2.1 递归实现

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int N=1e3+10;
int m[20]={1,1,2,6,4,2,2,4,2,8,4,4,8,4,6,8,8,6,8,2};
char buf[N],bufLen;
int num[N],numLen;
int lastdigit(){
    if(numLen==0) return 1;
    int res=m[num[1]%2*10+num[0]]%10;//对20取余
    for(int c=0,i=numLen-1;i>=0;i--){
        c=c*10+num[i],num[i]=c/5,c%=5;//除5
    }
    if(!num[numLen-1]) numLen--;
    return lastdigit()*res%10;
}
int main(){
    while(~scanf("%s",buf)){
        bufLen=strlen(buf);
        for(int i=0;i<bufLen;i++) num[i]=buf[bufLen-1-i]-'0';
        numLen=bufLen;
        printf("%d\n",lastdigit());
    }
    return 0;
}

2.2.2 非递归实现

#include<cstdio>
#include<cstring>
using namespace std;
const int N=1e5+10;
char str[N],a[N];
int m[20]={1,1,2,6,4,2,2,4,2,8,4,4,8,4,6,8,8,6,8,2};
int lastdigit(char *buf){
    int len=strlen(buf),ret=1;
    if(len==1) return m[buf[0]-'0'];
    for(int i=0;i<len;i++) a[i]=buf[len-1-i]-'0';
    for(;len;len-=!a[len-1]){
        ret=ret*m[a[1]%2*10+a[0]]%5;//对20取余
        for(int c=0,i=len-1;i>=0;i--){
            c=c*10+a[i],a[i]=c/5,c%=5;//除5
        }
    }
    return ret+ret%2*5;
}
int main(){
    while(~scanf("%s",str)){
        printf("%d\n",lastdigit(str));
    }
    return 0;
}