Atcoder AGC019简要题解

最新推荐文章于 2022-12-20 14:54:48 发布
最新推荐文章于 2022-12-20 14:54:48 发布
阅读量382 收藏
点赞数
分类专栏： contest
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_35649707/article/details/84333388
版权

传送门

Reverse and Compare

发现两个串如果翻转后一样且中心不一样，则必定有一个是回文串。然后就只用统计两端字母不同的串的个数了。

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
 
const int N=2e5+50;
char ch[N];
int n,cnt[26],s;
LL ans;
int main() {
	scanf("%s",ch+1); n=strlen(ch+1);
	for(int i=1;i<=n;i++)
		ans+=s-cnt[ch[i]-'a'], ++s, ++cnt[ch[i]-'a'];
	cout<<ans+1<<'\n';
}

Fountain Walk

这道题想了半天，直到我看见每行每列最多一个喷泉。。。

然后就是求个最多能经过多少个了，也就是求个LIS。

#include <bits/stdc++.h>
using namespace std;
 
const int RLEN=1<<18|1;
inline char nc() {
	static char ibuf[RLEN],*ib,*ob;
	(ib==ob) && (ob=(ib=ibuf)+fread(ibuf,1,RLEN,stdin));
	return (ib==ob) ? -1 : *ib++;
}
inline int rd() {
	char ch=nc(); int i=0,f=1;
	while(!isdigit(ch)) {if(ch=='-')f=-1; ch=nc();}
	while(isdigit(ch)) {i=(i<<1)+(i<<3)+ch-'0'; ch=nc(); }
	return i*f;
}
 
const int N=2e5+50;
struct P {int x,y;} a[N];
int sx,sy,tx,ty,n,d[N],top;
int main() {
	sx=rd(), sy=rd(), tx=rd(), ty=rd(), n=rd();
	for(int i=1;i<=n;i++) {
		a[i].x=rd(), a[i].y=rd();
		if(a[i].x<min(sx,tx) || a[i].x>max(sx,tx) || a[i].y<min(sy,ty) || a[i].y>max(sy,ty)) --i, --n;
	}
	if(sx>tx) {
		sx=-sx; tx=-tx;
		for(int i=1;i<=n;i++) a[i].x=-a[i].x;
	}
	if(sy>ty) {
		sy=-sy; ty=-ty;
		for(int i=1;i<=n;i++) a[i].y=-a[i].y;
	}
	sort(a+1,a+n+1,[](const P &a,const P &b) {return a.x<b.x;});
	for(int i=1;i<=n;i++) {
		if((!top) || d[top]<a[i].y) d[++top]=a[i].y;
		else d[lower_bound(d+1,d+top+1,a[i].y)-d]=a[i].y;
	} 
	const double PI=acos(-1.0);
	double ans=(tx-sx+ty-sy)*100.0+top*(5.0*PI-20);
	if(top==min(tx-sx+1,ty-sy+1)) ans+=5.0*PI;
	printf("%.20f\n",ans);
}

Shift and Flip

枚举最后哪里开始相同，然后 O ( n log ⁡ n ) O(n \log n) O(nlogn)做个贪心，复杂度 O ( n 2 log ⁡ n ) O(n^2 \log n) O(n2logn)。

#include <bits/stdc++.h>
using namespace std;
typedef pair <int,int> pii;
const int RLEN=1<<18|1;
inline char nc() {
	static char ibuf[RLEN],*ib,*ob;
	(ib==ob) && (ob=(ib=ibuf)+fread(ibuf,1,RLEN,stdin));
	return (ib==ob) ? -1 : *ib++;
}
inline int rd() {
	char ch=nc(); int i=0,f=1;
	while(!isdigit(ch)) {if(ch=='-')f=-1; ch=nc();}
	while(isdigit(ch)) {i=(i<<1)+(i<<3)+ch-'0'; ch=nc();}
	return i*f;
}
 
