码迷，mamicode.com
首页 > 其他好文 > 详细

【AtCoder】AGC018

时间：2019-03-09 23:48:05      阅读：151      评论：0      收藏：0      [点我收藏+]

标签：using   [1]   ++   get   ras   mes   amp   不同   str   

A - Getting Difference

我们肯定可以得到这些数的gcd，然后判断每个数减整数倍的gcd能否得到K

#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define pdi pair<db,int>
#define mp make_pair
#define pb push_back
#define enter putchar('\n')
#define space putchar(' ')
#define eps 1e-8
#define mo 974711
#define MAXN 200005
//#define ivorysi
using namespace std;
typedef long long int64;
typedef double db;
template<class T>
void read(T &res) {
    res = 0;char c = getchar();T f = 1;
    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,K;
int a[MAXN],g;
int gcd(int a,int b) {
    return b == 0 ? a : gcd(b,a % b);
}
void Solve() {
    read(N);read(K);
    for(int i = 1 ; i <= N ; ++i) {
    read(a[i]);
    g = gcd(a[i],g);
    }
    for(int i = 1 ; i <= N ; ++i) {
    if(a[i] >= K && (a[i] - K) % g == 0) {
        puts("POSSIBLE");
        return;
    } 
    }
    puts("IMPOSSIBLE");
}
int main() {
#ifdef ivorysi
    freopen("f1.in","r",stdin);
#endif
    Solve();
    return 0;
}

B - Sports Festival

对于初始的喜爱度答案，如果我们不把含有个数最多的比赛ban掉，答案永远不会下降，所以我们试试一次次ban掉人数最多的比赛，取每次最小值

#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define pdi pair<db,int>
#define mp make_pair
#define pb push_back
#define enter putchar('\n')
#define space putchar(' ')
#define eps 1e-8
#define mo 974711
#define MAXN 200005
//#define ivorysi
using namespace std;
typedef long long int64;
typedef double db;
template<class T>
void read(T &res) {
    res = 0;char c = getchar();T f = 1;
    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,M;
int a[305][305],ans;
int cnt[305],pos[305];
bool vis[305];
void Solve() {
    read(N);read(M);
    for(int i = 1 ; i <= N ; ++i) {
    for(int j = 1 ; j <= M ; ++j) {
        read(a[i][j]);
    }
    }
    for(int i = 1 ; i <= N ; ++i) {cnt[a[i][1]]++;pos[i] = 1;}
    for(int i = 1 ; i <= M ; ++i) ans = max(ans,cnt[i]);
    for(int k = 1 ; k < M ; ++k) {
    int t = 1;
    for(int j = 1 ; j <= M ; ++j) {
        if(cnt[j] > cnt[t]) t = j;
    }
    vis[t] = 1;
    for(int i = 1 ; i <= N ; ++i) {
        cnt[a[i][pos[i]]]--;
        while(vis[a[i][pos[i]]]) {++pos[i];}
        cnt[a[i][pos[i]]]++;
    }
    int tmp = 0;
    for(int j = 1 ; j <= M ; ++j) {
        tmp = max(tmp,cnt[j]);
    }
    ans = min(ans,tmp);
    }
    out(ans);enter;
}
int main() {
#ifdef ivorysi
    freopen("f1.in","r",stdin);
#endif
    Solve();
    return 0;
}

C - Coins

把每个金币变成\((A_{i} - C_{i},B_{i} - C{i},0)\)，然后总贡献加上所有\(C[i]\)
按照\(A[i] - B[i]\)排序，然后枚举一个中间点K，在K里面选\(B[i] - C[i]\)最大的Y个，在后面选\(A[i] - C[i]\)最大的\(X\)个，取最小值

#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define pdi pair<db,int>
#define mp make_pair
#define pb push_back
#define enter putchar('\n')
#define space putchar(' ')
#define eps 1e-8
#define mo 974711
#define MAXN 200005
//#define ivorysi
using namespace std;
typedef long long int64;
typedef double db;
template<class T>
void read(T &res) {
    res = 0;char c = getchar();T f = 1;
    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 X,Y,Z;
int64 A[MAXN],B[MAXN],C[MAXN],ans,f[MAXN],b[MAXN];
int id[MAXN];
multiset<int64> S;
void Solve() {
    read(X);read(Y);read(Z);
    for(int i = 1 ; i <= X + Y + Z ; ++i) {
    read(A[i]);read(B[i]);read(C[i]);
    ans += C[i];
    A[i] -= C[i];B[i] -= C[i];
    id[i] = i;
    }
    sort(id + 1,id + X + Y + Z + 1,[](int a,int b){return A[a] - B[a] < A[b] - B[b];});
    int64 all = 0;
    for(int i = 1 ; i <= X + Y + Z ; ++i) {
    S.insert(B[id[i]]);
    all += B[id[i]];
    if(S.size() > Y) {
        all -= *(S.begin());
        S.erase(S.begin());
    }
    f[i] = all;
    }
    S.clear();all = 0;
    for(int i = X + Y + Z ; i >= 1 ; --i) {
    S.insert(A[id[i]]);
    all += A[id[i]];
    if(S.size() > X) {
        all -= *(S.begin());
        S.erase(S.begin());
    }
    b[i] = all;
    }
    int64 tmp = f[Y] + b[Y + 1];
    for(int i = Y + 1; i <= Y + Z ; ++i) {
    tmp = max(tmp,f[i] + b[i + 1]);
    }
    out(ans + tmp);enter;
}
int main() {
#ifdef ivorysi
    freopen("f1.in","r",stdin);
#endif
    Solve();
    return 0;
}

D - Tree and Hamilton Path

熟练选手直接分析每条边能否取到最大贡献
每条边最大贡献是二倍的它断开后分成子树的较小的子树大小
如果是哈密顿回路，这个贡献显然可以通过求一个重心构造出来
哈密顿路分析一下，只和起点有关，和起点逐步往子树扩大的方向走，扩大到至少大小为\(\frac{N}{2}\)时所经过的边权和，就是我们要扣除的
如果有一条边可以把子树断成\(\frac{N}{2}\)，那么减少的一定是这条边，否则其他情况都要经过这条边
否则求一个重心，找和重心相连的边最小的，其他情况一定包含一条和重心相连的边

#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define pdi pair<db,int>
#define mp make_pair
#define pb push_back
#define enter putchar('\n')
#define space putchar(' ')
#define eps 1e-8
#define mo 974711
#define MAXN 200005
//#define ivorysi
using namespace std;
typedef long long int64;
typedef double db;
template<class T>
void read(T &res) {
    res = 0;char c = getchar();T f = 1;
    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);
}
struct node {
    int to,next;int64 val;
}E[MAXN * 2];
int N,sumE,head[MAXN],siz[MAXN],son[MAXN],dep[MAXN],fa[MAXN];
int64 ans;
void add(int u,int v,int64 c) {
    E[++sumE].to = v;
    E[sumE].next = head[u];
    E[sumE].val = c;
    head[u] = sumE;
}
void dfs(int u) {
    siz[u] = 1;
    dep[u] = dep[fa[u]] + 1;
    for(int i = head[u] ; i ; i = E[i].next) {
    int v = E[i].to;
    if(v != fa[u]) {
        fa[v] = u;
        dfs(v);
        siz[u] += siz[v];
        son[u] = max(siz[v],son[u]);
        ans += min(siz[v],N - siz[v]) * 2 * E[i].val;
    }
    }
    son[u] = max(son[u],N - siz[u]);
}
void Solve() {
    read(N);
    int a,b;int64 c;
    for(int i = 1 ; i < N ; ++i) {
    read(a);read(b);read(c);
    add(a,b,c);add(b,a,c);
    }
    dfs(1);
    int t = 1;
    for(int i = 2 ; i <= N ; ++i) {
    if(son[i] < son[t] || (son[i] == son[t] && dep[i] > dep[t])) t = i;
    }
    if(N % 2 == 0 && son[t] == N / 2) {
    for(int i = head[t] ; i ; i = E[i].next) {
        if(E[i].to == fa[t]) ans -= E[i].val;
    }
    }
    else {
    int64 k = 0x7fffffff;
    for(int i = head[t] ; i ; i = E[i].next) {
        k = min(k,E[i].val);
    }
    ans -= k;
    }
    out(ans);enter;
}
int main() {
#ifdef ivorysi
    freopen("f1.in","r",stdin);
#endif
    Solve();
    return 0;
}

E - Sightseeing Plan

问网格图中某一矩形区域到顶点的和
居然是可以\(O(1)\)统计的，感觉非常长知识

具体怎么\(O(1)\)呢

问\((0,0)\)点到\((x,y)\)点的方案数，是\(\binom{x + y}{y}\)

记为\(W(x,y)\)

就是\(\sum_{i = 0}^{y}W(x,i) = W(x + 1,y)\)

就是在选\(x,i\)往右走一步，然后就一直往上走

\(\sum_{i = 0}^{x}\sum_{j = 0}^{y} W(i,j) = W(x + 1,y + 1) - 1\)

因为横竖同理，我们可以把每个都合成\(W(i + 1,y)\)，少一个\(W(0,y)\)，所以-1

这样我们求一个区域的值

\(\sum_{i = x_{1}}^{x_{2}}\sum_{j = y_{1}}^{y_{2}} W(x_{2} + 1,y_{2} + 1) - W(x_{1},y_{2} + 1) - W(x{2} + 1,y_{2}) + W(x_{1},y_{1})\)

#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define pdi pair<db,int>
#define mp make_pair
#define pb push_back
#define enter putchar('\n')
#define space putchar(' ')
#define eps 1e-8
#define mo 974711
#define MAXN 1000005
//#define ivorysi
using namespace std;
typedef long long int64;
typedef double db;
template<class T>
void read(T &res) {
    res = 0;char c = getchar();T f = 1;
    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 = 1000000007;
const int V = 2000005;
int X[7],Y[7],ans,fac[V + 5],invfac[V + 5];

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]));
}
int W(int x,int y) {
    return C(x + y,x);
}
int way1(int x,int y) {
    int s1 = x - X[2],s2 = x - X[1];
    int t1 = y - Y[2],t2 = y - Y[1];
    return inc(inc(W(s2 + 1,t2 + 1),W(s1,t1)),MOD - inc(W(s2 + 1,t1),W(s1,t2 + 1)));
}
int way2(int x,int y) {
    int s1 = X[5] - x,s2 = X[6] - x;
    int t1 = Y[5] - y,t2 = Y[6] - y;
    return inc(inc(W(s2 + 1,t2 + 1),W(s1,t1)),MOD - inc(W(s2 + 1,t1),W(s1,t2 + 1)));
}
void Solve() {
    fac[0] = 1;
    for(int i = 1 ; i <= V ; ++i) fac[i] = mul(fac[i - 1],i);
    invfac[V] = fpow(fac[V],MOD - 2);
    for(int i = V - 1 ; i >= 0 ; --i) invfac[i] = mul(invfac[i + 1],i + 1);
    for(int i = 1 ; i <= 6 ; ++i) read(X[i]);
    for(int i = 1 ; i <= 6 ; ++i) read(Y[i]);
    for(int i = X[3] ; i <= X[4] ; ++i) {
    update(ans,mul(MOD - (i + Y[3] - 1),mul(way1(i,Y[3] - 1),way2(i,Y[3]))));
    update(ans,mul(i + Y[4],mul(way1(i,Y[4]),way2(i,Y[4] + 1))));
    }
    for(int j = Y[3] ; j <= Y[4] ; ++j) {
    update(ans,mul(MOD - (X[3] - 1 + j),mul(way1(X[3] - 1,j),way2(X[3],j))));
    update(ans,mul(j + X[4],mul(way1(X[4],j),way2(X[4] + 1,j))));
    }
    out(ans);enter;
}
int main() {
#ifdef ivorysi
    freopen("f1.in","r",stdin);
#endif
    Solve();
    return 0;
}

