POJ 1150 The Last Non-zero Digit 数论+容斥
POJ 1150 The Last Non-zero Digit 数论+容斥
题目地址: id=1150" rel="nofollow" style="color:rgb(0,136,204);text-decoration:none;">POJ 1150
题意:
求排列P(n, m)后面第一个非0的数。
分析:
为了熟悉题目中的理论。我先做了俩0基础的题目: id=1401" rel="nofollow" style="color:rgb(0,136,204);text-decoration:none;">POJ 1401 pid=954" rel="nofollow" style="color:rgb(0,136,204);text-decoration:none;">NYOJ 954
1401 && ZOJ 2202 Factorial 阶乘N!的末尾零的个数
954 求N!二进制末尾几个0
这题想了一下。十进制末尾几个0能够转化为几个5因子。二进制最后一位非0能够转化为2因子。可是10进制就不行了,并且这个不是单单N!。而是n!/m!。orz。
实在太弱仅仅能研究题解,推荐SCAU_Lyon巨巨的题解和abilitytao巨巨的题解。
然后我坐在电脑前面看了一晚上题解,最终。我发现我太弱了,我好不easy有点看懂了。
大概讲一下吧。
我们要把排列转化成前面两题的形式。
P(n, m) = n!/(n-m)!,我们让m = n - m直接按n!/m!做,也就是求m∗(m+1)∗(m+2)∗....∗n了。
这时候,我们仅仅要统计在[m, n]这个范围里面的数的最后一位即可了,假设直接去暴力会跪。由于我们不仅要统计每一个数的最后一位,另一些2和5因子混在数中,我们还要消掉那些因子,然后取末位。
于是要把2和5抽出来,然后就仅仅剩3,7,9,统计抽完后的3,7,9,然后还得记得:由于这里面2可能会比5多,所以要把多出来的2算出来。
仅仅要知道怎样统计n!,就能统计排列。
个中计算有非常厉害的技巧,表示orz。
一、计算n!中消除2,5后末位为x的公式为:
f(n, x) = n/10 + (n%10>=x) + f(n/2, x) + f(n/2, 5) - f(n/10, x).
解释下:
1. x限定3,7,9。
2. n/10 + (n%10>=x)也就是:每十个数以内末数字是3,7,9在没有除去2和5两种因子前都仅仅会出现一次。
3. 加上抽出2和5后的子问题。(这里跟前面两题原理一样)
4. 抽出2和5的时候,会多抽出了一次10,这时候就要用容斥定理减去。
二、然后计算多余的2就比較easy理解了。就跟前面俩题一样。求出2的个数和5的个数。减一下就是n!中多出来的2的个数了。
三、然后我们就能得到[m,n]中抽出2,5后末位为3,7,9的个数,以及多出来的2的个数了。
这时候假设直接while(cnt--)去算可能会超时+爆范围。
我们能够发现2^n,3^n,7^n,9^n的末位都是有规律可寻的。于是就能直接算了。
具体细节见代码。
代码:
/*
* Author: illuz <iilluzen[at]gmail.com>
* File: 1150.cpp
* Create Date: 2014-05-26 22:28:45
* Descripton:
*/
#include <cstdio>
#include <cstring>
const int N = 2e5;
int cnt3, cnt7, cnt9, cnt2;
int n, m;
int rec[10][N];
// 计算n!中消除2,5后末位为x的数量
int f(int n, int k) {
if (n < 1)
return 0;
if (n < N && rec[k][n] != -1)
return rec[k][n];
int ret = n / 10 + (n % 10 >= k) + f(n / 2, k) + f(n / 5, k) - f(n / 10, k);
if (n < N)
rec[k][n] = ret;
return ret;
}
// 多出来的2的个数
int more(int n) {
int num2 = 0, num5 = 0;
int t = n;
while (t != 0) {
num2 += t / 2;
t /= 2;
}
while (n != 0) {
num5 += n / 5;
n /= 5;
}
return num2 - num5;
}
int main()
{
memset(rec, -1, sizeof(rec));
while (~scanf("%d%d", &n, &m)) {
if (m == 0) {
puts("1");
continue;
}
m = n - m;
cnt2 = more(n) - more(m);
cnt3 = f(n, 3) - f(m, 3);
cnt7 = f(n, 7) - f(m, 7);
cnt9 = f(n, 9) - f(m, 9);
// printf("%d %d %d %d\n", cnt2, cnt3, cnt7, cnt9);
// 2 4 8 6
if (cnt2-- == 0)
cnt2 = 1;
else if (cnt2 % 4 == 0)
cnt2 = 2;
else if (cnt2 % 4 == 1)
cnt2 = 4;
else if (cnt2 % 4 == 2)
cnt2 = 8;
else
cnt2 = 6;
// 1 3 9 7
if (cnt3 % 4 == 0)
cnt3 = 1;
else if (cnt3 % 4 == 1)
cnt3 = 3;
else if (cnt3 % 4 == 2) cnt3 = 9;
else
cnt3 = 7;
// 1 7 9 3
if (cnt7 % 4 == 0)
cnt7 = 1;
else if (cnt7 % 4 == 1)
cnt7 = 7;
else if (cnt7 % 4 == 2)
cnt7 = 9;
else
cnt7 = 3;
// 1 9
if (cnt9 % 2 == 0)
cnt9 = 1;
else
cnt9 = 9;
printf("%d\n", cnt2 * cnt3 * cnt7 * cnt9 % 10);
}
return 0;
}
POJ 1150 The Last Non-zero Digit 数论+容斥的更多相关文章
- Coprime (单色三角形+莫比乌斯反演(数论容斥))
这道题,先说一下单色三角形吧,推荐一篇noip的论文<国家集训队2003论文集许智磊> 链接:https://wenku.baidu.com/view/e87725c52cc58bd631 ...
- hdu1695:数论+容斥
题目大意: 求x属于[1,b]和 y属于[1,d]的 gcd(x,y)=k 的方案数 题解: 观察发现 gcd()=k 不好处理,想到将x=x/k,y=y/k 后 gcd(x,y)=1.. 即问题转化 ...
- CROC 2016 - Elimination Round (Rated Unofficial Edition) F - Cowslip Collections 数论 + 容斥
F - Cowslip Collections http://codeforces.com/blog/entry/43868 这个题解讲的很好... #include<bits/stdc++.h ...
- 数论 + 容斥 - HDU 4059 The Boss on Mars
The Boss on Mars Problem's Link Mean: 给定一个整数n,求1~n中所有与n互质的数的四次方的和.(1<=n<=1e8) analyse: 看似简单,倘若 ...
- 数论 + 容斥 - HDU 1695 GCD
problem's Link mean 给定五个数a,b,c,d,k,从1~a中选一个数x,1~b中选一个数y,使得gcd(x,y)=k. 求满足条件的pair(x,y)数. analyse 由于b, ...
- HDU - 2204 Eddy's爱好 (数论+容斥)
题意:求\(1 - N(1\le N \le 1e18)\)中,能表示成\(M^k(M>0,k>1)\)的数的个数 分析:正整数p可以表示成\(p = m^k = m^{r*k'}\)的形 ...
- 【bzoj2393】Cirno的完美算数教室 数论容斥
Description ~Cirno发现了一种baka数,这种数呢~只含有2和⑨两种数字~~ 现在Cirno想知道~一个区间中~~有多少个数能被baka数整除~ 但是Cirno这么天才的妖精才不屑去数 ...
- HDU 5514 Frogs (数论容斥)
题意:有n只青蛙,m个石头(围成圆圈).第i只青蛙每次只能条ai个石头,问最后所有青蛙跳过的石头的下标总和是多少? 析:首先可以知道的是第 i 只青蛙可以跳到 k * gcd(ai, m),然后我就计 ...
- UOJ #129 / BZOJ 4197 / 洛谷 P2150 - [NOI2015]寿司晚宴 (状压dp+数论+容斥)
题面传送门 题意: 你有一个集合 \(S={2,3,\dots,n}\) 你要选择两个集合 \(A\) 和 \(B\),满足: \(A \subseteq S\),\(B \subseteq S\), ...
随机推荐
- 通过实体类生成建表SQL语句实现方法
import java.io.File; import java.io.FileOutputStream; import java.lang.reflect.Field; import java.ut ...
- T2988 删除数字【状压Dp+前缀和优化】
Online Judge:从Topcoder搬过来,具体哪一题不清楚 Label:状压Dp+前缀和优化 题目描述 给定两个数A和N,形成一个长度为N+1的序列,(A,A+1,A+2,...,A+N-1 ...
- Leetcode279. Perfect Squares完全平方数
给定正整数 n,找到若干个完全平方数(比如 1, 4, 9, 16, ...)使得它们的和等于 n.你需要让组成和的完全平方数的个数最少. 示例 1: 输入: n = 12 输出: 3 解释: 12 ...
- java虚拟机(十一)--GC日志分析
GC相关:java虚拟机(六)--垃圾收集器和内存分配策略 java虚拟机(五)--垃圾回收机制GC 打印日志相关参数: -XX:+PrintGCDetails -XX:PrintGCTimestam ...
- php报错:Notice: iconv(): Wrong charset, conversion from `GBK' to `UTF8' is not allowed
代码写错了 $raw_data = iconv("GBK", "UTF8", $raw_data); 改成 $raw_data = iconv("GB ...
- layer时间插件
引入: <link rel="stylesheet" href="<{$cdnsite}>/default/common/layui/css/layui ...
- 微信小程序上传报错:以下文件没有被打包上传: · .gitignore
简单粗暴的办法就是:找到gitignore文件,把该文件删除掉即可. 在使用Git的过程中,我们喜欢有的文件比如日志,临时文件,编译的中间文件等不要提交到代码仓库,这时就要设置相应的忽略规则,来忽略这 ...
- HDU4004 二分答案
第一道二分答案...今天看了大牛的博客,突然发现有个叫“二分枚举答案”的方法好像很牛,于是便搜了些资料..发现并不是很难,可能是我了解的只是冰山一脚罢了...加油ACMer!!!! #include& ...
- DP学习之路(1) 01背包
动态规划是算法中一门很重要的思想,其通过对每一步的假设规划,不停的寻找最优最有利的解决方案,然后一步一步求解出来. 而01背包是其中最基本的一种dp思想,其题目一般为给定一个容量为V的背包,然后有n件 ...
- VS2010中Cocos2d-x中文乱码问题
不罗嗦,直接进入主题,VS2010的默认编码是"GB2312",那么以创建一个label为例,当我们使用 CCLabelTTF::create(const char *label, ...