const int N=2e3+50;
int n,pre[N],suf[N];
char a[N],b[N];
inline int dis(int a,int b) {
	if(a<=b) return b-a;
	else return b+n-a;
}
inline int solve(int p) {
	vector <pii> vec;
	for(int i=0;i<n;i++) {
		int tar=(i+p)%n;
		if(a[tar]^b[i]) 
			vec.push_back(pii(dis(pre[tar],tar),dis(tar,suf[tar])));
	}
	if(!vec.size()) return min(dis(0,p),dis(p,0));
	int ans=1e9;
	sort(vec.begin(),vec.end(),[](const pii &a,const pii &b) {return a.first<b.first || (a.first==b.first && a.second>b.second);});
	for(int i=vec.size()-1,mx=0;~i;i--) {
		mx=max(mx,vec[i].second);
		int mx2=(i==0) ? 0 : vec[i-1].first;
		ans=min(ans,mx2+mx2+mx+min(dis((p+mx)%n,0),dis(0,(p+mx)%n)));
		ans=min(ans,mx+mx+mx2+min(dis((p-mx2+n)%n,0),dis(0,(p-mx2+n)%n)));
	} return ans+vec.size();
}
int main() {
	scanf("%s%s",a,b);
	n=strlen(a);
	int s1=0, s2=0;
	for(int i=0;i<n;i++)
		s1+=a[i]-'0', s2+=b[i]-'0';
	if(s1 && !s2) {return puts("-1"),0;}
	
	if(b[0]=='1') pre[0]=0;
	else {
		for(int i=n-1;~i;i--)
			if(b[i]=='1') {pre[0]=i; break;}
	}
	for(int i=1;i<n;i++)
		pre[i]=(b[i]=='1' ? i : pre[i-1]);
	if(b[n-1]=='1') suf[n-1]=n-1;
	else {
		for(int i=0;i<n;i++) 
			if(b[i]=='1') {suf[n-1]=i; break;} 
	}
	for(int i=n-2;~i;i--)
		suf[i]=(b[i]=='1' ? i : suf[i+1]);
	int ans=1e9;
	for(int i=0;i<n;i++)
		ans=min(ans,solve(i));
		
	reverse(a,a+n); reverse(b,b+n);
	if(b[0]=='1') pre[0]=0;
	else {
		for(int i=n-1;~i;i--)
			if(b[i]=='1') {pre[0]=i; break;}
	}
	for(int i=1;i<n;i++)
		pre[i]=(b[i]=='1' ? i : pre[i-1]);
	if(b[n-1]=='1') suf[n-1]=n-1;
	else {
		for(int i=0;i<n;i++) 
			if(b[i]=='1') {suf[n-1]=i; break;} 
	}
	for(int i=n-2;~i;i--)
		suf[i]=(b[i]=='1' ? i : suf[i+1]);
	for(int i=0;i<n;i++)
		ans=min(ans,solve(i));
	
	cout<<ans<<'\n';
}

Shuffle and Swap

发现合法的就是那些不相等的位置中间插入一些相等的位置。

把这个指数生成函数搞出来做个FFT快速幂即可，时间复杂度 O ( n log ⁡ 2 n ) O(n \log^2 n) O(nlog2n)。

#include <bits/stdc++.h>
using namespace std;
 
const int N=1e6+50, mod=998244353;
inline int add(int x,int y) {return (x+y>=mod) ? (x+y-mod) : (x+y);}
inline int dec(int x,int y) {return (x-y<0) ? (x-y+mod) : (x-y);}
inline int mul(int x,int y) {return (long long)x*y%mod;}
inline int power(int a,int b,int rs=1) {for(;b;b>>=1,a=mul(a,a)) if(b&1) rs=mul(rs,a); return rs;}
inline int sgn(int x) {return (x&1) ? (mod-1) : 1;}
 
