yzoj P1122 阶乘 题解

T组数据,给出N,求出N!最右边非零的数。

对于30%的数据,N <= 30,T<=10。

对于全部的数据,N <= 10^2009,T<=30。

一道数学题

解析

N!/(10^x)最后一位数字即是结果。10^x进行拆分,变成5^x* 2^x。怎么除以5^x呢,好办,乘的时候含有5的倍数的一项全部不乘进去,再递归此过程。即

1 2 3 4 (15) 6 7 8 9 (25)

11 12 13 14 (35)16 17 18 19 (45)

21 22 23 24 (55) 26 27 28 29 (65) ...

再递归处理

1 2 3 4 (1*5) 6

而且可以发现,不算5倍数一项进去,每十个一组,最末尾结果是一样的:6,且6*6还是6!

接下来的问题,怎么除2^x呢?

分析发现:

2^0 = 1

2^1 = 2

2^2 = 4

2^3 = 8

2^4 =16=(6)

且2,4,6,8乘以6最末尾还是原来的数。所以,可以分析看x%4是多少,若为1,说明原来的结果相当于多乘了一个2(其他的那些2刚好是4的倍数,不影响结果),所以,我们本来是要除以2的,但是很显然对结果再补乘3个2(2^3)即可消除影响,得到正确答案。其他情况同理。至此,思路理顺。

代码

#include<bits/stdc++.h>
using namespace std;
int T,n,ans,mod,rest,x;
bool flag;
const int v[10]={1,1,2,6,4,4,4,8,4,6};
char s[3000];
int a[3000];
const int k[4]={6,8,4,2};
int main(){
    scanf("%d",&T);
    while(T--){
        scanf("%s",s);
        memset(a,0,sizeof(a));
        n=strlen(s);
        for(int i=1;i<=n;++i){
            a[i]=s[n-i]-'0';
        }
        ans=1;
        mod=0;
        flag=0;
        while(1){
            rest=0;
            ans=(ans*v[a[1]])%10;
            if(n>1) ans=(ans*6)%10;
            for(int i=n;i>=1;i--){
                x=rest*10+a[i];
                rest=x%5;a[i]=x/5;
            }
            while((n>0)&&(a[n]==0)) n--;
            if(n==0) break;
            flag=true;
            mod=(mod+a[2]*10+a[1])%4;
        }
        if(flag) ans=(ans*k[mod])%10;
        printf("%d\n",ans);
    }
    return 0;
}

原文地址:https://www.cnblogs.com/donkey2603089141/p/11416508.html

时间: 2024-10-10 17:10:33

yzoj P1122 阶乘 题解的相关文章

yzoj P2345 战争 题解

