【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();
}