两道阶乘题做法

其实蛮简单的。。。

【fac】

求n!末尾第一位非0数。

【pku1150】

求p(n,m)末尾第一位非0数。

n非常大。。。

第一题n达到10^100次方。。。这个显然无法模拟。。。如果直接做无论空间或者时间都无法支持。。

第二题n达到2000000，显然也是无法模拟的

我们注意到 零必然是由2和5支持的、5的个数比2的个数少

同时

1 2 3 4 5 6 7 8 9 10 11 12 13 14

f（n）为n!末尾第一位非0数。

我们注意到 1*2*3*4=24 末尾为4.。。能被4整初。。提出一个2支持5.

同时 6*7*8*9 末尾也为4 、、 能被4乘除。。提出一个2支持5.

n提供了n/5个4 和n/5个5.。。如果mod 5<>0其余的特殊处理一下

每个4提供一个2给5

一共有2^(n/5)...然后5的倍数提出5后变成 1 2  3 4 5 6 7 8 9、、

f(n)=2^（n/5)*1*2*3*4*5、、、*p【n mod 5)】（p为特殊处理部分）

f(n)=2^(n/5)*f(n/5)

递归处理即可

n除以5要高精。。其实可以乘2除以10就可以。。节约代码。、、

pku1150

这题有点不一样。。无法直接除。。必需换其他方法做。。。

我们可以将2，5提出、、末尾非0数就是1 3 7 9个数决定的了这个只要知道个数就很好求了。。

2和5的个数我们可以在log（n）时间内求出来

首先1 2 3 4 5 6 7 8 。将其分成奇数偶数组。。。

2 4 6 8 10.

1 3 5 7 9、、、、、、

将偶数组2 4 6 8 10. 12 的2提出后为1 2 3 4 5 6 7.。。。。。

f（n）:=f(n/2)+g(n)

g（n）表示n中几个提出2和5的后尾数为k。。

1 3 5 7 9 11 13 15 17 19 21 23 25.。。

我们注意到其中还有5（5，15，25，35.。。。））

所以再将5提出。。

变为1 3 5 7 9 11 13.。。。。依次递归下去

g(n)=g(n/5)+n mod 10+ord（k>=n mod 10)

p(n,m)=n!/(n-m)!

将2，3，5，7，9为尾数的个数求出来。。然后相减。。根据同余求出尾数。

{$m 100000000}//bzoj1308 type row=array[0..101]of integer; const p:array[0..9]of integer=(1,1,2,6,4,1,6,2,6,4); var a:row; s:string; i,tot:longint; f2:Array[1..4]of integer; function dv(a:Row):row; var i:longint; begin fillchar(dv,sizeof(dv),0); for i:=1 to a[0] do dv[i]:=a[i]+a[i]; for i:=1 to a[0] do begin dv[i+1]:=dv[i+1]+dv[i] div 10; dv[i]:=dv[i] mod 10; end; if dv[a[0]+1]>0 then dv[0]:=a[0]+1 else dv[0]:=a[0]; for i:=1 to dv[0] do dv[i]:=dv[i+1]; dv[dv[0]]:=0; dec(dv[0]); end; function work(a:row):longint; var b:row; x:integer; begin if (a[1]=0)and(a[0]=0) then exit(1); b:=dv(a); if (b[1]=0)And(b[0]=0) then x:=1 else x:=f2[(b[2]*10+b[1]) mod 4]; work:=(work(b)*x*p[a[1]])mod 10; end; begin while not seekeof do begin readln(s); if s='0' then begin writeln(0); continue; end; fillchar(a,sizeof(a),0); f2[0]:=6;f2[1]:=2;f2[2]:=4;f2[3]:=8; tot:=0; for i:=length(s) downto 1 do begin inc(tot); a[tot]:=ord(s[i])-48; end; a[0]:=length(s); writeln(work(a)); end; end.

//pku1150 var f2,f3,f7,f9:Array[0..4]of longint; n,m,a3,a5,a2,a7,a9:int64; i:longint; function count(x,k:int64):int64; begin count:=0; while x>0 do begin inc(count,x div k); x:=x div k; end; end; function g(X,k:int64):int64; begin if x=0 then exit(0); g:=x div 10+ord(x mod 10>=k)+g(X div 5,k); end; function work(X,k:int64):int64; var now:int64; begin if x=0 then exit(0); work:=work(x div 2,k)+g(x,k); end; begin f2[0]:=1;f3[0]:=1;f7[0]:=1;f9[0]:=1; for i:=1 to 4 do begin f2[i]:=f2[i-1]*2 mod 10; f3[i]:=f3[i-1]*3 mod 10; f7[i]:=f7[i-1]*7 mod 10; f9[i]:=f9[i-1]*9 mod 10; end; f2[0]:=6; while not seekeof do begin readln(n,m); a2:=count(n,2)-count(n-m,2); a5:=count(n,5)-count(n-m,5); a3:=work(n,3)-work(n-m,3); a7:=work(n,7)-work(n-m,7); a9:=work(n,9)-work(n-m,9); if a5>a2 then writeln(5) else if a2=a5 then writeln(f9[a9 mod 4]*f7[a7 mod 4]*f3[a3 mod 4] mod 10) else writeln(f9[a9 mod 4]*f7[a7 mod 4]*f2[(a2-a5) mod 4]*f3[a3 mod 4] mod 10); end; end.