纯数论 30分:纯暴力,直接模拟判断t秒后,判断hp是否小于0 60分: atk>=h,就是一炮一个,那么军队会在min(n,t)秒之后停止攻击,那么总伤害就是a[n+(n-1) +(n-2)+........(n-min(n,t)+1) 等差数列求和d=a(2n-min(n,t)+1)*(min(n,t))/2; 如果d>=hp Yes,否则No: 100分 考虑前t秒军队伤害总和,然后和hp比较大小 先算几下打死一个士兵,设m下 那么前m秒伤害为nma m+1到2m秒伤害为m(n-1)a

[USACO07OCT]障碍路线 &amp; yzoj P1130 拐弯 题解

题意 给出n* n 的图,A为起点,B为终点,* 为障碍,.可以行走,问最少需要拐90度的弯多少次,无法到达输出-1. 解析 思路:构造N * M * 4个点,即将原图的每个点分裂成4个点.其中点(i,j,k)表示在(i,j)时人的方向是k,然后对于两个点(i,j,k)和(i,j,kk),如果k和kk是两个旋转90度能转换的方向,就连一条边权为1的边,而对于(i,j,k)和(i+dx[ k],j+dy[k],k)连一条边权为0的边,表示从(i,j)在方向为k的情况下能向k方向走一步到达(i+dx

yzoj P2371 爬山 题解

背景 其实 Kano 曾经到过由乃?,当然这名字?看?主就是 Yuno 嘛.当年 Kano 看见了由乃?,内?突然涌出了?股杜甫会当凌绝顶,?览众??的 豪?,于是毅然决定登?.但是 Kano 总是习惯性乱丢垃圾,增重环卫?? 的负担,Yuno 并不想让 Kano 登?,于是她果断在?上设置了结界-- 题意 Yuno 为了?便登?者,在?上造了 N 个营地,编号从 0 开始.当结界发动时,每当第 $ i(> 0) $ 号营地内有?,那么他将被传送到第 $ A_i (< i) $ 号营 地,如此

【题解】阶乘因子

题目描述 桐桐刚刚学习了自然数N的阶乘:阶乘(N!)被定义成从1到N的所有整数的乘积,例如5!=5×4×3×2×1=120.随着数N的增大,N!增长的非常快,5!=120,10!=3628800.桐桐想到了一种方法来列举那么大的数:不是直接列出该数,而是按照顺序列举出该数中各个质数因子出现的次数.如825可描述为(0 1 2 0 1),意思是对825分解质因数,这些质数因子中有0个2,1个3,2个5,0个7,1个11.请你编一个程序,读入N值,帮助桐桐按顺序输出N!所包含的质数因子的个数. 输入

yzoj P2344 斯卡布罗集市 题解

共t条街对于每一条街上有n个店铺(n可能不相同),每次只能选两端第一个营业的店铺采购,采购第i个店铺会获得幸福度ai,采购完后,这个店铺和它相邻的店铺便会关门,问最大幸福度? 考场想了一下dp,一开始想一维但发现不好处理,二维参数也没有想出来,于是便开始了我的暴力瞎搞之旅,我随手写了几个例子发现对于n为奇数无论怎么采购,幸福度是固定的为a1, a3 , a5........an于是便可以直接累加,然而对于n为偶数又怎么去处理呢?我也随手写了n为6,8,10的数据发现其实就只有三种情况. 第一种情

yzoj P2349 取数 题解

题意 1到n个自然数中选k个自然数要求两两不相邻,问有多少种方法,模m eg(1 3 5 ) 又是一道打表规律题,正常解法dp可以通过前缀和优化到O(N* K).另外我们可以重新定义F[I,J]表示从1到I中选择J个不连续数的方案数.通过考虑I选还是不选来进行状态转移. 1.如果不选I:则方案数为F[I-1,J]; 2.如果选I:由于不能选相邻两个数,所以I-1不能选,则剩余的J-1个数只能在1到I-2中选,即F[I-2,J-1]; f[i][j]=f[i-1][j]+f[i-2][j-1] (

yzoj 2372 小B的数字 题解

题意 判断是否存在一个序列 $ b_i $ 使得 $ \prod_{i = 1}^{n} b_i ?| b_i^{a_i}$ 恒成立,其中 $ b_i $ 中的每个数都是2的正整数次幂. 样例输入 3 2 3 2 3 3 3 3 2 1 10 样例输出 YES YES NO 数据范围 对于 100% 的数据有 $ n \leq 10^5,a_i \leq 10,T \leq 10$ 解析 首先拿到这道题,考场一看就知道不是规律题就是数学公式题,事实上是的. 我们可以设 $ b_i=2^{x_i}

【51NOD-0】1058 N的阶乘的长度

[算法]数学 [题解]n!的位数相当于ans=log10(n!)上取整,然后就可以拆出来加了. 可以用log10(i)或log(i)/log(10) 阶乘好像有个斯特林公式-- #include<cstdio> #include<cmath> using namespace std; int main() { int n; scanf("%d",&n); long double ans=0; for(int i=1;i<=n;i++)ans+=lo

BZOJ 4416 【SHOI2013】 阶乘字符串

题目链接:阶乘字符串 又是一道不会做的题--看了题解后我被吓傻了-- 首先我们可以有一个显然的\(O(2^nn)\)的做法.我们先预处理出\(g_{i,j}\)表示字符串中\(i\)号位置开始第一个\(j\)字符出现在什么位置.然后就可以用\(f_S\)表示使得\(S\)集合内字符的排列全都出现的最小长度,然后就可以递推了. 然后--翻了一波题解,发现当\(n>21\)的时候无解--听说合法的串长应该是\(n^2\)级别的,所以当\(n>21\)的时候就无解了--然后就可以\(O(2^nn)\