F - Two Trees

显然如果两个点在两棵树里儿子的奇偶性不同，一定误解

否则把两棵树建出一个无向图，然后如果两个点都有偶数个儿子，那么连一条边，再建一个点，把两棵树的树根都连在上面

求一个欧拉回路，那么如果从第一棵树到第二棵树，那么这个点就是1，否则就是-1

这样的道理是什么呢。。。可以从底向上归纳，显然两棵树的叶子点有两条边，出度和入度相等，给叶子上填上相应的数，然后叶子指向父亲的相当于把+1或-1贡献给父亲，从底向上归纳也行

在学完选修2-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 MAXN 100005
#define eps 1e-12
//#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;
int A[MAXN],B[MAXN],s[MAXN][2],ans[MAXN];
struct node {
    int to,next;
}E[MAXN * 10];
int head[MAXN * 2],sumE = 1;
bool vis[MAXN * 10];
void add(int u,int v) {
    E[++sumE].to = v;
    E[sumE].next = head[u];
    head[u] = sumE;
}
void euler_road(int u) {
    for(int &i = head[u] ; i ; i = E[i].next) {
        int w = i;
        if(!vis[w]) {
            vis[w] = vis[w ^ 1] = 1;
            int v = E[w].to;
            if(u != 2 * N + 1 && v != 2 * N + 1 && abs(u - v) == N) {
                if(u < v) ans[u] = 1;
                else ans[v] = -1;
            }
            euler_road(v);
        }
    }
}
void Solve() {
    read(N);
    for(int i = 1 ; i <= N ; ++i) {
        read(A[i]);
        if(A[i] != -1) {
            add(i,A[i]);add(A[i],i);
            s[A[i]][0]++;
        }
        else {
            add(i,2 * N + 1);add(2 * N + 1,i);
        }
    }
    for(int i = 1 ; i <= N ; ++i) {
        read(B[i]);
        if(B[i] != -1) {
            add(i + N,B[i] + N);add(B[i] + N,i + N);
            s[B[i]][1]++;
        }
        else {
            add(i + N,2 * N + 1);add(2 * N + 1,i + N);
        }
    }
    for(int i = 1 ; i <= N ; ++i) {
        if((s[i][0] ^ s[i][1]) & 1) {
            puts("IMPOSSIBLE");
            return;
        }
        if(s[i][0] % 2 == 0) {
            add(i,i + N);add(i + N,i);
        }
    }
    puts("POSSIBLE");
    euler_road(1);
    for(int i = 1 ; i <= N ; ++i) {
        out(ans[i]);space;
    }
    enter;
}
int main() {
#ifdef ivorysi
    freopen("f1.in","r",stdin);
#endif
    Solve();
}

【AtCoder】AGC018

标签：using   [1]   ++   get   ras   mes   amp   不同   str   

原文地址：https://www.cnblogs.com/ivorysi/p/10503504.html
踩
(0)
赞
(0)
   
举报
评论 一句话评论（0）
登录后才能评论！
分享档案
更多>
2021年07月29日 (22)
2021年07月28日 (40)
2021年07月27日 (32)
2021年07月26日 (79)
2021年07月23日 (29)
2021年07月22日 (30)
2021年07月21日 (42)
2021年07月20日 (16)
2021年07月19日 (90)
2021年07月16日 (35)
周排行
mamicode.com排行更多图片
更多
友情链接
兰亭集智   国之画   百度统计   站长统计   阿里云   chrome插件   新版天听网
关于我们 - 联系我们 - 留言反馈
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码！