Atcoder Tenka1 Programmer Contest 2019 E - Polynomial Divisors
题意:
给出一个多项式,问有多少个质数\(p\)使得\(p\;|\;f(x)\),不管\(x\)取何值
思路:
首先所有系数的\(gcd\)的质因子都是可以的。
再考虑一个结论,如果在\(\bmod p\)意义下,多项式中存在\((x^p - x)\)这个因式,那么这个质数\(p\)也是可以的
显然\(p \leq n\),那么我们只要枚举每个\(\leq n\)的质数,做模\(p\)意义下的多项式除法,判断余数是否为\(0\)即可。
证明:
- 充分性:考虑\(p\;|\;f(x)\),即\(f(x) = kp\),即在\(\bmod p\)意义下,\(f(x) = 0\),根据欧拉定理,分两种情况讨论
- \(x < p\),又因为\(p\)是质数,那么显然有\((x, p) = 1\),那么\(x^{p - 1} \equiv 1 \pmod p\),有\(x^{p} - x \equiv 0 \pmod p\)
- \(x \geq p\),如果\(gcd(x, p)\)不为\(1\),那么显然有\(gcd(x, p) = p\),那么已经满足\(p\;|\;f(x)\),否则套用欧拉定理
- 必要性:如果\(p\;|\;f(x)\),那么\(0, 1, \cdots, p - 1\)必然为\(f(x)\)的一个根,那么\(f(x)\)有因式\(x(x - 1)(x - 2)\cdots(x - (p - 1))\)。我们考虑这个因式与\(x^p - x\)是等价的,如果不是等价的,那么作差之后,最高次变为\(p - 1\),而根的个数却有\(p\)个,显然矛盾
#include <bits/stdc++.h>
using namespace std;
#define N 10010
int n, a[N], b[N];
bool isprime(int x) {
for (int i = 2; 1ll * i * i <= x; ++i) {
if (x % i == 0) {
return false;
}
}
return true;
}
int gcd(int a, int b) {
return b ? gcd(b, a % b) : a;
}
bool ok(int p) {
if (a[0] % p) {
return false;
}
for (int i = 0; i <= n; ++i) {
b[i] = a[i];
}
for (int i = n; i >= p - 1; --i) {
(b[i - (p - 1)] += b[i]) %= p;
b[i] = 0;
}
for (int i = 0; i <= n; ++i) {
if (b[i] % p) {
return false;
}
}
return true;
}
int main() {
while (scanf("%d", &n) != EOF) {
int G = 0;
for (int i = 0; i <= n; ++i) {
scanf("%d", a + i);
G = gcd(G, abs(a[i]));
}
reverse(a, a + 1 + n);
vector <int> res;
for (int i = 2; 1ll * i * i <= G; ++i) {
if (G % i == 0) {
res.push_back(i);
while (G % i == 0) {
G /= i;
}
}
}
if (G > 1) {
res.push_back(G);
}
for (int i = 2; i <= n; ++i) {
if (isprime(i) && ok(i)) {
res.push_back(i);
}
}
sort(res.begin(), res.end());
res.erase(unique(res.begin(), res.end()), res.end());
for (auto it : res) {
printf("%d\n", it);
}
// puts("------------");
}
return 0;
}
Atcoder Tenka1 Programmer Contest 2019 E - Polynomial Divisors的更多相关文章
- Atcoder Tenka1 Programmer Contest 2019 题解
link 题面真简洁 qaq C Stones 最终一定是连续一段 . 加上连续一段 # .直接枚举断点记录前缀和统计即可. #include<bits/stdc++.h> #define ...
- Atcoder Tenka1 Programmer Contest 2019题解
传送门 \(C\ Stones\) 最后肯定形如左边一段白+右边一段黑,枚举一下中间的断点,预处理一下前缀和就可以了 int main(){ // freopen("testdata.in& ...
- Atcoder Tenka1 Programmer Contest 2019
C 签到题,f[i][0/1]表示以i结尾最后一个为白/黑的最小值,转移显然. #include<bits/stdc++.h> using namespace std; ; ]; char ...
- Atcoder Tenka1 Programmer Contest 2019 D Three Colors
题意: 有\(n\)个石头,每个石头有权值,可以给它们染'R', 'G', 'B'三种颜色,如下定义一种染色方案为合法方案: 所有石头都染上了一种颜色 令\(R, G, B\)为染了'R', 染了'G ...
- 【AtCoder】Tenka1 Programmer Contest 2019
Tenka1 Programmer Contest 2019 C - Stones 题面大意:有一个01序列,改变一个位置上的值花费1,问变成没有0在1右边的序列花费最少多少 直接枚举前i个都变成0即 ...
- Atcoder Tenka1 Programmer Contest D: IntegerotS 【思维题,位运算】
http://tenka1-2017.contest.atcoder.jp/tasks/tenka1_2017_d 给定N,K和A1...AN,B1...BN,选取若干个Ai使它们的或运算值小于等于K ...
- Atcoder Tenka1 Programmer Contest C C - 4/N
http://tenka1-2017.contest.atcoder.jp/tasks/tenka1_2017_c 我怀疑我是不是智障.... 本来一直的想法是能不能构造出答案,把N按奇偶分,偶数好办 ...
- Tenka1 Programmer Contest 2019
C:即要使前一部分为白色后一部分为黑色,枚举分割点前缀和计算答案取min即可. #include<bits/stdc++.h> using namespace std; #define l ...
- Tenka1 Programmer Contest 2019 D - Three Colors
Three Colors 思路:dp 设sum为所有边的总和 不能组成三角形的情况:某条边长度>=ceil(sum/2),可以用dp求出这种情况的方案数,然后用总方案数减去就可以求出答案. 注意 ...
随机推荐
- Junit单元测试初识
写过单元测试的小童鞋对于Junit一定不陌生,可小白我,刚刚开始接触,这里就把我的测试实验,做一下记录,以便以后方便查看.学习使用JUnit4,既然使用最新版本了,就不要再考虑老版本是如何使用的了,J ...
- Centos 7.x系统安装后的初始化配置
1.配置IP.网关,编辑/etc/sysconfig/network-scripts/ifcfg-eno16777736 TYPE=Ethernet BOOTPROTO=none //默认为dhcp ...
- 【C#】浏览器源代码
using System; using System.Collections.Generic; using System.ComponentModel; using System.Data; usin ...
- java读取写入oracle的blob字段工具类
import com.hzunitech.fxgk.sys.model.UtFileData;import com.jfinal.kit.PathKit;import com.jfinal.plugi ...
- iOS - 音乐播放器之怎么获取音乐列表
方法一: 这个方法是通过获取到沙盒路径,来得到音乐的路径(使用这个方法需要把音乐放进沙盒) NSFileManager *manager = [NSFileManager defaultManager ...
- Shell if else
语句1if [ expression ]then ...fi 语句2if [ expression ]then ...else ...fi 语句3if [ expression 1 ]then ... ...
- 找到多个与名为“xxx”的控制器匹配的类型。如果为此请求(“{controller}/{action}/{id}”)提供服务的路由没有指定命名空间以搜索与此请求相匹配的控制器,则会发生这种情况。
一次在建MVC 项目的进行开发的时候,因为后来想到了一个更好的项目名称,就把 Web项目的名称重命名 改了, 然后 程序集名称,默认命名空间,都改成新的了,刚建立的项目本身也不大,运行起来,总是报 & ...
- HDU 3038 - How Many Answers Are Wrong - [经典带权并查集]
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3038 Time Limit: 2000/1000 MS (Java/Others) Memory Li ...
- POJ 3067 - Japan - [归并排序/树状数组(BIT)求逆序对]
Time Limit: 1000MS Memory Limit: 65536K Description Japan plans to welcome the ACM ICPC World Finals ...
- Crontab '2>&1 &' 含义
Crontab '2>&1 &' 含义 - dosttyy - 博客园 https://www.cnblogs.com/dosttyy/p/4810026.html