UVa1635 - Irrelevant Elements
通过观察发现其规律符合杨辉三角
需要注意的是最后ai的系数是C(i-1,n-1)
那么,问题就可以变成判断C(0,n-1),C(1,n-1)。。。。C(n-1,n-1)哪些是m的倍数
只需要计算出m的唯一分解式中各个素因子在C(i-1,n-1)中的指数即可完成判断
然而为了节省时间,实际上我们只需算出m的每一个素因子在C(i-1,n-1)项中 含有几个即可
即我们将c(i-1,n-1)依次除以m的每一个素因子,直到无法整出,即可得出该项素因子的个数
紫薯上给出一个公式C(k,n)=(n-k+1)*c(k-1,n)/k
可以用这个公式递推出所要求出的素因子的个数
具体怎么做呢?
说明:f[i].first中存的时m的素因子,f[i].second中存的是其对应的个数
d[i]用来存c(i-1,n-1)中素因子的个数
int d[50]={0};
for (int k=;k<=n;k++){
int x=n-k+;
for (int i=;i<f.size();i++){
int p=f[i].first;
d[i]+=x/p;
x%=p;
}
x=k-;
for (int i=;i<f.size();i++){
int p=f[i].first;
d[i]-=x/p;
x%=p;
}
int flag=;
for (int i=;i<f.size();i++){
if (d[i]<f[i].second) {flag=;break;}
}
if (flag) cout<<k<<" ok"<<endl;else cout<<"no"<<endl;
}
注意精度问题,所以采用分开除的方法。例外注意是利用的d[i]这个数组的生存周期是在总的循环里实现的递推
这里给出网上某大神的优秀代码
转:http://blog.csdn.net/u014800748/article/details/43927205
#define _CRT_SECURE_NO_WARNINGS
#include<iostream>
#include<algorithm>
#include<string>
#include<sstream>
#include<set>
#include<vector>
#include<stack>
#include<map>
#include<queue>
#include<deque>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<ctime>
#include<functional>
using namespace std; #define N 100005
int fac[][];//一张表,fac[i][0]存放素因数,fac[i][1]存放其指数
int fac_c[];
int a[N]; void factor(int m)//分解m
{
int&num = fac[][];//fac[0][0]是表头,存放总的个数,用引用比较方便
num = ;
for (int i = ; i*i <= m;i++)
if (m%i == )
{
fac[++num][] = i;
fac[num][] = ;
do
{
fac[num][]++;
m /= i;
} while (m%i == );//将i除干净
}
if (m > )//如果分解到最后m仍然大于1,说明它是一个素数。注意:如果只是判断素因子有哪些,可以没有此处判断,否则必须有此步
{
fac[++num][] = m;
fac[num][] = ;
}
} bool check(int n, int j)//按照递推公式来计算第j项,检查唯一分解式的指数
{
int num = fac[][];
int a = n - j;//其实是((n-1)+j+1)化简后的结果
int b = j;
for (int i = ; i <= num; i++)
{
int p = fac[i][];
int&q = fac_c[i];
for (; a%p == ; a /= p, q++);//为了提高效率,只用检验m的分解式中的素因数即可
for (; b%p == ; b /= p, q--);
}
for (int i = ; i <= num;i++)
if (fac[i][] > fac_c[i])
return false;
return true;
} int main()
{
//freopen("test.txt", "r", stdin);
int n, m;
while (cin >> n >> m)
{
int cnt = ;
factor(m);
memset(fac_c, , sizeof(fac_c));
for (int i = ; i < n;i++)//直接检查1到n-1项(从0开始)
if (check(n, i))
a[cnt++] = i + ;
printf("%d\n", cnt);
for (int i = ; i < cnt; i++)
printf("%s%d", i == ? "" : " ", a[i]);
printf("\n");
}
return ;
}
UVa1635 - Irrelevant Elements的更多相关文章
- UVA1635 Irrelevant Elements —— 唯一分解定理 + 二项式定理
题目链接:https://vjudge.net/problem/UVA-1635 (紫书320) 题解: 1.根据二项式定理, 可得递推公式: C(n,k) = (n-k+1)/k * C(n, k- ...
- UVA1635 Irrelevant Elements(唯一分解定理 + 组合数递推)
http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=51196 紫书P320; 题意:给定n个数a1,a2····an,依次求出相邻 ...
- Irrelevant Elements UVA - 1635 二项式定理+组合数公式+素数筛+唯一分解定理
/** 题目:Irrelevant Elements UVA - 1635 链接:https://vjudge.net/problem/UVA-1635 题意:給定n,m;題意抽象成(a+b)^(n- ...
- POJ2167 Irrelevant Elements
Time Limit: 5000MS Memory Limit: 65536KB 64bit IO Format: %lld & %llu Description Young cryp ...
- POJ 2167 Irrelevant Elements 质因数分解
Irrelevant Elements Time Limit: 5000MS Memory Limit: 65536K Total Submissions: 2231 Accepted: 55 ...
- 【UVa1635】Irrelevant Elements - 唯一分解定理
题意 给你 \(n\) 个数,每次求出相邻两个数的和组成新数列.经过 \(n-1\) 次操作后,得到一个数.求这个数 \(mod \ m\) 与哪些项无关. 如:当 \(m=2 \ , \ n=2\) ...
- Irrelevant Elements, ACM/ICPC NEERC 2004, UVa1635
这种题目最重要的是思路了清晰 #include <cstdio> #include <cstring> ;//sqrt(n)+1 is enough ][]; ]; int a ...
- Irrelevant Elements UVA-1635 (二项式定理)
vjudge链接 原题链接 乍一看似乎没什么思路,但是写几个简单的例子之后规律就变得很明显. 比如当 n=5 时,每一步计算后的结果如下: a1 a1+a2 a1+2a2+a3 a1+3a2+3a3+ ...
- UVA 1635 Irrelevant Elements
https://vjudge.net/problem/UVA-1635 题意:n个数,每相邻两个求和,最后变成1个数,问这个数除m的余数与第几个数无关 n个数使用次数分别为C(n-1,i) i∈[0, ...
随机推荐
- 201521123096《Java程序设计》第八周学习总结
1. 本周学习总结 1.1 以你喜欢的方式(思维导图或其他)归纳总结集合与泛型相关内容. 2. 书面作业 本次作业题集集合 1.List中指定元素的删除(题目4-1) 1.1 实验总结 实验中使用了s ...
- 201521123031 《Java程序设计》第4周学习总结
---恢复内容开始--- 1. 本周学习总结 1.1 尝试使用思维导图总结有关继承的知识点. 1.2 使用常规方法总结其他上课内容. (1)父类只能有一个,即单继承,子类继承父类的全部成员(属性和方法 ...
- 201521123040《Java程序设计》第2周学习总结
1. 本周学习总结 本周学习目录: <1>回顾(Java程序) <2>掌握各种数据类型的使用 <3>掌握运算符和表达式的使用 <4>枚举类.String ...
- 201521123089 《Java程序设计》第1周学习总结
1. 本周学习总结 学习了java的发展史以及JDK,JVM,JRE; 2. 书面作业 1.为什么java程序可以跨平台运行?执行java程序的步骤是什么? java有JVM,实在虚拟机上运行的. 2 ...
- 201521123027 <java程序设计>第九周学习总结
1.本周学习总结 1.1 以你喜欢的方式(思维导图或其他)归纳总结异常相关内容. 2.书面作业 Q1.常用异常 题目5-1 1.1 截图你的提交结果(出现学号) 1.2 自己以前编写的代码中经常出现什 ...
- MarkDown 例子
一个例子: 例子开始 1. 本章学习总结 今天主要学习了三个知识点 封装 继承 多态 2. 书面作业 Q1. java HelloWorld命令中,HelloWorld这个参数是什么含义? 今天学了一 ...
- 为bookStore添加权限【动态代理和注解】
前言 目前为止,我们已经学习了动态代理技术和注解技术了.于是我们想要为之前的bookStore项目添加权限控制-.. 只有用户有权限的时候,后台管理才可以进行相对应的操作-.. 实现思路 之前我们做权 ...
- 易语言 【寻找文本】命令的bug
最近在重写易语言模块的时候,在取子文本操作时老是出错,经常出现一些奇怪的问题,一开始以为是代码问题,可是找半天硬是找不到问题所在. 于是进入了找bug模式,这么几行代码,看了我半个小时,左改右改,总感 ...
- OpenShift中的持续交付
上一文中讲述了如何在AWS下搭建OpenShift集群.这篇文章将目光转向如何在OpenShift中实现CI/CD以及产品环境的部署. 持续交付 如果要打造一个持续交付的流水线,首先要考虑多环境的问题 ...
- 记一次 node.js 的populate用法
最近在学习node.js做一个简单的博客的网站,用的express框架和mongodb的数据库.之前没有接触过这个数据库,所有在一开始写的时候遇到了一些问题,如何初始化model类型,又如何实现简单的 ...