紫书 例题10-6 无关的元素 UVa1635
来源:互联网 发布:淘宝商城二手市场 编辑:程序博客网 时间:2024/05/21 22:29
题意:给定n个数a1,a2····an,依次求出相邻两个数值和,将得到一个新数列,重复上述操作,最后结果将变为一个数,问这个数除以m的余数与那些数无关?例如n=3,m=2时,第一次得到a1+a2,a2+a3,在求和得a1+2*a2+a3,它除以2的余数和a2无关。1<=n<=10^5, 2<=m<=10^9
分析:
通过试验可以发现,本题等价于求解C(n-1,i)的组合数中有哪些是m的倍数,可以利用唯一分解定理来判断:事先分解m,随后利用递推式计算每一项中包含m的素因数的指数即可。
代码如下:
////Created by BLUEBUFF 2016/1/11//Copyright (c) 2016 BLUEBUFF.All Rights Reserved//#pragma comment(linker,"/STACK:102400000,102400000")//#include <ext/pb_ds/assoc_container.hpp>//#include <ext/pb_ds/tree_policy.hpp>//#include <ext/pb_ds/hash_policy.hpp>//#include <bits/stdc++.h>#include <set>#include <map>#include <queue>#include <stack>#include <cmath>#include <cstdio>#include <time.h>#include <cstdlib>#include <cstring>#include <complex>#include <sstream> //isstringstream#include <iostream>#include <algorithm>using namespace std;//using namespace __gnu_pbds;typedef long long LL;typedef unsigned long long uLL;typedef pair<int, LL> pp;#define REP1(i, a, b) for(int i = a; i < b; i++)#define REP2(i, a, b) for(int i = a; i <= b; i++)#define REP3(i, a, b) for(int i = a; i >= b; i--)#define CLR(a, b) memset(a, b, sizeof(a))#define MP(x, y) make_pair(x,y)template <class T1, class T2>inline void getmax(T1 &a, T2 b) { if (b>a)a = b; }template <class T1, class T2>inline void getmin(T1 &a, T2 b) { if (b<a)a = b; }const int maxn = 100010;const int maxm = 1e5+5;const int maxs = 10;const int maxp = 1e3 + 10;const int INF = 1e9;const int UNF = -1e9;const int mod = 1e9 + 7;const int rev = (mod + 1) >> 1; // FWT//const double PI = acos(-1);//headint fac1[100][2]; //fac[i][0]存放素因数, fac[i][1]存放其指数int fac2[100], a[maxn];void factor(int m){ int &num = fac1[0][0];//fac[0][0]是表头,存放总的个数,用引用比较方便 num = 0; for(int i = 2; i * i <= m; i++){ if(m % i == 0){ fac1[++num][0] = i; fac1[num][1] = 0; do{ fac1[num][1]++; m /= i; }while(m % i == 0); //将i除干净 } } if(m > 1)//如果分解到最后m仍然大于1,说明它是一个素数。注意:如果只是判断素因子有哪些,可以没有此处判断,否则必须有此步 { fac1[++num][0] = m; fac1[num][1] = 1; }}bool check(int n, int j)//按照递推公式来计算第j项,检查唯一分解式的指数{ int num = fac1[0][0]; int a = n - j, b = j; for(int i = 1; i <= num; i++){ int p = fac1[i][0]; int &q = fac2[i]; for(; a % p == 0; a /= p, q++);//为了提高效率,只用检验m的分解式中的素因数即可 for(; b % p == 0; b /= p, q--); } for(int i = 1; i <= num; i++){ if(fac1[i][1] > fac2[i]){ return false; } } return true;}int main(){ int n, m; while(scanf("%d%d", &n, &m) != EOF) { int cnt = 0; factor(m); memset(fac2, 0, sizeof(fac2)); for(int i = 1; i < n; i++){ if(check(n, i)){ a[cnt++] = i + 1; } } printf("%d\n", cnt); for(int i = 0; i < cnt; i++){ printf("%d%c", a[i], i == (cnt - 1) ? '\n' : ' '); } } return 0;}
0 0
- 例题10-6 无关的元素 UVa1635
- 紫书 例题10-6 无关的元素 UVa1635
- 【原创】【组合数学】Uva1635 无关的元素
- UVa1635
- 例10-6 uva1635(唯一分解定理)
- 无关元素
- 【学校OJ】组合数学 无关的元素
- 无关的元素(排列与组合) By ACReaper
- 2017-10-31C#引申的一个例题——不重复抽取元素的组合
- 紫书 例题 10-13 危险的组合
- 元素的比较次数与元素初始排序无关的排序方法
- 算法入门经典第二版第10章uva1635
- iOS 10打印 无关的日志
- 证明你真正富有的8个元素,其中7个与金钱无关
- 算法竞赛入门经典:第十章 数学概念与方法 10.7无关的元素
- 唯一分解式,二项式定理(无关的元素,uva 1635)
- 紫书 例题 10-2 不爽的裁判 UVa12169 ex_gcd
- 紫书 例题10-25 约瑟夫的数论问题 UVa1363
- Spring MVC集成Tiles2
- 多态原理探究_证明vptr指针的存在
- ubuntu服务器双ip(一条公网一条内网 )配置:
- Python mysql实践小记
- View设置了圆角背景长按出现默认边角
- 紫书 例题10-6 无关的元素 UVa1635
- 在 XML 中有 5 个预定义的实体引用
- ES 深入搜索04_Proximity Matching
- bootstrap 单选/开关控件
- css-包含块框外边距叠加
- tomcat 不能远程访问
- 哈夫曼树以及哈夫曼编码
- iOS 获取视频的任意一帧
- ES 深入搜索05_Partial Matching