【AtCoder】AGC028 (A-E)题解
A - Two Abbreviations
如果有最小答案的话这个答案一定是N和M的lcm
我们考虑一下什么情况下
\(k \frac{L}{N} = h \frac{L}{M}\)且\(k,g\)互质
显然是在\(k = \frac{N}{gcd(N,M)},h = \frac{M}{gcd(N,M)}\)的时候成立
我们只要不断枚举\(k * i,h * i\)判断两个串里这个位置的数是否相同即可
代码
#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define mp make_pair
#define pb push_back
#define enter putchar('\n')
#define space putchar(' ')
//#define ivorysi
#define MAXN 100005
typedef long long int64;
using namespace std;
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);
}
int64 N,M;
char S[MAXN],T[MAXN];
int64 gcd(int64 a,int64 b) {
return b == 0 ? a : gcd(b,a % b);
}
int64 lcm(int64 a,int64 b) {
return a * b / gcd(a,b);
}
void Solve() {
read(N);read(M);
scanf("%s",S + 1);scanf("%s",T + 1);
int64 c = lcm(N,M);
int64 a = c / M,b = c / N;
for(int i = 0 ; i <= N ; ++i) {
if(i * a + 1 > N || i * b + 1 > M) break;
if(S[i * a + 1] != T[i * b + 1]) {puts("-1");return;}
}
out(c);enter;
}
int main() {
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Solve();
return 0;
}
B - Removing Blocks
我计数计的挺麻烦的
我是对于每一段统计一下这一段会被哪些序列统计到
显然对于长度为K的一个段，除掉最靠前的一段和最靠后的一段，中间的每个段，如果这一段被取了，那么这个操作序列必定包含一个子序列，就是这一段两侧被提前取了，中间是K个数的任意排列
然后我们对于每个长度特殊处理最前一段和最后一段（因为他们的端点只有一个），对于中间的序列计算排列个数，而中间序列的代价总和可以通过前缀和的前缀和快速算出来
代码
#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define mp make_pair
#define pb push_back
#define enter putchar('\n')
#define space putchar(' ')
//#define ivorysi
#define MAXN 100005
typedef long long int64;
using namespace std;
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;
int N,A[MAXN],fac[MAXN],invfac[MAXN],inv[MAXN],sum[MAXN],sum_of_sum[MAXN];
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 C(int n,int m) {
if(n < m) return 0;
return mul(fac[n],mul(invfac[m],invfac[n - m]));
}
void Solve() {
read(N);
for(int i = 1 ; i <= N ; ++i) {
read(A[i]);
sum[i] = inc(sum[i - 1],A[i]);
sum_of_sum[i] = inc(sum_of_sum[i - 1],sum[i]);
}
inv[1] = 1;
for(int i = 2 ; i <= N ; ++i) {
inv[i] = mul(inv[MOD % i],MOD - MOD / i);
}
fac[0] = invfac[0] = 1;
for(int i = 1 ; i <= N ; ++i) {
fac[i] = mul(fac[i - 1],i);
invfac[i] = mul(invfac[i - 1],inv[i]);
}
int ans = 0;
for(int i = 1 ; i <= N ; ++i) {
if(i == N) update(ans,mul(sum[N],fac[N]));
else {
update(ans,mul(inc(sum[i],inc(sum[N],MOD - sum[N - i])),mul(mul(fac[i],C(N,i + 1)),fac[N - i - 1])));
if(i < N - 1) {
int t = inc(sum_of_sum[N - 1],MOD - sum_of_sum[i]);
update(t,MOD - sum_of_sum[N - i - 1]);
update(ans,mul(t,mul(mul(mul(2,fac[i]),C(N,i + 2)),fac[N - i - 2])));
}
}
}
out(ans);enter;
}
int main() {
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Solve();
return 0;
}
C - Min Cost Cycle
显然如果最优的话，答案至少为这2×N个数里前N小的数的和
我们把所有前N小的数标出来，称为1，后N个数称为0
那么一个点有4种情况 01 10 00 11
后两种个数相同
如果我们有00（我们一定会有相同个数的11）的话，那么一定可行，为什么，因为我们把01首尾相接成一个01，把10首尾相接成一个10，然后01 00 10 11就可以了
如果我们只有10或只有01，也一定可行
唯一取不到这个解的是有10，并且有01，而且不存在00
那么这个时候我们枚举两个1相撞时1的值（如果这个1和另一个1相撞的时候它不是最小值也没关系，最后会被更新掉），同时计算一下两个0相撞时0能取到的最大值，可以用一个set维护，因为我们需要删掉枚举的1所带的那个0，除了它所在的类型只有1个的情况
代码
#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define mp make_pair
#define pb push_back
#define enter putchar('\n')
#define space putchar(' ')
//#define ivorysi
#define MAXN 100005
typedef long long int64;
using namespace std;
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;
int64 p[MAXN][2],sum,ans,t3;
int L[MAXN],R[MAXN],tot,ty[MAXN],cnt[10];
pii pos[MAXN * 2];
multiset<int64> S;
bool cmp(pii a,pii b) {
return p[a.fi][a.se] < p[b.fi][b.se];
}
void Solve() {
read(N);
for(int i = 1 ; i <= N ; ++i) {
read(p[i][1]);read(p[i][0]);
ans += p[i][1] + p[i][0];
pos[++tot] = mp(i,0);
pos[++tot] = mp(i,1);
}
sort(pos + 1,pos + tot + 1,cmp);
for(int i = 1 ; i <= tot / 2 ; ++i) {
sum += p[pos[i].fi][pos[i].se];
if(!pos[i].se) L[pos[i].fi] = 1;
else R[pos[i].fi] = 1;
S.insert(p[pos[i].fi][pos[i].se]);
}
for(int i = 1 ; i <= N ; ++i) {
if(L[i] && !R[i]) {ty[i] = 1;cnt[1]++;}
if(!L[i] && R[i]) {ty[i] = 2;cnt[2]++;}
if(!L[i] && !R[i]) {ty[i] = 3;cnt[3]++;}
if(L[i] && R[i]) {ty[i] = 4;cnt[4]++;}
}
if(!cnt[1] || !cnt[2] || cnt[3]) {out(sum);enter;return;}
for(int i = tot / 2 + 1 ; i <= tot ; ++i) {
int u = pos[i].fi;
if(cnt[ty[u]] != 1) {
S.erase(S.find(p[u][pos[i].se ^ 1]));
ans = min(ans,sum + p[pos[i].fi][pos[i].se] - *(--S.end()));
S.insert(p[u][pos[i].se ^ 1]);
}
else ans = min(ans,sum + p[pos[i].fi][pos[i].se] - *(--S.end()));
}
out(ans);enter;
}
int main() {
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Solve();
return 0;
}
D - Chords
这题看起来很繁琐，但是实际上逻辑很简单
\(dp[i][j]\)表示一个联通块标记为\(i\)为最小的标号，\(j\)为最大的标号，这样的方案有多少个
首先要满足初始的时候\(i,j\)之间的点没有连到外部的点
然后计算方案是\(g(x)\)
x是\(i,j\)之间没有匹配的点的个数，\(g(x) = 1*3*5...(x - 3)(x - 1)\)
然后再减掉不合法的方案
\(dp[i][j] -= dp[i][k] * g(z)\)z是\(k + 1,j\)中没有匹配过的点
最后输出的时候还要乘上一个\(g(y)\)y表示所有除了\(i,j\)之外未匹配的点
代码
#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#define enter putchar('\n')
#define space putchar(' ')
#define fi first
#define se second
#define mp make_pair
#define MAXN 200005
//#define ivorysi
#define pii pair<int,int>
using namespace std;
typedef long long int64;
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) {putchar('-');x = -x;}
if(x >= 10) out(x / 10);
putchar('0' + x % 10);
}
const int MOD = 1000000007;
int N,K;
int A[305],B[305];
int cnt[605][605],dp[605][605],g[605];
bool vis[605][605];
bool in(int a,int l,int r) {
return a >= l && a <= r;
}
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 Solve() {
read(N);read(K);
for(int i = 1 ; i <= K ; ++i) {
read(A[i]);read(B[i]);
}
g[0] = 1;
for(int i = 1 ; i <= 2 * N ; ++i) {
if(i & 1) g[i] = 0;
else g[i] = mul(g[i - 2],i - 1);
}
for(int i = 1 ; i <= 2 * N ; ++i) {
for(int j = i + 1 ; j <= 2 * N ; ++j) {
for(int h = 1 ; h <= K ; ++h) {
if(in(A[h],i,j) && in(B[h],i,j)) {
++cnt[i][j];
}
else if((!in(A[h],i,j) && in(B[h],i,j)) || (in(A[h],i,j) && !in(B[h],i,j))) {
vis[i][j] = 1;
break;
}
}
}
}
for(int d = 2 ; d <= 2 * N ; d += 2) {
for(int i = 1 ; i <= 2 * N ; ++i) {
int j = i + d - 1;
if(j > 2 * N) break;
if(vis[i][j]) continue;
int x = j - i + 1 - 2 * cnt[i][j];
dp[i][j] = g[x];
for(int h = i + 1 ; h < j ; h += 2) {
if(!dp[i][h]) continue;
x = j - h - 2 * cnt[h + 1][j];
dp[i][j] = inc(dp[i][j],MOD - mul(dp[i][h],g[x]));
}
}
}
int ans = 0;
for(int i = 1 ; i <= 2 * N ; ++i) {
for(int j = i + 1 ; j <= 2 * N ; j += 2) {
int y = 2 * N - (j - i + 1 - 2 * cnt[i][j]) - 2 * K;
ans = inc(ans,mul(dp[i][j],g[y]));
}
}
out(ans);enter;
}
int main() {
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Solve();
}
E - High Elements
我们需要已知部分\(X\)和\(Y\)的情况下，知道后面的数能不能使得最终情况合法
我们设\(X\)高元素的个数为\(C_x\)，\(Y\)高元素个数为\(C_y\)
我们还需要两个高元素序列\(a,b\)满足\(C_x + |a| = C_y + |b|\)
我们可以保证\(a,b\)至少有一个，其中所有的高元素都是p原来的高元素，证明这个，我们只要找到两个合法的序列，然后每次减少两个序列中高元素，直到不能减了
这样的话我们设\(a\)全是P原来的高元素组成的序列
然后我们决定一下\(b\)，\(b\)确定了\(a\)就自动是未选的所有高元素了
设\(b\)中的高元素有\(k\)个是P中的高元素，\(m\)是非高元素，剩下的序列中有Q个是高元素
那么我们满足
\(C_x + Q - k = C_y + |b|\)
\(C_x + Q - k = C_y + k + m\)
\(2k + m = C_x - C_y + Q\)
也就是我们设所有高元素长度为2，非高元素长度为1，求一个上升子序列，看看长度是否正好为一个常数
我们只需要分奇偶性找一个最长的升子序列，和这个常数比较一下就行
需要用线段树维护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 MAXN 200005
//#define ivorysi
#define space putchar(' ')
#define enter putchar('\n')
using namespace std;
typedef long long int64;
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;}
if(x >= 10) {
out(x / 10);
}
putchar('0' + x % 10);
}
int N,P[MAXN],val[MAXN],ans[MAXN],curx,cury;
struct node {
int maxv,L,R,lc,rc;
}tr[MAXN * 20];
int rt[2],Ncnt;
void build(int &u,int L,int R,int v) {
u = ++Ncnt;
tr[u].L = L;tr[u].R = R;
tr[u].maxv = v;
if(L == R) return;
int mid = (L + R) >> 1;
build(tr[u].lc,L,mid,v);
build(tr[u].rc,mid + 1,R,v);
}
void Change(int u,int pos,int v) {
if(tr[u].L == tr[u].R) {tr[u].maxv = v;return;}
int mid = (tr[u].L + tr[u].R) >> 1;
if(pos <= mid) Change(tr[u].lc,pos,v);
else Change(tr[u].rc,pos,v);
tr[u].maxv = max(tr[tr[u].lc].maxv,tr[tr[u].rc].maxv);
}
int Query(int u,int l,int r) {
if(tr[u].L == l && tr[u].R == r) return tr[u].maxv;
int mid = (tr[u].L + tr[u].R) >> 1;
if(r <= mid) return Query(tr[u].lc,l,r);
else if(l > mid) return Query(tr[u].rc,l,r);
else return max(Query(tr[u].lc,l,mid),Query(tr[u].rc,mid + 1,r));
}
void Solve() {
read(N);
for(int i = 1 ; i <= N ; ++i) read(P[i]);
int cur = 0,Q = 0;
for(int i = 1 ; i <= N ; ++i) {
cur = max(cur,P[i]);
if(cur == P[i]) {val[i] = 2;++Q;}
else val[i] = 1;
}
build(rt[0],1,N + 1,0);build(rt[1],1,N + 1,-2 * N);
for(int i = N ; i >= 1 ; --i) {
int tmp[2];
for(int j = 0 ; j < 2 ; ++j) tmp[j] = Query(rt[j],P[i] + 1,N + 1);
for(int j = 0 ; j < 2 ; ++j) {
if(tmp[j] + val[i] > 0) {
int t = (tmp[j] + val[i]) & 1;
Change(rt[t],P[i],tmp[j] + val[i]);
}
}
}
int len[2] = {0,0},mx_num[2] = {0,0};
for(int i = 1 ; i <= N ; ++i) {
ans[i] = -1;
Change(rt[0],P[i],0);Change(rt[1],P[i],-2 * N);
if(val[i] == 2) --Q;
for(int k = 0 ; k < 2 ; ++k) {
int nxt_len[2],nxt_mx_num[2];
memcpy(nxt_len,len,sizeof(len));
memcpy(nxt_mx_num,mx_num,sizeof(mx_num));
nxt_mx_num[k] = max(nxt_mx_num[k],P[i]);
if(nxt_mx_num[k] == P[i]) ++nxt_len[k];
if(i == N) {
if(nxt_len[0] == nxt_len[1]) {ans[i] = k;break;}
continue;
}
bool flag = 0;
for(int j = 0 ; j < 2 ; ++j) {
int t = nxt_len[j] - nxt_len[j ^ 1] + Q;
if(t >= 0 && Query(rt[t & 1],nxt_mx_num[j ^ 1] + 1,N + 1) >= t) {flag = 1;break;}
}
if(flag) {
ans[i] = k;
memcpy(len,nxt_len,sizeof(len));
memcpy(mx_num,nxt_mx_num,sizeof(mx_num));
break;
}
}
if(ans[i] < 0) {
puts("-1");return;
}
}
for(int i = 1 ; i <= N ; ++i) putchar('0' + ans[i]);
enter;
}
int main() {
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Solve();
}
【AtCoder】AGC028 (A-E)题解的更多相关文章
- AtCoder ExaWizards 2019 简要题解
AtCoder ExaWizards 2019 简要题解 Tags:题解 link:https://atcoder.jp/contests/exawizards2019 很水的一场ARC啊,随随便便就 ...
- AtCoder Beginner Contest 154 题解
人生第一场 AtCoder,纪念一下 话说年后的 AtCoder 比赛怎么这么少啊(大雾 AtCoder Beginner Contest 154 题解 A - Remaining Balls We ...
- AtCoder Beginner Contest 153 题解
目录 AtCoder Beginner Contest 153 题解 A - Serval vs Monster 题意 做法 程序 B - Common Raccoon vs Monster 题意 做 ...
- AtCoder Beginner Contest 177 题解
AtCoder Beginner Contest 177 题解 目录 AtCoder Beginner Contest 177 题解 A - Don't be late B - Substring C ...
- AtCoder Beginner Contest 184 题解
AtCoder Beginner Contest 184 题解 目录 AtCoder Beginner Contest 184 题解 A - Determinant B - Quizzes C - S ...
- AtCoder Beginner Contest 173 题解
AtCoder Beginner Contest 173 题解 目录 AtCoder Beginner Contest 173 题解 A - Payment B - Judge Status Summ ...
- AtCoder Beginner Contest 172 题解
AtCoder Beginner Contest 172 题解 目录 AtCoder Beginner Contest 172 题解 A - Calc B - Minor Change C - Tsu ...
- AtCoder Beginner Contest 169 题解
AtCoder Beginner Contest 169 题解 这场比赛比较简单,证明我没有咕咕咕的时候到了! A - Multiplication 1 没什么好说的,直接读入两个数输出乘积就好了. ...
- AtCoder Beginner Contest 148 题解
目录 AtCoder Beginner Contest 148 题解 前言 A - Round One 题意 做法 程序 B - Strings with the Same Length 题意 做法 ...
- AtCoder Grand Contest 017 题解
A - Biscuits 题目: 给出 \(n\) 个物品,每个物品有一个权值. 问有多少种选取方式使得物品权值之和 \(\bmod\space 2\) 为 \(p\). \(n \leq 50\) ...
随机推荐
- MT【103】二阶递推找规律
评:如果直接找$a_n$的二阶递推式:$a_{n+2}-2\sqrt{2}a_{n+1}-a_n=0$有根号,不利于估计尾数.
- springboot中radis配置和使用【进阶二】
1.yml文件配置 #redis哨兵模式配置 redis: namespace: xxx:xxx:redis connection: cacheRedis: database: 0 timeout: ...
- [SDOI2017]硬币游戏
考虑生成函数来做 g(x)函数就是0+0*x+...+1*x^s+...+|∑|^(n-s)x^n 就是最后s位必须填这个串,但是前面随便填的方案数 然后枚举之前出现了哪个串(包括自己),如果没有相交 ...
- 【CF522A】Reposts
题目大意:给定一个有向图,求图中最长路. 题解:直接拓扑排序后按照拓扑序枚举即可.处理时应将字符串通过 map 映射成一个点,同时注意字符串大小写转换,C++ string 中没有提供直接大小写转换的 ...
- java Random.nextInt()方法
转: java Random.nextInt()方法 lic int nextInt(int n) 该方法的作用是生成一个随机的int值,该值介于[0,n)的区间,也就是0到n之间的随机int值,包含 ...
- C# 编码规范、命名规则
1 规范目的 ……………………………………………………… 3 2 适用范围 ……………………………………………………… 3 3 代码注释 ……………………………………………………… 3 3.1 ...
- Spark记录-Scala异常处理与文件I/O
Scala的异常处理类似许多其他语言(如Java).它不是以正常方式返回值,方法可以通过抛出异常来终止. 但是,Scala实际上并没有检查异常. 当您想要处理异常时,要像Java一样使用try {.. ...
- HTML5页面开发的基础性模板
分享一个HTML5页面开发的基础性模板,包含了两个版本: 开发版本 注释版本 开发版本 <!DOCTYPE html> <html> <head> <meta ...
- javascript私有方法揭示为公有方法
揭示模式可用于将私有方法暴露为公共方法.当为了对象的运转而将所有功能放置在一个对象中以及想尽可能地保护该对象,这种揭示模式就非常有用. 板栗: var myarray; (function(){ va ...
- CSS3 - chrome，傲游，360极速浏览器不支持小于12px的字号的解决办法
google流量器chrome,傲游,360极速浏览器都是基于webkit内核浏览器,默认不支持小于font-size小于12px 的字号,即定义font-size小于12px时会发现字体大小依然是1 ...