namespace FFT {
	const int G=3;
	int A[N],B[N],w[N],pos[N],k;
	inline void init(int n) {
		for(k=1;k<=n;k<<=1);
		memset(A,0,sizeof(int)*k);
		memset(B,0,sizeof(int)*k);
		for(int i=1;i<k;i++) pos[i]=(i&1) ? ((pos[i>>1]>>1)^(k>>1)) : (pos[i>>1]>>1);
	}
	inline void dft(int *a) {
		for(int i=1;i<k;i++)
			if(pos[i]>i) swap(a[pos[i]],a[i]);
		for(int bl=1;bl<k;bl<<=1) {
			int tl=bl<<1, wn=power(G,(mod-1)/tl);
			w[0]=1; for(int i=1;i<bl;i++) w[i]=mul(w[i-1],wn);
			for(int bg=0;bg<k;bg+=tl)
				for(int j=0;j<bl;j++) {
					int &t1=a[bg+j], &t2=a[bg+j+bl], t=mul(t2,w[j]);
					t2=dec(t1,t); t1=add(t1,t);
				}
		}
	}
	inline void func() {
		dft(A); dft(B);
		for(int i=0;i<k;i++) B[i]=mul(B[i],A[i]);
		dft(B); const int inv=power(k,mod-2);
		for(int i=0;i<k;i++) B[i]=mul(B[i],inv);
		reverse(B+1,B+k); 
	}
}
struct combin {
	int fac[N],ifac[N];
	combin() {
		fac[0]=1;
		for(int i=1;i<N;i++) fac[i]=mul(fac[i-1],i);
		ifac[0]=ifac[1]=1;
		for(int i=2;i<N;i++) ifac[i]=mul(mod-mod/i,ifac[mod%i]);
		for(int i=2;i<N;i++) ifac[i]=mul(ifac[i-1],ifac[i]); 
	}
	inline int C(int a,int b) {return mul(fac[a],mul(ifac[b],ifac[a-b]));}
} C;
 
struct poly {
	vector <int> a;
	poly(int d=0,int t=0) {a.resize(d+1); a[d]=t;}
	inline int& operator [](const int &i) {return a[i];}
	inline const int& operator [](const int &i) const {return a[i];}
	inline int deg() const {return a.size()-1;}
	inline poly extend(int k) {poly c=*this; c.a.resize(k); return c;}
	friend inline poly operator *(const poly &a,const poly &b) {
		poly c(a.deg()+b.deg()); FFT::init(c.deg());
		for(int i=0;i<=a.deg();i++) FFT::A[i]=a[i];
		for(int i=0;i<=b.deg();i++) FFT::B[i]=b[i];
		FFT::func();
		for(int i=0;i<=c.deg();i++) c[i]=FFT::B[i];
		return c;
	}
	inline void pt() {
		for(int i=0;i<=deg();i++) cerr<<a[i]<<' '; cerr<<'\n';
	}
};
 
int n,m,l;
char a[N],b[N];
int main() {
	scanf("%s%s",a,b);
	l=strlen(a);
	for(int i=0;i<l;i++)
		if(a[i]=='1' &&a[i]==b[i]) ++n;
		else if(a[i]=='1') ++m;
	poly g(n+1,0); 
	for(int i=1;i<=n+1;i++) g[i]=C.ifac[i];
	poly f(n+m,0); f[0]=1;
	for(int b=m;b;b>>=1,g=(g*g).extend(n+m+1))
		if(b&1) f=(f*g).extend(n+m+1);
	int ans=0;
	for(int i=m;i<=n+m;i++) {
		int res=mul(mul(f[i],C.fac[i]),C.C(n+m,i));
		res=mul(res,mul(C.fac[n+m-i],C.fac[n+m-i]));
		res=mul(res,C.C(n,i-m));
		res=mul(res,C.fac[i-m]);
		ans=add(ans,res);
	}
	cout<<mul(ans,C.fac[m]);
}

Yes or No

考虑一个方案，发现会被计算 max ⁡ { n , m } + 对 角 线 经 过 次 数 2 \max\{n,m\}+\frac{对角线经过次数}{2} max{n,m}+2对角线经过次数​次，那我们对对角线上每个点算一算会被经过多少次就行了。

#include <bits/stdc++.h>
using namespace std;

const int N=1e6+50, mod=998244353;
inline int add(int x,int y) {return (x+y>=mod) ? (x+y-mod) : (x+y);}
inline int dec(int x,int y) {return (x-y<0) ? (x-y+mod) : (x-y);}
inline int mul(int x,int y) {return (long long)x*y%mod;}
inline int power(int a,int b,int rs=1) {for(;b;b>>=1,a=mul(a,a)) if(b&1) rs=mul(rs,a); return rs;}
inline int inv(int a) {return power(a,mod-2);}

