[NOIP2017模拟]Factorial Surplus Tail

题目描述
FST 作为 OIer ，经常会遇到和阶乘有关的问题，但是一个数的阶乘末尾总是会有很多 0 ，FST 认为这很不美观，但是 FST 觉得如果 0 的个数是偶数的话，还是可以接受的。

所以就有这样一个问题，FST 想知道 0！，1！，2！… … (n-1)!，n! 中有多少数的末尾 0 个数是偶数。（注意0！是1，0算偶数）

输入格式
读入有若干行，每行一个正整数 n ，最后一行是一个 -1 。

输出格式
对于每个 n 输出一行，为 0！，1！，2！… … (n-1)! ，n! 中末尾 0 个数是偶数的个数。

样例数据
输入
2
3
10
-1
输出
3
4
6

备注
【数据范围】
测试点 1、2：n≤10；数据组数=1；
测试点 3、4：n≤10000；数据组数=10；
测试点 5、6、7、8：n≤$10^9$；数据组数=$10^5$；
测试点 9、10：n≤$10^{18}$；数据组数=$10^5$；

分析：
算法一： 最暴力的暴力 20%
直接把阶乘算出来数 0 的个数。

算法二： 较优美的暴力 40%
因为质因子 2 肯定比质因子 5 多（这个就不需要证明了吧） ，所以相当于统计质因
子 5 的个数，所以暴力可以写地更简单。

算法三： 分段打表 60%
我们可以每隔 100W 个数打个表，这样暴力复杂度为 O($10^6$)。

算法四： 类 DP 算法 100%
第一次 DP，预处理计算[0，$5^k$）的答案和 0 个数的奇偶变化，可以从 $5^{k-1}$
转移。
可以观察下规律： （0 表示偶数，1 表示奇数）
k=0： 0
k=1： 00000
k=2： 0000011111000001111100000
k=3： 0000011111000001111100000000001111100000… …
有 2 种可能：
k是偶数， 那么 k-1是奇数， 0 个数的奇偶会变， $ans_k=3*ans_{k-1}+2*(5^{k-1}-ans_{k-1})$
k是奇数，那么 k-1是偶数，0 个数的奇偶不会变，$ans_k=ans_{k-1}*5$

第二次 DP，把 n+1 转化为 5 进制，从高位开始做。
然后我们可以分类讨论，例如当前位的数是 3，当前位是奇数位，现在的奇偶性为偶，那么转移答案为 $ans_k*2+(5^k-ans_k)*1$，因为我们预处理的答案相当于从偶数开始，然后奇数位会改奇偶性，分配给 3 段就是($ans_k$) ($ans_k$取反) ($ans_k$)。
（可以不管它，我都看不懂这是什么鸟语，看下面模拟可能还好一点）
每次复杂度 O($log_5n$)

代码：
这个题不好理解，更不好解释，代码下方我会模拟一次

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<ctime>
#include<cmath>
#include<algorithm>
#include<cctype>
#include<iomanip>
#include<queue>
#include<set>
using namespace std;

long long getlong()
{
    long long sum=0,f=1;
    char ch;
    for(ch=getchar();(ch<'0'||ch>'9')&&ch!='-';ch=getchar());
    if(ch=='-')
    {
        f=-1;
        ch=getchar();
    }
    for(;ch>='0'&&ch<='9';ch=getchar())
        sum=(sum<<3)+(sum<<1)+ch-48;
    return sum*f;
}

long long n,mi[30],ans,num[30][2],tmp;

int main()
{
    freopen("fstagain.in","r",stdin);
    freopen("fstagain.out","w",stdout);

    mi[0]=1;
    for(int i=1;i<=25;++i)//预处理5的次方
        mi[i]=mi[i-1]*5;

    num[0][0]=1;//num[i][0/1]表示[0,5^i)中有多少个0/1
    for(int i=1;i<=25;++i)//根据规律预处理
        if(i%2==1)
        {
            num[i][0]=num[i-1][0]*5;
            num[i][1]=num[i-1][1]*5;
        }
        else
        {
            num[i][0]=3*num[i-1][0]+2*num[i-1][1];
            num[i][1]=3*num[i-1][1]+2*num[i-1][0];
        }

    n=getlong();
    while(n!=-1)
    {
        ans=0;
        tmp=0;
        for(int i=25;i>=0;--i)//第二次dp，结合下方模拟理解
        {
            for(int j=1;j<=n/mi[i];++j)
            {
                ans+=num[i][tmp];
                if(i%2==1)
                    tmp^=1;
            }

            n%=mi[i];
        }

        cout<<ans+num[0][tmp]<<endl;//因为我们的区间左闭右开，是[0,n-1]，还差一个数，最后加上
        n=getlong();
    }

    return 0;
}

模拟
例如，现在的n为74（五进制244，是为了方便看加几次满足的），直接进入第二个dp，i=2时，加上两个$5^2$时满足的个数；i=1时，n还余下24，i为奇数，要取正取反取正取反；i=0时，n还余下4，i为偶数，i=1时刚好取反于是到这一步是取正了，加上四个1。最后还有一个数，由于i=0时没有取正取反，所以直接单独再加一个1。

本题结。