【AtCoder】ARC100 题解
C - Linear Approximation
找出\(A_i - i\)的中位数作为\(b\)即可
题解
#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);
}
int64 A[MAXN];
int N;
void Solve() {
read(N);
for(int i = 1 ; i <= N ; ++i) {
read(A[i]);
A[i] -= i;
}
sort(A + 1,A + N + 1);
int64 t = A[N / 2 + 1],ans = 0;
for(int i = 1 ; i <= N ; ++i) {
ans += abs(A[i] - t);
}
out(ans);enter;
}
int main() {
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Solve();
}
D - Equal Cut
枚举2 和 3中间的位置,两边都必须切成绝对值相差最小才能使总体绝对值相差最小,切的位置不断右移,直接两个指针扫一遍即可
#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 200005
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 a[MAXN],sum[MAXN];
int64 get_abs(int l1,int r1,int l2,int r2) {
return abs((sum[r1] - sum[l1 - 1]) - (sum[r2] - sum[l2 - 1]));
}
void Solve() {
read(N);
for(int i = 1 ; i <= N ; ++i) {
read(a[i]);sum[i] = sum[i - 1] + a[i];
}
int l = 1,p = 2,r = p + 1;
int64 ans = sum[N];
while(p <= N - 2) {
while(l + 1 < p && get_abs(1,l,l + 1,p) > get_abs(1,l + 1,l + 2,p)) ++l;
r = max(r,p + 1);
while(r + 1 < N && get_abs(p + 1,r,r + 1,N) > get_abs(p + 1,r + 1,r + 2,N)) ++r;
int64 m[] = {sum[l],sum[p] - sum[l],sum[r] - sum[p],sum[N] - sum[r]};
int64 tmp = 0;
for(int i = 0 ; i <= 3 ; ++i) {
for(int j = i + 1 ; j <= 3 ; ++j) {
tmp = max(tmp,abs(m[j] - m[i]));
}
}
ans = min(ans,tmp);
++p;
}
out(ans);enter;
}
int main() {
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Solve();
return 0;
}
E - Or Plus Max
我们转化一下问题,处理出每个i or j正好是k的子集的答案,再处理成前缀max即可
这样的话类似FMT的更新每个子集里的最大和次大即可
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <cmath>
#include <queue>
#include <ctime>
#include <map>
#include <set>
#define fi first
#define se second
#define pii pair<int,int>
//#define ivorysi
#define mp make_pair
#define pb push_back
#define enter putchar('\n')
#define space putchar(' ')
#define MAXN 100005
using namespace std;
typedef long long int64;
typedef double db;
typedef unsigned int u32;
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 - '0' + c;
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;
pii f[(1 << 18) + 5];
int ans[(1 << 18) + 5];
pii Merge(pii a,pii b) {
if(a.fi < b.fi) swap(a,b);
return mp(a.fi,max(a.se,b.fi));
}
void Solve() {
read(N);
for(int i = 0 ; i < (1 << N) ; ++i) {
read(f[i].fi);f[i].se = 0;
}
for(int i = 1 ; i < (1 << N) ; i <<= 1) {
for(int j = 0 ; j < (1 << N) ; ++j) {
if(j & i) {
f[j] = Merge(f[j],f[j ^ i]);
}
}
}
for(int i = 1 ; i < (1 << N) ; ++i) {
ans[i] = f[i].fi + f[i].se;
ans[i] = max(ans[i - 1],ans[i]);
out(ans[i]);enter;
}
}
int main() {
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Solve();
}
F - Colorful Sequences
我不会数数啊QAQ
先求出所有的序列里M这一段出现的次数的总和
答案是\((N - M + 1)K^{N - M}\)
然后求M这一段出现在不多彩的序列里次数的总和
如果M已经是多彩的了,那么答案是0
如果M不是多彩的且没有重复的数字
那么求所有N长的序列里M长含有不同数字的连续子段有多少个,答案除上\(\frac{K!}{(K - M)!}\)
那么记录dp[i][j]作为第i个,然后前j个数都是互不相同的数,j+1开始出现重复
更新的时候从dp[i - 1][h]更新
\(\left\{\begin{matrix}
1 & h \geq j \\
K - h & h = j - 1\\
0 & h < j - 1
\end{matrix}\right.\)
然后用前缀和处理可以做到\(O(NK)\)
用cnt[i][j]表示i长的序列里,倒数j个数都是互不相同的数,M长含有不同数字的连续子段有多少个
在每遇到一个dp[i][j]的j>=M就累加上
剩下转移方式类似
如果M是多彩的且有重复数字
那么就记录F为前缀最多到哪是互不相同的,B为后缀最多到哪是互不相同的
枚举M所在位置用类似的dp转移
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <vector>
#include <cstring>
//#define ivorysi
#define fi first
#define se second
#define MAXN 25005
#define enter putchar('\n')
#define space putchar(' ')
typedef long long ll;
using namespace std;
template <class T>
void read(T &res) {
res = 0;char c = getchar();T f = 1;
while(c < '0' || c > '9') {
c = getchar();
if(c == '-') f = -1;
}
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);
}
const int MOD = 1000000007;
int N,K,M;
int A[MAXN],fac[MAXN],invfac[MAXN],inv[MAXN];
int F,B,L,vis[405];
int dp[MAXN][405],cnt[MAXN][405],sum[405],sum_cnt[405],f[MAXN],b[MAXN];
int mul(int a,int b) {
return 1LL * a * b % MOD;
}
int inc(int a,int b) {
return a + b >= MOD ? a + b - MOD : a + b;
}
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 Init() {
read(N);read(K);read(M);
for(int i = 1 ; i <= M ; ++i) read(A[i]);
inv[1] = 1;
for(int i = 2 ; i <= K ; ++i) inv[i] = mul(inv[MOD % i],MOD - MOD / i);
fac[0] = invfac[0] = 1;
for(int i = 1 ; i <= K ; ++i) {
fac[i] = mul(fac[i - 1],i);
invfac[i] = mul(invfac[i - 1],inv[i]);
}
F = 0;B = 0;
memset(vis,0,sizeof(vis));
for(int i = 1 ; i <= M ; ++i) {
if(!vis[A[i]]) {
++F;
vis[A[i]] = 1;
}
else break;
}
memset(vis,0,sizeof(vis));
for(int i = M ; i >= 1 ; --i) {
if(!vis[A[i]]) {
++B;
vis[A[i]] = 1;
}
else break;
}
memset(vis,0,sizeof(vis));
int l = 0;
for(int i = 1 ; i <= M ; ++i) {
l = max(l,vis[A[i]]);
L = max(L,i - l);
vis[A[i]] = i;
}
}
void Process(int st,int *a) {
memset(dp,0,sizeof(dp));
dp[0][st] = 1;
memset(sum,0,sizeof(sum));
for(int i = st ; i <= K ; ++i) sum[i] = 1;
a[0] = 1;
for(int i = 1 ; i <= N ; ++i) {
for(int j = 1 ; j < K ; ++j) {
dp[i][j] = inc(dp[i][j],mul(dp[i - 1][j - 1],(K - j + 1)));
dp[i][j] = inc(dp[i][j],inc(sum[K],MOD - sum[j - 1]));
}
for(int j = 1 ; j <= K ; ++j) {
sum[j] = inc(sum[j - 1],dp[i][j]);
}
a[i] = sum[K - 1];
}
}
void Solve() {
if(L == K) {
out(mul(N - M + 1,fpow(K,N - M)));enter;
}
else if(F == M) {
dp[0][0] = 1;
int ans = mul(N - M + 1,fpow(K,N - M)),tmp = 0;
for(int i = 1 ; i <= N ; ++i) {
for(int j = 1 ; j < K ; ++j) {
dp[i][j] = inc(dp[i][j],mul(dp[i - 1][j - 1],(K - j + 1)));
dp[i][j] = inc(dp[i][j],inc(sum[K],MOD - sum[j - 1]));
cnt[i][j] = inc(cnt[i][j],mul(cnt[i - 1][j - 1],(K - j + 1)));
cnt[i][j] = inc(cnt[i][j],inc(sum_cnt[K],MOD - sum_cnt[j - 1]));
if(j >= M) cnt[i][j] = inc(cnt[i][j],dp[i][j]);
}
for(int j = 1 ; j <= K ; ++j) {
sum[j] = inc(sum[j - 1],dp[i][j]);
sum_cnt[j] = inc(sum_cnt[j - 1],cnt[i][j]);
}
}
for(int j = 0 ; j < K ; ++j) {
tmp = inc(tmp,cnt[N][j]);
}
tmp = mul(tmp,fpow(mul(fac[K],invfac[K - M]),MOD - 2));
ans = inc(ans,MOD - tmp);
out(ans);enter;
}
else {
Process(F,f);Process(B,b);
int ans = mul(N - M + 1,fpow(K,N - M));
for(int i = 1 ; i <= N - M + 1 ; ++i) {
int j = i + M - 1;
ans = inc(ans,MOD - mul(f[i - 1],b[N - j]));
}
out(ans);enter;
}
}
int main() {
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Init();
Solve();
}
【AtCoder】ARC100 题解的更多相关文章
- AtCoder ExaWizards2019题解
AtCoder ExaWizards2019题解 AtCoder (因为代码直接用模板写的,可能有点冗长) A.Regular Triangle 给你三根棍子的长度,问你能否用他们组成等边三角形. 什 ...
- KYOCERA Programming Contest 2021(AtCoder Beginner Contest 200) 题解
KYOCERA Programming Contest 2021(AtCoder Beginner Contest 200) 题解 哦淦我已经菜到被ABC吊打了. A - Century 首先把当前年 ...
- AT2370 Piling Up
https://www.luogu.org/jump/atcoder/2370 题解 答案不是\(2^{2m}\)因为每轮的第一次取球可能会不够. 我们可以设\(dp[i][j]\)表示到了第\(i\ ...
- Triple Shift
来源:Atcoder ARC 136 B - Triple Shift (atcoder.jp) 题解:这道题我们不可能去硬模拟(大多数这种题都不能这样去模拟的),然后我们就要去发现特性, 发现把 a ...
- 重修 Slope Trick(看这篇绝对够!)
Slope Trick 算法存在十余载了,但是我没有找到多少拍手叫好的讲解 blog,所以凭借本人粗拙的理解来写这篇文章. 本文除标明外所有图片均为本人手绘(若丑见谅),画图真的不容易啊 qwq(无耻 ...
- AtCoder Regular Contest 094 (ARC094) CDE题解
原文链接http://www.cnblogs.com/zhouzhendong/p/8735114.html $AtCoder\ Regular\ Contest\ 094(ARC094)\ CDE$ ...
- AtCoder ExaWizards 2019 简要题解
AtCoder ExaWizards 2019 简要题解 Tags:题解 link:https://atcoder.jp/contests/exawizards2019 很水的一场ARC啊,随随便便就 ...
- AtCoder Grand Contest 017 题解
A - Biscuits 题目: 给出 \(n\) 个物品,每个物品有一个权值. 问有多少种选取方式使得物品权值之和 \(\bmod\space 2\) 为 \(p\). \(n \leq 50\) ...
- AtCoder Regular Contest 077 被虐记&题解
直到\(7:58\)才知道今天\(8:00\)有\(AtCoder\)的菜鸡来写题解啦. C - pushpush 题目: 给定一个长为\(n\)的序列,第\(i\)次操作做如下的事 : 将\(a_i ...
随机推荐
- 洛谷P3275 [SCOI2011]糖果(差分约束,最长路,Tarjan,拓扑排序)
洛谷题目传送门 差分约束模板题,等于双向连0边,小于等于单向连0边,小于单向连1边,我太蒻了,总喜欢正边权跑最长路...... 看遍了讨论版,我是真的不敢再入复杂度有点超级伪的SPFA的坑了 为了保证 ...
- 杭高OI20190125 (genies出题)
/* 当一个人先从自己的内心开始奋斗,他就开始迈向了成功 ——genies (朝阳的二愣子) */ HGOI寒假赛第一场,欢迎来自各种学校的各式各样的巨老233333 感觉自己好渺小.还是NOIP ( ...
- HGOI20180813 (NOIP2018 提高组 Day2 模拟试题)
省常中省选提高Day2 继续 第一题就考了贪心,正解95pts的贪心策略第一印象是想到的,但是被自己否定掉了qwq,然后打了 不是正解的贪心,样例5没过(可怜)思路如下:先找出每个门对应可以通过的人数 ...
- 解决 Microsoft Excel has stopped working
安装了project 2013后, excel 2013 一打开文件就报这个错,解决方法如下: Step1: Try to repair Office 2013 installation and ch ...
- BZOJ 1031 [JSOI2007]字符加密Cipher 后缀数组教程
1031: [JSOI2007]字符加密Cipher Description 喜欢钻研问题的JS同学,最近又迷上了对加密方法的思考.一天,他突然想出了一种他认为是终极的加密办法:把需要加密的信息排成一 ...
- Eclipse 报 "The builder launch configuration could not be found" 错误
Eclipse 报 "The builder launch configuration could not be found" 错误的解决办法 标签: eclipseEclipse ...
- (转)Python中的模块循环导入问题
本文转自: https://wiki.woodpecker.org.cn/moin/MiscItems/2008-11-25 问题 cleven <shenglipang@gmail.com&g ...
- Bugfree——CentOS6.8搭建测试环境
参考资料:http://blog.csdn.net/qq_29227939/article/details/52295917 BugFree基于PHP和MySQL开发,是免费且开发源代码的缺陷管理系统 ...
- postman提取接口的返回值及动态设置变量(一)
一.提取接口返回值 1.当返回值是返回JSON时 let json = JSON.parse(responseBody); // responseBody是包含整个返回内容的字符串 let foo ...
- Spring Mvc + Maven + yuicompressor 使用 profile 来压缩 javascript ,css 文件; (十)
profile相关知识点: 在开发项目时,设想有以下场景: 你的Maven项目存放在一个远程代码库中(比如github),该项目需要访问数据库,你有两台电脑,一台是Linux,一台是Mac OS X, ...