【AtCoder】Tenka1 Programmer Contest 2019
Tenka1 Programmer Contest 2019
C - Stones
题面大意:有一个01序列,改变一个位置上的值花费1,问变成没有0在1右边的序列花费最少多少
直接枚举前i个都变成0即可
#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define mp make_pair
#define pb push_back
#define space putchar(' ')
#define enter putchar('\n')
#define eps 1e-10
#define MAXN 200005
//#define ivorysi
using namespace std;
typedef long long int64;
typedef unsigned int u32;
typedef double db;
template<class T>
void read(T &res) {
res = 0;T f = 1;char c = getchar();
while(c < '0' || c > '9') {
if(c == '-') f = -1;
c = getchar();
}
while(c >= '0' && c <= '9') {
res = res * 10 +c - '0';
c = getchar();
}
res *= f;
}
template<class T>
void out(T x) {
if(x < 0) {x = -x;putchar('-');}
if(x >= 10) {
out(x / 10);
}
putchar('0' + x % 10);
}
int N;
char s[MAXN];
int sum[MAXN];
void Solve() {
read(N);
scanf("%s",s + 1);
for(int i = 1 ; i <= N ; ++i) {
sum[i] = sum[i - 1];
if(s[i] == '#') sum[i]++;
}
int ans = N;
for(int i = 0 ; i <= N ; ++i) {
ans = min(ans,sum[i] + N - i - (sum[N] - sum[i]));
}
out(ans);enter;
}
int main() {
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Solve();
}
D - Three Colors
大意:将数分到三个标号为123的集合,要求三个集合能构成三角形
题解:也就是\(R + B + G = S\)满足\(G,R,B < \frac{S}{2}\)
不合法的情况就是有一个集合大于等于\(\frac{S}{2}\)的方案,对于等于\(\frac{S}{2}\)再用一个dp去重即可
#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define mp make_pair
#define pb push_back
#define space putchar(' ')
#define enter putchar('\n')
#define eps 1e-10
#define MAXN 200005
//#define ivorysi
using namespace std;
typedef long long int64;
typedef unsigned int u32;
typedef double db;
template<class T>
void read(T &res) {
res = 0;T f = 1;char c = getchar();
while(c < '0' || c > '9') {
if(c == '-') f = -1;
c = getchar();
}
while(c >= '0' && c <= '9') {
res = res * 10 +c - '0';
c = getchar();
}
res *= f;
}
template<class T>
void out(T x) {
if(x < 0) {x = -x;putchar('-');}
if(x >= 10) {
out(x / 10);
}
putchar('0' + x % 10);
}
const int MOD = 998244353;
int N;
int a[305],S;
int dp[305][90005],f[90005];
int inc(int a,int b) {
return a + b >= MOD ? a + b - MOD : a + b;
}
int mul(int a,int b) {
return 1LL * a * b % MOD;
}
void update(int &x,int y) {
x = inc(x,y);
}
int fpow(int x,int c) {
int res = 1,t = x;
while(c) {
if(c & 1) res = mul(res,t);
t = mul(t,t);
c >>= 1;
}
return res;
}
void Solve() {
read(N);
for(int i = 1 ; i <= N ; ++i) {
read(a[i]);
S += a[i];
}
dp[0][0] = 1;
f[0] = 1;
for(int i = 1 ; i <= N ; ++i) {
for(int j = S ; j >= 0 ; --j) {
update(dp[i][j],mul(2,dp[i - 1][j]));
if(j >= a[i]) {
update(dp[i][j],dp[i - 1][j - a[i]]);
update(f[j],f[j - a[i]]);
}
}
}
int ans = fpow(3,N);
for(int j = 0 ; j <= S ; ++j) {
if(j * 2 >= S) {
update(ans,MOD - mul(3,dp[N][j]));
}
}
if(S % 2 == 0) {
update(ans,mul(f[S / 2],3));
}
out(ans);enter;
}
int main() {
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Solve();
}
E - Polynomial Divisors
大意:给定多项式求多项式取所有整数值时都能被整除的质数
盲猜肯定要么是所有系数的gcd的质因数要么就小于最高次项,然后判解不会就直接暴力,以为数据造不满然后被卡死了,atc的数据还是强
事实上判解的话只要判这个多项式在modp意义下是否是多项式\(x^{p} - x\)的倍数即可
充分性显然(费马小定理)
必要性合法的式子会含有\(x(x - 1)(x - 2)..(x - (p - 1))\),然后这个东西在modp意义下是\(x^{p} - x\)相等,如果不是全等的话,必然有一个小于等于p - 1的多项式是有p个根,这不存在
#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define mp make_pair
#define pb push_back
#define space putchar(' ')
#define enter putchar('\n')
#define eps 1e-10
#define MAXN 10005
//#define ivorysi
using namespace std;
typedef long long int64;
typedef unsigned int u32;
typedef double db;
template<class T>
void read(T &res) {
res = 0;T f = 1;char c = getchar();
while(c < '0' || c > '9') {
if(c == '-') f = -1;
c = getchar();
}
while(c >= '0' && c <= '9') {
res = res * 10 +c - '0';
c = getchar();
}
res *= f;
}
template<class T>
void out(T x) {
if(x < 0) {x = -x;putchar('-');}
if(x >= 10) {
out(x / 10);
}
putchar('0' + x % 10);
}
set<int> S;
int g;
int N,a[MAXN],b[MAXN],c[MAXN];
int prime[MAXN],tot;
bool nonprime[MAXN];
int gcd(int a,int b) {
return b == 0 ? a : gcd(b,a % b);
}
void Solve() {
read(N);
for(int i = N ; i >= 0 ; --i) {
read(a[i]);
g = gcd(g,abs(a[i]));
}
for(int i = 2 ; i <= g / i ; ++i) {
if(g % i == 0) {
S.insert(i);
while(g % i == 0) g /= i;
}
}
if(g > 1) S.insert(g);
for(int i = 2 ; i <= 10000 ; ++i) {
if(!nonprime[i]) {
prime[++tot] = i;
}
for(int j = 1 ; j <= tot ; ++j) {
if(1LL * i * prime[j] > 10000) break;
nonprime[i * prime[j]] = 1;
if(i % prime[j] == 0) break;
}
}
for(int i = 1 ; i <= tot ; ++i) {
if(prime[i] > N) break;
int p = prime[i];
if(S.find(p) != S.end()) continue;
memset(b,0,sizeof(b));
memset(c,0,sizeof(c));
bool flag = 1;
for(int i = N ; i >= 0 ; --i) {
b[i] = (a[i] + c[i]) % p;
if(i >= p) {
c[1 + i - p] += b[i];
c[1 + i - p] %= p;
}
else if(b[i]){
flag = 0;break;
}
}
if(flag) S.insert(p);
}
for(auto t : S) {
out(t);enter;
}
}
int main() {
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Solve();
}
F - Banned X
大意:找到一个长度为N序列只含0,1,2,求其中任意一段连续的和不为x的答案
首先统计所有和小于x的答案
然后就是如果x为奇数,我只选2且序列总和大于x的答案
然后就是,因为我每一段的和都不是x,每次增加是1或2,那么肯定存在一个和是\(X - 1\)的前缀和
此时我只能选2,如果序列开头第一个不为0的数是2,那么我可以删掉这个2,往后继续加一个2,直到我删掉时遇到了一个1
所以我们枚举i个2+1个1的序列,枚举构成X - 1所用的序列长度,之后只能往后加不超过i个2即可
#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define mp make_pair
#define pb push_back
#define space putchar(' ')
#define enter putchar('\n')
#define eps 1e-10
#define MAXN 10005
//#define ivorysi
using namespace std;
typedef long long int64;
typedef unsigned int u32;
typedef double db;
template<class T>
void read(T &res) {
res = 0;T f = 1;char c = getchar();
while(c < '0' || c > '9') {
if(c == '-') f = -1;
c = getchar();
}
while(c >= '0' && c <= '9') {
res = res * 10 +c - '0';
c = getchar();
}
res *= f;
}
template<class T>
void out(T x) {
if(x < 0) {x = -x;putchar('-');}
if(x >= 10) {
out(x / 10);
}
putchar('0' + x % 10);
}
const int MOD = 998244353;
int N,X;
int dp[3005][6005],fac[6005],invfac[6005];
int f[3005];
int inc(int a,int b) {
return a + b >= MOD ? a + b - MOD : a + b;
}
int mul(int a,int b) {
return 1LL * a * b % MOD;
}
void update(int &x,int y) {
x = inc(x,y);
}
int fpow(int x,int c) {
int res = 1,t = x;
while(c) {
if(c & 1) res = mul(res,t);
t = mul(t,t);
c >>= 1;
}
return res;
}
int C(int n,int m) {
if(n < m) return 0;
return mul(fac[n],mul(invfac[m],invfac[n - m]));
}
void Solve() {
read(N);read(X);
dp[0][0] = 1;
fac[0] = 1;
for(int i = 1 ; i <= 2 * N ; ++i) fac[i] = mul(fac[i - 1],i);
invfac[2 * N] = fpow(fac[2 * N],MOD - 2);
for(int i = 2 * N - 1 ; i >= 0 ; --i) {
invfac[i] = mul(invfac[i + 1],i + 1);
}
int ans = 0;
for(int i = 1 ; i <= N ; ++i) {
for(int j = 2 * i ; j >= 1 ; --j) {
if(j >= 2) update(dp[i][j],dp[i - 1][j - 2]);
if(j >= 1) update(dp[i][j],dp[i - 1][j - 1]);
}
}
for(int i = 0 ; i <= N ; ++i) {
for(int j = 0 ; j <= 2 * i ; ++j) {
if(j < X) update(ans,mul(dp[i][j],C(N,i)));
}
}
for(int i = 1 ; i < N ; ++i) {
int t = i * 2 + 1;
if(t >= X) break;
for(int j = 0 ; j <= N ; ++j) {
if(i + j + 2 <= N) {
update(f[i + j + 2],dp[j][X - 1 - t]);
if(i + j + i + 2 <= N) update(f[i + j + i + 2],MOD - dp[j][X - 1 - t]);
}
}
}
for(int i = 1 ; i <= N ; ++i) f[i] = inc(f[i],f[i - 1]);
for(int i = 1 ; i <= N ; ++i) {
update(ans,mul(f[i],C(N,i)));
}
if(X & 1) {
for(int i = 1 ; i <= N ; ++i) {
if(i * 2 > X) update(ans,C(N,i));
}
}
out(ans);enter;
}
int main() {
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Solve();
}
【AtCoder】Tenka1 Programmer Contest 2019的更多相关文章
- 【AtCoder】Tenka1 Programmer Contest
C - 4/N 列出个方程枚举解一下 #include <bits/stdc++.h> #define fi first #define se second #define pii pai ...
- 【AtCoder】Tenka1 Programmer Contest(C - F)
C - Align 考的时候,我大胆猜了结论,就是一小一大一小一大这么排 证明的话,由于我们总是要加上相邻的最大值而减去最小值,我们就让最大值都保持在前面 如果长度为奇数,要么就是大小大小大,要么是小 ...
- 【AtCoder】AISing Programming Contest 2019
本来以为是1199rated的..仔细一看发现是1999,所以就做了一下 这场涨分很轻松啊...为啥又没打 等pkuwc考完我一定打一场atcoder(咕咕咕,咕咕咕,咕咕咕咕咕咕咕~) 但是其实我思 ...
- 【AtCoder】Yahoo Programming Contest 2019
A - Anti-Adjacency K <= (N + 1) / 2 #include <bits/stdc++.h> #define fi first #define se se ...
- 【AtCoder】KEYENCE Programming Contest 2019
A - Beginning 这个年份恐怕需要+2 #include <bits/stdc++.h> #define fi first #define se second #define p ...
- Atcoder Tenka1 Programmer Contest 2019 题解
link 题面真简洁 qaq C Stones 最终一定是连续一段 . 加上连续一段 # .直接枚举断点记录前缀和统计即可. #include<bits/stdc++.h> #define ...
- Atcoder Tenka1 Programmer Contest 2019
C 签到题,f[i][0/1]表示以i结尾最后一个为白/黑的最小值,转移显然. #include<bits/stdc++.h> using namespace std; ; ]; char ...
- Atcoder Tenka1 Programmer Contest 2019 E - Polynomial Divisors
题意: 给出一个多项式,问有多少个质数\(p\)使得\(p\;|\;f(x)\),不管\(x\)取何值 思路: 首先所有系数的\(gcd\)的质因子都是可以的. 再考虑一个结论,如果在\(\bmod ...
- Atcoder Tenka1 Programmer Contest 2019 D Three Colors
题意: 有\(n\)个石头,每个石头有权值,可以给它们染'R', 'G', 'B'三种颜色,如下定义一种染色方案为合法方案: 所有石头都染上了一种颜色 令\(R, G, B\)为染了'R', 染了'G ...
随机推荐
- 【转】C++对成员访问运算符->的重载
运算符->的重载比较特别,它只能是非静态的成员函数形式,而且没有参数. 1.如果返回值是一个原始指针,那么就将运算符的右操作数当作这个原始指针所指向类型的成员进行访问: 2.如果返回值是另一个类 ...
- v-on 绑定自定义事件
每个 Vue 实例都实现了事件接口,即: 使用 $on(eventName) 监听事件 使用 $emit(eventName) 触发事件 Vue 的事件系统与浏览器的 EventTarget API ...
- web@css高级选择器(after,befor用法),基本css样式
1.高阶选择器:子代后代,相邻通用兄弟,交集并集,属性,伪类,伪元素子代后代选择器 div>p{} div p{}相邻通用兄弟 div+p{} div~p{}理解:div同学的同桌p di ...
- centos7怎么永久修改hosname
centos7怎么永久修改hosname 其实,一般来说安装好虚拟机之后,一般都会进行修改hostname,之前也是在修改的时候,遇到过问题,但是没有深究,今天在修改的时候,好好研究了一下,之前看到好 ...
- C/C++中容易造成内存溢出的函数
1.strcpy() strcpy()函数将源字符串复制到缓冲区.没有指定要复制字符的具体数目.复制字符的数目直接取决于源字符串中的数目.如果源字符串碰巧来自用户输入,且没有专门限制其大小,则有可能会 ...
- php如何实现图片点击下载,并保存本地?-----本例子为二维码的生成图片,并支持点击下载
### 今天因为工作需要,完成了一个二维码的生成图片,并支持点击下载的 ### 控制器文件,相关代码 // 生成二维码 $url = action('Apih5\\VersionController@ ...
- 【原创】Linux基础之命令行浏览器links
有时服务器环境受限,比如在内网环境不能暴露端口从外网访问,用curl看html代码比较累,这时可以使用命令行浏览器来查看相关页面 links 官方:http://links.twibright.com ...
- ELK Packetbeat 部署指南
http://www.ttlsa.com/elk/elk-packetbeat-deployment-guide/
- Bootstrap 固定底部导航栏菜单
直接上代码: <!DOCTYPE html> <html> <head> <meta http-equiv="content-type" ...
- 对象Object
功能分类 1. 创建对象 把各对数自身拥有的可枚举属性复制到第一个对象并返回:obj = Object.assign(o1, o2, o3),o1=obj ...