struct combin {
	int fac[N],ifac[N];
	combin() {
		fac[0]=1;
		for(int i=1;i<N;i++) fac[i]=mul(fac[i-1],i);
		ifac[0]=ifac[1]=1;
		for(int i=2;i<N;i++) ifac[i]=mul(mod-mod/i,ifac[mod%i]);
		for(int i=2;i<N;i++) ifac[i]=mul(ifac[i-1],ifac[i]);
	}
	inline int C(int a,int b) {return mul(fac[a],mul(ifac[b],ifac[a-b]));}
} C;
#define x1 dqoaspjdioq
#define x2 dnasuihquih
#define y1 adsjiohejqu
#define y2 dnasuigqyuw
inline int path(int x1,int y1,int x2,int y2) {
	return C.C(x2-x1+y2-y1,x2-x1);
}
int n,m;
int main() {
	cin>>n>>m; if(n>m) swap(n,m);
	int ans=0;
	for(int i=1;i<=n;i++)
		ans=add(ans,mul(path(0,0,i,i),path(i,i,n,m)));
	ans=mul(ans,inv(2));
	ans=mul(ans,inv(path(0,0,n,m)));
	cout<<add(m,ans);
}
DZYO
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
AtCoder abc352 A题题解视频
05-05
AtCoder abc352 A题题解视频
AtCoder abc352-C题题解
最新发布
05-05
AtCoder abc352-C题题解
【Atcoder】AGC019 B-F简要题解
远行客
03-21 202
D.套路枚举 E.转图论DP F.+-1转网格图DP期望
AtCoder AGC019E Shuffle and Swap (DP、FFT、多项式求逆、多项式快速幂)
suncongbo's blog
09-11 208
题目链接 https://atcoder.jp/contests/agc019/tasks/agc019_e 题解 tourist的神仙E题啊做不来做不来……这题我好像想歪了啊= =…… 首先我们可以考虑，什么样的操作序列才是合法的？ 有用的位置只有两种，一种是两个序列在这个位置上都是1, 称作11型，另一种是一个0一个1, 称作01型。设两种位置分别有\(A\)个和\(2B\)个。 考...
[AGC019E]Shuffle and Swap
weixin_33790053的博客
06-11 133
[AGC019E]Shuffle and Swap 题目大意： 给出两个长度为\(n(n\le10000)\)的\(01\)串\(A_{1\sim n}\)和\(B_{1\sim n}\)。两个串均有\(k\)个'1'。令\(a_{1\sim k}\)和\(b_{1\sim k}\)分别表示\(A\)和\(B\)中所有'1'出现的位置。将\(a\)和\(b\)等概率随机排列，按\(1\sim k\...
[agc019e]Shuffle and Swap
不来也不去的一只失忆蝴蝶
10-11 635
前言计数水平不行了。 居然不会n^2去dp这个模型。题意太麻烦了看网站。DP假设有x个公用1和y个非公用1。 用图论来理解。 最后一定形成y条链，若干个环。 其中链中的边在序列中要按顺序，环则随意。 环可以先不管。 考虑dp[i,j]表示目前做出来两个i+j的序列，有i个公用1和j个非公用1，区分编号，有顺序，形成j条链的方案数。 有dp[i,j]=dp[i-1,j]*i*j+dp[i
AtCoder abc352-B题题解
05-05
AtCoder abc352-B题题解
AtCoder Regular Contest 102 Task C 题解（含代码）
01-16
【内测版】思维专项 1.0 Task O 一句话题解 by LJ24TangHaoLong「唐浩滝」 题面： 求a, b, c不大于n且满足a+b, b+c, a+c为k倍数的三元组(a, b, c)的个数 原题： AtCoder Regular Contest 102 Task C
Atcoder题解与视频集(2022.05.23).pdf
05-23
Atcoder 题解与视频集 Atcoder 是一个在线编程平台，提供了大量的编程题目和比赛项目，涵盖了算法、数据结构、数学、计算机科学等多个领域。以下是 Atcoder 题解与视频集的知识点总结： 1. Atcoder 的基本概念：...
atcoder:提交给AtCoder的代码
04-08
存档AtCoder接受的（AC）代码。
Atcoder ARC102测试数据
10-12
日本竞赛网站Atcoder的ARC102场的测试数据，题面可以在Atcoder找到。
[agc019f]Yes or No
ZLTJohn的博客
02-16 461
题目大意 给出n+m个问题，其中n个答案是yes，m个no，每一次回答问题，你选择回答yes或no，你能知道当前局面未回答问题的yes和no的数量。问最后你期望回答对多少个问题。 n,m≤5e5,在模998244353意义下计算答案。 解题思路 设剩下a个yes，b个no的局面是(a,b) 首先手玩一下这个问题，很容易列出一个dp转移式，然后发现最优策略：对于(a,b)，我们肯定是猜数量...
AGC019 E.Shuffle and Swap-DP+NTT
Fizzmy的博客
03-04 469
传送门 题意： 给出两个01串A,b，记aiaia_i表示A中1的出现位置，bibib_i表示B中1的出现位置，将a数组和b数组打乱后依次次交换AaiAaiA_{a_i}和AbiAbiA_{b_i}，求有几种方式使得A=B 字符串长度&lt;=10000 Solution： 我们可以把答案拆分成两步： 1.枚举a和b的匹配 2.打乱匹配顺序 假设我们已经完成了操作1，我们来...
【DP】AGC019E Shuffle and Swap
氧化钠的博客
10-09 368
分析： 可以尝试构造每一个合法数列。 显然，如果某个位置，Ai=1A_i=1Ai​=1且Bi=1B_i=1Bi​=1，那么这个A中的1被拿走后，必须有一个1填回来。此类点定义为公共点。 如果某个位置，Ai=1且Bi=0A_i=1且B_i=0Ai​=1且Bi​=0，那么这个A中的1被拿走后，不能有其他1填进来。此类点定义为非公共点。对应的，称Ai=0且Bi=1A_i=0且B_i=1Ai​=0且Bi​...
【概率与期望】AGC019 F Yes or No
氧化钠的博客
10-09 448
分析： 首先，很容易搞出来一个N^2的DP做法，显然会T，但是对发现正解有很大帮助： 将这个DP转换为一个网格，就变为：从起始点(n,m)出发，到达目标点(0,0)的期望路径长度。 首先，有一个很厉害的结论：因为每次都是按照最优策略选答案，所以不妨设n&gt;=m，那么一定能答对n次。证明很显然，每次当n&gt;m时我们都会n的那一项，尽管答案可能会错，但错了之后，n会更大，我们更会选n。直到最终...
AGC019F Yes or No
LazyCrazyCat的博客
01-04 181
题意 有N+M个问题，其中N个答案是YES，M个是NO 每次你可以回答YES或NO，问答对个数的期望 当然还有个条件，当你回答了某个问题后，你会得知这道题的正确答案。所以问你采取最优方案时的期望 题解 我们先考虑假如不会知道答案 当N&amp;gt;M时，显然我们就应该一直答Yes，这样我们答对的期望就是(N+M)NN+M=N(N+M)\frac{N}{N+M}=N(N+M)N+MN​=N 反之亦然，所以...
如何使用Python爬虫向AtCoder提交代码？
weixin_42589700的博客
12-20 140
为了使用 Python 爬虫向 AtCoder 提交代码，您需要执行以下步骤： 安装所需的库。要使用 Python 爬虫向 AtCoder 提交代码，您需要安装至少两个库：requests 和 BeautifulSoup。你可以使用 pip install requests beautifulsoup4 命令安装这些库。 登录 AtCoder。要提交代码，您需要先登录 AtCoder。可以使用...
AtCoder Beginner Contest 134 简要题解
05-23
AtCoder Beginner Contest 134 是一场 AtCoder 的入门级比赛，以下是每道题的简要题解： A - Dodecagon 题目描述：已知一个正十二边形的边长，求它的面积。 解题思路：正十二边形的内角为 $150^\circ$，因此可以...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值