Atcoder 水题选做
为什么是水题选做呢？因为我只会水题啊 (
为什么是$Atcoder$呢？因为暑假学长来讲课的时候讲了三件事：不要用洛谷，不要用dev-c++，不要用单步调试。$bzoj$太难了,$Topcoder$整了好久没学会用,中国人要长期打$codeforces$那作息简直反人类,所以就选到$atcoder$啦.而且$AT$的英语感觉比较好懂,或许是因为都是亚洲人,思维比较同步?
那么哪些题目会被放到这个地方呢？$ABC$的后两题，也就是$ARC$的前两题(如果同时举办的话).至于$ARC$的后两题以及$AGC$的题,以后可能会单开一篇来写吧.似乎每当$ABC$和$ARC$一起办的时候就稍微难一点,只办$ABC$的时候题就非常水.
做了两套题后感觉有些$C$的思路还是比较妙的，最好评的是即使是$D$也不会有很高的代码难度。前段时间掉进了大数据结构&&实现麻烦的题目的大坑,现在发现自己的思维好差啊...以后要多做一点思维题.
107 C:https://arc101.contest.atcoder.jp/tasks/arc101_a
题意概述：数轴上放有$n$根蜡烛,初始时站在$0$坐标,移动的速度是一个定值,求最少用多少时间可以点亮$k$根蜡烛.
尺取法+分类讨论,如果一段区间横跨$0$坐标,那就先走离$0$近的部分,再折回来点亮另一边的,如果与$0$没有公共部分,直接一次走过去即可.
1 # include <cstdio> 2 # include <iostream> 3 # include <cstring> 4 # include <string> 5 # include <cmath> 6 # include <algorithm> 7 # define R register int 8 9 using namespace std; 10 11 const int maxn=100005; 12 int n,k,p1,p2,equ; 13 int cnt,a[maxn],ans; 14 15 int ab (int x) 16 { 17 if(x<0) return -x; 18 return x; 19 } 20 21 int ask (int x,int y) 22 { 23 if(x<0&&y<0) return -min(x,y); 24 if(x>0&&y>0) return max(x,y); 25 return ab(x)+ab(y)+min(ab(x),ab(y)); 26 } 27 28 int main() 29 { 30 scanf("%d%d",&n,&k); 31 for (R i=1;i<=n;++i) 32 { 33 scanf("%d",&a[i]); 34 if(a[i]<0) equ=i; 35 } 36 p1=p2=equ; 37 if(p2<k) p2=p1=k; 38 cnt=1; 39 while (cnt<k&&p1>1) p1--,cnt++; 40 if(cnt==k) ans=ask(a[p1],a[p2]); 41 while (p2<n) 42 { 43 p1++; 44 p2++; 45 ans=min(ask(a[p1],a[p2]),ans); 46 } 47 printf("%d",ans); 48 return 0;
107 D:https://arc101.contest.atcoder.jp/tasks/arc101_b
题意概述:给定一个长度为$n$的序列,对它的所有连续子区间拿出来求中位数(如果是偶数长度就取较大的那个),再将这些中位数排起来求中位数,问这个中位数是几.$n<=10^5$
中位数真是一个神奇的东西...
看到中位数就可以试一下二分答案加上离散化,这道题也不例外.首先可以算出来最终的序列一共有多少个数,所以有多少个数字应该比中位数小也就可以算出来了.二分答案之后怎么做呢?可以首先根据与$ans$的大小关系把序列离散化成$-1,1$,求一个前缀和,现在的问题是怎么求出有多少个连续序列的中位数小于$ans$,容易发现如果$s_i-s_{j-1}>0$,即$s_i>s_{j-1}$,那么$s_j-s_i$就是一个合法的区间了,使用树状数组统计即可.因为树状数组的下标不能为负数,所以统一加上一个$n+1$平移到正数.注意因为二分时可能分到不属于这个序列的数,所以二分时应当取$min$.时间复杂度$O(Nlog^2N)$
1 # include <cstdio> 2 # include <cstring> 3 # include <iostream> 4 # define R register int 5 # define ll long long 6 7 using namespace std; 8 9 const int maxn=100005; 10 int n; 11 int a[maxn],m,l,r,mid,ans,c[maxn],t[maxn<<1]; 12 13 void ins (int x) { for (R i=x;i<=2*n+2;i+=(i&(-i))) t[i]++; } 14 int ask (int x) { int ans=0; for (R i=x;i;i-=(i&(-i))) ans+=t[i]; return ans; } 15 16 bool check (int ans) 17 { 18 ll cnt=0; 19 memset(t,0,sizeof(t)); 20 for (R i=1;i<=n;++i) 21 if(a[i]<=ans) c[i]=c[i-1]+1; 22 else c[i]=c[i-1]-1; 23 ins(n+1); 24 for (R i=1;i<=n;++i) 25 { 26 cnt+=ask(c[i]+n); 27 ins(c[i]+n+1); 28 } 29 ll b=1LL*n*(n+1)/2; 30 b/=2; 31 return cnt>=b+1; 32 } 33 34 int main() 35 { 36 scanf("%d",&n); 37 for (R i=1;i<=n;++i) 38 scanf("%d",&a[i]),m=max(m,a[i]); 39 l=0,r=m,ans=m; 40 while(l<=r) 41 { 42 mid=(l+r)>>1; 43 if(check(mid)) 44 ans=min(ans,mid),r=mid-1; 45 else 46 l=mid+1; 47 } 48 printf("%d",ans); 49 return 0; 50 }
109 C：https://abc109.contest.atcoder.jp/tasks/abc109_c
题意概述：给出$n$个城市的坐标以及起始坐标,每次可以移动$d$的距离（或正或负），求最大的$d$使得可以从起点出发经过所有点。
首先$d$必然是任意点到初始点的距离的因子,否则到不了,然后...答案就是所有点离初始点距离的最大公约数了。(这道$C$是挺水的)
1 # include <cstdio> 2 # include <iostream> 3 4 using namespace std; 5 6 int n,x,ans; 7 int a[100005]; 8 9 int gcd (int a,int b) 10 { 11 return b?gcd(b,a%b):a; 12 } 13 14 int main() 15 { 16 scanf("%d%d",&n,&x); 17 for (int i=1;i<=n;++i) 18 { 19 scanf("%d",&a[i]); 20 a[i]-=x; 21 if(a[i]<0) a[i]=-a[i]; 22 } 23 ans=a[1]; 24 for (int i=2;i<=n;++i) 25 ans=gcd(ans,a[i]); 26 printf("%d",ans); 27 return 0; 28 }
109 D:https://abc109.contest.atcoder.jp/tasks/abc109_d
题意概述：给定一个$H*W$的矩阵，每个格子中有一些硬币，每次可以从一个格子向与他四联通的任意一个格子移动一个硬币(只能移动一次)，使得最终的局面中有偶数个硬币的格子最多。
我感觉这个题挺好的...因为只能移动一次,可以将每个格子的硬币数量$\%2$表示奇偶性,现在问题是怎样移动使得为$0$的格子最多.四联通是一个比较麻烦的东西,因为可能导致处理完的格子再变化.但是可以发现只需要考虑向下和向右移动的情况.偶+偶=偶,奇+奇=偶,偶+奇=奇.从偶数格子里往外移动一定是不合算的,不过两个奇数格子互相传递,中间经过一个偶数格子还是可以的.从一个奇数格子往偶数格子移动一定不会使解更劣:如果后来又将这个多余的硬币移到了奇数格子就可以使解加一,即使不可以,也不过是调换了这两个格子的奇偶性.这是这道题的解法就基本成型了:从上到下遍历棋盘，如果遇到奇数就往右边移动一个，如果已经是最后一列就考虑往下移动，这样构造出来的答案一定是最优的。
1 # include <cstdio> 2 # include <iostream> 3 4 using namespace std; 5 6 int h,w,f=1,j,x,lasx,lasy; 7 int g[502][502]; 8 int H,a[250005],b[250005],c[250005],d[250005]; 9 10 int main() 11 { 12 scanf("%d%d",&h,&w); 13 for (int i=1;i<=h;++i) 14 for (int j=1;j<=w;++j) 15 scanf("%d",&g[i][j]),g[i][j]%=2; 16 for (int i=1;i<=h;++i) 17 { 18 if(i%2) j=1; else j=w; 19 while (1<=j&&j<=w) 20 { 21 if(x) 22 { 23 a[++H]=lasx,b[H]=lasy,c[H]=i,d[H]=j; 24 if(g[i][j]) x=0; 25 } 26 else if(g[i][j]) x=1; 27 lasx=i; 28 lasy=j; 29 j+=f; 30 } 31 f=-f; 32 } 33 printf("%d\n",H); 34 for (int i=1;i<=H;++i) 35 printf("%d %d %d %d\n",a[i],b[i],c[i],d[i]); 36 return 0; 37 }
110 C:https://abc110.contest.atcoder.jp/tasks/abc110_c
题意概述：给出两个字符串，要求构造一个字符的映射关系，使得$S$=$T$.
这题的难度一半在于读题，题面里那个置换非常有迷惑性,题目说的是可以任意将两种字母进行整体交换,其实就是要找到一个映射关系.
似乎就是$NOIP$潜伏者那道题?总之就是开两个桶记录两个集合之间的对应关系,如果出现矛盾就...有矛盾啊.
1 # include <cstdio> 2 # include <iostream> 3 # include <queue> 4 # include <algorithm> 5 # include <cstring> 6 # include <string> 7 # define R register int 8 # define ll long long 9 10 using namespace std; 11 12 const int maxn=200005; 13 char s[maxn],t[maxn]; 14 int a[30],b[30],len,ans=0; 15 16 int main() 17 { 18 scanf("%s",s+1); 19 scanf("%s",t+1); 20 len=strlen(s+1); 21 for (R i=1;i<=len;++i) 22 { 23 if(a[ s[i]-'a'+1 ]==0) a[ s[i]-'a'+1 ]=t[i]-'a'+1; 24 else if(a[ s[i]-'a'+1 ]!=t[i]-'a'+1) ans=1; 25 if(b[ t[i]-'a'+1 ]==0) b[ t[i]-'a'+1 ]=s[i]-'a'+1; 26 else if(b[ t[i]-'a'+1 ]!=s[i]-'a'+1) ans=1; 27 } 28 if(ans) printf("No"); 29 else printf("Yes"); 30 return 0; 31 }
111 D:https://abc110.contest.atcoder.jp/tasks/abc110_d
题意概述:给出$n,m$,问满足下面式子的数列有多少个.$n<=10^5,m<=10^9$
$10^9$基本上就钦定是根号级别的算法了,线性是不可能的,$log$还不如开到$long$ $long$,大胆猜测首先要分解$m$。
题目就被转化成了,将$m$的质因数排成一排,可以任意合并,求个数小于等于$n$的方案数.虽然题目说的是严格等于$n$,然而不够时添$1$即可.
然后就会发现这个转化毫无意义~,而且还提高了编写的难度.
可以用另一种思路来做这个题,考虑设置$n$个盒子,每个盒子里存一个数,一开始都是$1$,每次往盒子里面放数就相当于把这个数乘进去.此时我们有了一个非常有利于解题的模型,对于每种质因数单独考虑,利用乘法原理统计答案.对于每种质因数应当怎么做呢?等价于将$n$个无差别物体放进$m$个有差别盒子中,允许有空盒子,设置$m$个虚点后插板法即可.
1 # include <cstdio> 2 # include <iostream> 3 # include <queue> 4 # include <cstring> 5 # include <string> 6 # define mod 1000000007 7 # define R register int 8 # define ll long long 9 10 using namespace std; 11 12 const int maxn=100100; 13 int n,m,fac=1,h; 14 int p[maxn],c[maxn]; 15 ll f[maxn],inv[maxn]; 16 ll ans=1; 17 18 ll qui (ll a,ll b) 19 { 20 ll s=1; 21 while (b) 22 { 23 if(b&1LL) s=s*a%mod; 24 a=a*a%mod; 25 b=b>>1LL; 26 } 27 return s%mod; 28 } 29 30 ll C (int a,int b) 31 { 32 return f[b]*qui(f[a],mod-2)%mod*qui(f[b-a],mod-2)%mod; 33 } 34 35 int main() 36 { 37 scanf("%d%d",&n,&m); 38 f[0]=1; 39 for (R i=1;i<=n+50;++i) 40 f[i]=1LL*i*f[i-1]%mod; 41 for (R i=2;1LL*i*i<=m;++i) 42 { 43 if(m%i==0) p[++h]=i; 44 while (m%i==0) c[h]++,m/=i; 45 } 46 if(m!=1) p[++h]=m,c[h]=1; 47 for (R i=1;i<=h;++i) 48 ans=(ans*C(n-1,n+c[i]-1))%mod; 49 printf("%lld",ans%mod); 50 return 0; 51 }
111 C:https://abc111.contest.atcoder.jp/tasks/arc103_a
2018.10.6：一个月没做发现只缺了两场,开心~
题意概述：给定一个长度为$n$的序列(保证$n$是偶数),可以随意改动序列中的数字,要求改成ABABABAB的形式,($A!=B$),求最少需要改动的次数.
自以为是道很水的题,然而$WA$了两次,要注意细节啊...其实思路挺简单的,对于奇数部分和偶数部分分开考虑,找到出现次数最多的那个数字,之后将不是这个数字的数字都变成这个就好了...吗?注意$A!=B$,这时我的做法就处理不了了.其实加上这个限制也没有麻烦多少,再记录两个次大值,如果出现冲突就找一个次大值换上,这道题就做完了.最近见了不少记录次大值的题目,过两天凑上几个一起写一下.
1 # include <cstdio> 2 # define R register int 3 4 using namespace std; 5 6 const int maxn=100005; 7 int n,M,ans,cnt,a[maxn],v[maxn],b1,b2,mb1,mb2,mc1,mc2,c1,c2; 8 9 inline int min (int a,int b) { return b>=a?a:b; } 10 inline int max (int a,int b) { return b>=a?b:a; } 11 12 int main() 13 { 14 scanf("%d",&n); 15 for (R i=1;i<=n;++i) 16 scanf("%d",&a[i]),M=max(M,a[i]); 17 for (R i=1;i<=n;i+=2) 18 v[ a[i] ]++; 19 for (R i=1;i<=M;++i) 20 if(v[i]>=mb1) 21 b2=b1,mb2=mb1,mb1=v[i],b1=i; 22 else 23 if(v[i]>mb2) mb2=v[i],b2=i; 24 for (R i=1;i<=n;i+=2) 25 v[ a[i] ]--; 26 for (R i=2;i<=n;i+=2) 27 v[ a[i] ]++; 28 for (R i=1;i<=M;++i) 29 if(v[i]>=mc1) 30 c2=c1,mc2=mc1,mc1=v[i],c1=i; 31 else 32 if(v[i]>mc2) mc2=v[i],c2=i; 33 if(b1!=c1) 34 ans=n-mb1-mc1; 35 else 36 ans=min(n-mb1-mc2,n-mb2-mc1); 37 printf("%d",ans); 38 return 0; 39 }
111 D:https://abc111.contest.atcoder.jp/tasks/arc103_b
题意概述：构造题。给定平面上的$n$个点，要求构造一条包含多个部分的机械臂,每个部分之间有一个关节,起点在$(0,0)$位置，每个关节处可以选择四种方向，要求终点处正好在给定点上。需要构造的部分：部分的数量，每段的长度，对于每一个点应采取的方向序列。无解输出$-1$
一看到这道题以为自己又可以开心的水一道水题了,结果是被水题给水了...后来发现这次是和$ARC$一起办的,所以难一些.首先我觉得每个部分如果都是$1$肯定是最好的,非常灵活,部分数量就是所有点中距离原点的曼哈顿距离最大的那个距离.如果有两个点的曼哈顿距离的奇偶性不同,那么这两个点就不可能被同时到达,为无解.这个做法听起来非常美好,$OLE$之后才发现题目上写了一行$m<=40$,也就是不能超过四十个部分...以后再也不用全屏翻译的功能了,把这句话都给翻译成乱码了.看了题解发现是一个二进制构造,其实这种拼凑数字的题好多都是二进制,但是这个题包装的比较复杂,就没看出来.判断无解那里我写的是正确的，的确要求奇偶性相同，构造方法如下：从$2^{31}$开始,把每个$2$的幂次方都放上,如果目标距离是偶数就再加一条长度为一的边.摆放的时候采取一种修修补补的策略,目前哪个方向离目标差的比较多就把下一条机械臂指向那个方向即可.至于证明,有一个博客写的非常好:https://blog.csdn.net/zxyoi_dreamer/article/details/82905515
1 # include <cstdio> 2 # include <iostream> 3 # include <cstring> 4 # include <cmath> 5 # include <algorithm> 6 # include <string> 7 # define ll long long 8 # define R register int 9 10 using namespace std; 11 12 const int maxn=1003; 13 int n,m,f=-1; 14 ll x[maxn],y[maxn]; 15 ll d[35],len; 16 17 inline ll ab (ll x) { if(x<0) return -x; return x; } 18 19 inline void solve (ll x,ll y) 20 { 21 for (R i=1;i<=m-((f&1)^1);++i) 22 { 23 if(ab(x)>ab(y)) 24 { 25 if(x>0) putchar('R'),x-=d[i]; 26 else putchar('L'),x+=d[i]; 27 } 28 else 29 { 30 if(y>0) putchar('U'),y-=d[i]; 31 else putchar('D'),y+=d[i]; 32 } 33 } 34 if(f%2==0) putchar('L'); 35 putchar('\n'); 36 } 37 38 int main() 39 { 40 scanf("%d",&n); 41 for (R i=1;i<=n;++i) 42 { 43 scanf("%lld%lld",&x[i],&y[i]); 44 len=ab(x[i])+ab(y[i]); 45 if(f==-1) f=len%2; 46 if(f%2!=len%2) 47 { 48 printf("-1"); 49 return 0; 50 } 51 } 52 for (R i=30;i>=0;--i) d[++m]=(1LL<<i); 53 if(f%2==0) d[++m]=1; 54 printf("%d\n",m); 55 for (R i=1;i<=m;++i) 56 { 57 printf("%lld",d[i]); 58 if(i==m) printf("\n"); 59 else printf(" "); 60 } 61 for (R i=1;i<=n;++i) 62 solve(x[i]+((f&1)^1),y[i]); 63 return 0; 64 }
112 C:https://abc112.contest.atcoder.jp/tasks/abc112_c
本来今天是要做生成树的题目的，但是重启电脑后发现用来列任务清单的软件没了...所以又来做$AT$了.
这场$AT$好像挺水的,竟然有$576$个$AK$选手,$\%\%\%$
题意概述:有一个金字塔,中心点位于$(C_x,C_y)$每个点的高度有这样一个计算公式:$max(H-|x-C_x|-|y-C_y|,0)$,给出$n$个点的坐标及其高度，试求出金字塔的中心点。
看到这道题觉得简直太妙啦，是不是要用到绝对值化简的性质呢？是不是可以二分套二分呢？这时我看到了数据范围：保证中心点的横纵坐标都是整数，且位于区间$[0,100]$,给出点的数量不会超过$100$...突然想起金牌爷之前跟我说的"做思维题做多了,每见到一个题就有五六种奇思妙想,其实暴力就行"...
注意...边界...有一个小细节:如果某一次的高度为$0$,其实也就是给金字塔的高度设定了一个上限,这里一定要记住保存一下这个边界,等到输出前再看一看高度是否超出了这个边界即可$AC$
1 # include <cstdio> 2 # include <iostream> 3 # define R register int 4 5 using namespace std; 6 7 const int maxn=102; 8 int n,w; 9 int x[maxn],y[maxn],h[maxn]; 10 11 inline int ab (int a) { if(a<0) return -a; return a; } 12 13 inline bool check (int cx,int cy) 14 { 15 int hh,X=2000000005; 16 w=-1; 17 for (R i=1;i<=n;++i) 18 { 19 if(h[i]==0) 20 { 21 if(w==-1) 22 { 23 X=min(X,ab(x[i]-cx)+ab(y[i]-cy)); 24 continue; 25 } 26 if(w-ab(x[i]-cx)-ab(y[i]-cy)>0) return false; 27 continue; 28 } 29 hh=ab(x[i]-cx)+ab(y[i]-cy)+h[i]; 30 if(w==-1) w=hh; 31 if(hh!=w) return false; 32 } 33 if(w==-1) return false; 34 if(w>X) return false; 35 return true; 36 } 37 38 int main() 39 { 40 scanf("%d",&n); 41 for (R i=1;i<=n;++i) 42 scanf("%d%d%d",&x[i],&y[i],&h[i]); 43 for (R i=0;i<=100;++i) 44 for (R j=0;j<=100;++j) 45 { 46 if(check(i,j)) 47 { 48 printf("%d %d %d",i,j,w); 49 return 0; 50 } 51 } 52 return 0; 53 }
112 D:https://abc112.contest.atcoder.jp/tasks/abc112_d
题意概述：把一个数分成$n$份,使这$n$个数的最大公约数最大,输出.
$AT$还真喜欢出构造题.然而我最不会的就是构造题了.
这道题几乎可以说是我自己做出来的第一道构造啦,留念$2018-10-08$ $21:38:31$
首先我们可以对这个式子进行一番变形:$k*a_1+k*a_2+……+k*a_n=M$，再提取一下公因数：
所以答案一定是$m$的一个因子,但是不是每一个因子都能成为答案,玩一下第一个样例就能发现.为什么呢？因为$m$除以$k$之后可能根本组不出$n$个数来...所以要判断一下除出来的数字是否大于$n$.这里不能再枚举质因数了,要枚举约数,两个注意事项:
·$1$和$n$都是$n$的约数;
·每枚举一个因数就要枚举$m$除以它的那个因数,可以将复杂度降到$\sqrt{M}$
1 # include <cstdio> 2 # include <iostream> 3 using namespace std; 4 5 int n,m,ans; 6 7 inline void check (int x) { if(m/x>=n) ans=max(ans,x); } 8 9 int main() 10 { 11 scanf("%d%d",&n,&m); 12 for (int i=1;i*i<=m;++i) 13 { 14 if(m%i) continue; 15 check(i); 16 check(m/i); 17 } 18 printf("%d",ans); 19 return 0; 20 }
