Codeforces Round #566 (Div. 2)

最新推荐文章于 2019-06-17 15:11:25 发布
最新推荐文章于 2019-06-17 15:11:25 发布
阅读量238 收藏
点赞数
分类专栏： Codeforces
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/syh0313/article/details/91639446
版权

A.Filling Shapes

  数学题

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <string>
#include <cmath>
#include <algorithm>
using namespace std;
int n;
long long p[40];
int main()
{
	scanf("%d",&n);
	p[0]=1; for (int i=1;i<=30;i++) p[i]=p[i-1]*2ll;
	if (n&1) return printf("0\n"),0;
	printf("%lld\n",p[n/2]);
return 0; 
}

B.Plus from Picture

  找出最高最低最左最右的点，然后暴力判断就好了

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <set>
#include <cmath>
#define pii pair<int,int>;
using namespace std;
const int maxn=1010;
int n,m,topt,xmin,xmax,ymin,ymax,a,b,c,d,e,f,g,h,sum,now;
char s[maxn];
bool ff[510][510],ok;
set<int>ss;
int main()
{
	scanf("%d%d",&n,&m); xmin=501; ymin=501;
	for (int i=1;i<=n;i++)
	{
		scanf("%s",s+1);
		for (int j=1;j<=m;j++)
		 if (s[j]=='*') 
		 {
		 	if (i>xmax) {xmax=i; g=i; h=j;}
		 	if (i<xmin) {xmin=i; e=i; f=j;}
		 	if (j>ymax) {ymax=j; a=i; b=j;}
		 	if (j<ymin) {ymin=j; c=i; d=j;}
		 	ff[i][j]=1;
		 	sum++;
		 }
	}
	if (a!=c || h!=f) return printf("NO\n"),0;
	for (int i=e;i<=g;i++) if (!ff[i][f]) return printf("NO\n"),0;else now++;
	for (int j=d;j<=b;j++) if (!ff[a][j]) return printf("NO\n"),0;else now++;
	now--;
	if (a==e || a==g) return printf("NO\n"),0;
	if (sum>now) return printf("NO\n"),0;
	printf("YES\n");	
return 0;
}

C.Beautiful Lyrics

   模拟排序，首先将这些字符串按含有元音数量从大到小排序，元音数量相等时按a,e,i,o,u排序，

   那么我们应该先找出最多的元音数量相同且后缀元音一样的对数S，然后再将剩下的元音数量相同的字符对分别与S配对，

   配完之后S中多余的再两两配对

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <set>
#include <cmath>
#define pii pair<int,int>;
using namespace std;
const int maxn=1010;
int n,m,topt,xmin,xmax,ymin,ymax,a,b,c,d,e,f,g,h,sum,now;
char s[maxn];
bool ff[510][510],ok;
set<int>ss;
int main()
{
	scanf("%d%d",&n,&m); xmin=501; ymin=501;
	for (int i=1;i<=n;i++)
	{
		scanf("%s",s+1);
		for (int j=1;j<=m;j++)
		 if (s[j]=='*') 
		 {
		 	if (i>xmax) {xmax=i; g=i; h=j;}
		 	if (i<xmin) {xmin=i; e=i; f=j;}
		 	if (j>ymax) {ymax=j; a=i; b=j;}
		 	if (j<ymin) {ymin=j; c=i; d=j;}
		 	ff[i][j]=1;
		 	sum++;
		 }
	}
	if (a!=c || h!=f) return printf("NO\n"),0;
	for (int i=e;i<=g;i++) if (!ff[i][f]) return printf("NO\n"),0;else now++;
	for (int j=d;j<=b;j++) if (!ff[a][j]) return printf("NO\n"),0;else now++;
	now--;
	if (a==e || a==g) return printf("NO\n"),0;
	if (sum>now) return printf("NO\n"),0;
	printf("YES\n");	
	/*sort(a+1,a+topt+1,cmp); 
	for (int i=1;i<=topt-1;i++)
	{
		if (a[i].x!=a[i-1].x && a[i].x!=a[i+1].x) yy=y;
	 	if (a[i].x==a[i+1].x) ss.insert(a[i].x);
	 	if (a[i+1].x-a[i].x>1) {f=1; break;}
	}
	if ((int)ss.size()>=2 || f) return printf("NO\n"),0;
	sort(a+1,a+topt+1,cmp1); ss.clear();
	for (int i=1;i<=topt-1;i++)
	{
		if (a[i].y==a[i+1].y || f) ss.insert(a[i].y);
		if (a[i+1].y-a[i].y>1) {f=1; break;}
	}
	if ((int)ss.size()>=2 || f) return printf("NO\n"),0;
	printf("YES\n");*/
return 0;
}

D.Complete Mirror

    1.首先找出直径，先判断直径的两端点是否可行(BFS判断)，若可行，直接输出。

    2.再判断直径上有多少个点，若有偶数个点，直接不行(直径两端点始终不能满足条件)，若有奇数个点，我们构造一棵以直径中点为root的树

    3.对于这棵树，考虑有now个深度<最大深度的叶子节点

       (1) 判断root可不可以，若可以，直接输出

       (2) 若now>=2 直接不行(任一now<maxdep的叶子节点和直径一端不能同时满足)

       (3) 若now==1，那么满足条件的根节点只可能是该节点，判断下改点是否可以就ok了

       (4) 若now==0，找一条从根到叶子节点的所有father都只有一个儿子，若无这条链，直接不行，若有则判断这条链的叶子节           点是否满足条件

zzytql  orz

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <string>
#include <algorithm>
using namespace std;
const int maxn=100010;
int n,st[maxn<<1],nt[maxn<<1],to[maxn<<1],topt,a[2],di,d[2][maxn],deg[maxn],dep[maxn],root;
int dd[maxn],q[maxn],head,tail,col[maxn],ma,kk[maxn];
int madep;
bool f[maxn];
void add(int x,int y)
{to[++topt]=y; nt[topt]=st[x]; st[x]=topt;}
void find(int x,int de,int id)
{
	f[x]=1; int p=st[x]; if (de>ma) {ma=de; a[id]=x;}
	while (p)
	{
		if (!f[to[p]]) find(to[p],de+1,id);
		p=nt[p];
	}
}
void calc(int x,int y,int num)
{
	if (x==y) {di=num; return;}
	f[x]=1; int p=st[x];
	while (p)
	{
		if (!f[to[p]]) calc(to[p],y,num+1);
		p=nt[p];
	}
}
void dfs(int x,int num,int id)
{
	f[x]=1; d[id][x]=num; int p=st[x];
	while (p)
	{
		if (!f[to[p]]) dfs(to[p],num+1,id);
		p=nt[p];
	}
}
void getdep(int x,int de)
{
	f[x]=1; dep[x]=de; int p=st[x];
	if (x==a[0]) madep=de;
	if (deg[x]==1) kk[x]=x;
	int nn=0,too;
	while (p)
	{
		if (!f[to[p]]) getdep(to[p],de+1),nn++,too=to[p];
		p=nt[p];
	}
	if (nn==1) kk[x]=kk[too];
}
bool ok(int x)
{
	for (int i=1;i<=n;i++) dd[i]=0,col[i]=0;
	q[1]=x; head=tail=1; dd[x]=1; col[1]=deg[x];
	while (head<=tail)
	{
		int x=q[head++]; 
		if (!col[dd[x]]) col[dd[x]]=deg[x];
		else if (col[dd[x]]!=deg[x]) return 0;
		int p=st[x];
		while (p)
		{
			if (!dd[to[p]]) {dd[to[p]]=dd[x]+1; q[++tail]=to[p];}
			p=nt[p];
		}
	}
return 1;
}
int main()
{	
	scanf("%d",&n); a[0]=a[1]=1;
	for (int i=1;i<n;i++)
	{
		int xx,yy; scanf("%d%d",&xx,&yy);
		add(xx,yy); add(yy,xx); deg[xx]++; deg[yy]++;
	}
	ma=0; find(1,0,0); memset(f,0,sizeof f); ma=0; find(a[0],0,1);
	memset(f,0,sizeof f); calc(a[0],a[1],1); 
	if (ok(a[0])) return printf("%d\n",a[0]),0;
	if (ok(a[1])) return printf("%d\n",a[1]),0;
	if (di%2==0) return printf("-1\n"),0;
	memset(f,0,sizeof f); dfs(a[0],0,0); memset(f,0,sizeof f); dfs(a[1],0,1);
	for (int i=1;i<=n;i++) if (d[0][i]==di/2 && d[1][i]==di/2) {root=i; break;}
	//printf("lala %d %d %d %d\n",a[0],a[1],di,root);
	memset(f,0,sizeof f); getdep(root,0); int now=0,lc;
	for (int i=1;i<=n;i++)
	 if (i!=a[0] && i!=a[1] && deg[i]==1 && dep[i]<madep) now++,lc=i;
	if (now>=2) return printf("-1\n"),0;
	if (now==1)
	{
		if (ok(lc)) printf("%d\n",lc);else printf("-1\n");
		return 0;
	}
	if (!now)
	{
		if (ok(root)) return printf("%d\n",root),0;
		if (!kk[root]) return printf("-1\n"),0;
		if (ok(kk[root])) return printf("%d\n",kk[root]),0;
		printf("-1\n"); return 0;
	}
return 0;
}

 E.Product Oriented Recurrence

   欧拉降幂+矩阵快速幂

   首先等式两边同时乘上c^n，化简可得：c^n*f(n)=c^(n-1)*f(n-1) * c^(n-2)*f(n-2) * c^(n-3)*f(n-3)

   设g(n)=c^n*f(n);     那么问题就转化为求g(n);

   不妨设g(n)=g(1)^x(n)*g(2)^y(n)*g(3)^z(n);

    所以x(n)=x(n-1)+x(n-2)+x(n-3) 这个就可以矩阵快速幂一波(注意此时矩阵是在幂指数上搞的，所以要用欧拉降幂就是说%(mo-1))，然后y,z同理，那么g(n)就可以搞出来，然后f(n)就可以搞出来

    qwq我的代码比较窒息，被大佬教育了

    

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <string>
using namespace std;
long long n,f1,f2,f3,c,g1,g2,g3,mo,x1,x2,x3,ans=1,mi;
const int N=3;  
long long tmp[N][N],now[N][N]; 
void multi(long long a[][N],long long b[][N],int n)  
{  
    memset(tmp,0,sizeof tmp);  
    for(int i=0;i<n;i++)  
     for(int j=0;j<n;j++)  
      for(int k=0;k<n;k++)  
        tmp[i][j]=(tmp[i][j]+a[i][k]*b[k][j]%mo)%mo;  
    for(int i=0;i<n;i++)  
     for(int j=0;j<n;j++)  
        a[i][j]=tmp[i][j];  
}  
long long res[N][N];  
void Pow(long long a[][N],long long n)  
{  
    memset(res,0,sizeof res);  
    for(int i=0;i<3;i++) res[i][i]=1;  
    while(n)  
    {  
        if(n&1)  
            multi(res,a,3);  
        multi(a,a,3); 
        n>>=1;  
    }  
}
long long po(long long a,long long b)
{
	if (b==0) return 1;
	if (b==1) return a%mo;
	long long c=po(a,b/2);
	if (b&1) return c*c%mo*a%mo;else return c*c%mo;
}
int main()
{
	scanf("%lld%lld%lld%lld%lld",&n,&f1,&f2,&f3,&c); mo=1e9+7;
	g1=f1*c%mo; g2=f2*c%mo*c%mo; g3=f3*c%mo*c%mo*c%mo;
	mo=1e9+6; x1=0; x2=0; x3=1;
	now[0][0]=1; now[0][1]=1; now[0][2]=0;
	now[1][0]=1; now[1][1]=0; now[1][2]=1;
	now[2][0]=1; now[2][1]=0; now[2][2]=0;
	Pow(now,n-3);
	mi=((x1*res[0][0]%mo+x2*res[1][0]%mo)%mo+x3*res[2][0]%mo)%mo;
	mo=1e9+7; ans=po(g1,mi);

	mo=1e9+6; x1=0; x2=1; x3=0;
	now[0][0]=1; now[0][1]=1; now[0][2]=0;
	now[1][0]=1; now[1][1]=0; now[1][2]=1;
	now[2][0]=1; now[2][1]=0; now[2][2]=0;
	Pow(now,n-3);
	mi=((x1*res[0][0]%mo+x2*res[1][0]%mo)%mo+x3*res[2][0]%mo)%mo;
	mo=1e9+7; ans=ans*po(g2,mi)%mo;

	mo=1e9+6; x1=1; x2=0; x3=0;
	now[0][0]=1; now[0][1]=1; now[0][2]=0;
	now[1][0]=1; now[1][1]=0; now[1][2]=1;
	now[2][0]=1; now[2][1]=0; now[2][2]=0;
	Pow(now,n-3);
	mi=((x1*res[0][0]%mo+x2*res[1][0]%mo)%mo+x3*res[2][0]%mo)%mo;
	mo=1e9+7; ans=ans*po(g3,mi)%mo;

	long long inv=po(po(c%mo,n),mo-2);
	printf("%lld\n",ans*inv%mo);
return 0;	
}

F.Maximum Sine

   神仙数学题,设在x=k处取得最大值，二分要求的x到k之间的最小距离，可用类欧判断是否可行，然后知道了最小距离dis后用exgcd反解出x   

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <algorithm>
using namespace std;
int T;
long long a,b,p,q;
long long S(long long k){return (k*(k+1)/2ll);}
long long f(long long a,long long b,long long c,long long n)
{
    if (!a) return 0;
    if (n<0) return 0;
    if (!n) return b/c;
    if(a>=c || b>=c) return ((a/c)*S(n)+(n+1)*(b/c)+f(a%c,b%c,c,n));
    long long m=(a*n+b)/c; return (m*n-f(c,c-b-1,a,m-1));
}
long long work(long long a,long long b,long long c,long long l,long long r)
{return f(a,b,c,r)-f(a,b,c,l-1);}
void exgcd(long long a,long long b,long long &x,long long &y) 
{if(!b) {x=1; y=0; return;} exgcd(b,a%b,y,x),y-=a/b*x;}
long long solve(long long p,long long q,long long t) 
{
  long long gd=__gcd(p,q); if(t%gd!=0) {return 1e18;}
  p/=gd; q/=gd; t/=gd; long long x,y; exgcd(p,q,x,y);
  x*=t; y*=t; long long k=(a-x)/q; x+=k*q;
  while(x>=a) {x-=q;} while(x<a) {x+=q;}
return x;
}
int main()
{
	scanf("%d",&T);
	while (T--)
	{
		scanf("%lld%lld%lld%lld",&a,&b,&p,&q); p<<=1; q<<=1;
		long long l=0,m=q/2,r=m,ans=r;
		while (l<=r)
		{
			long long mid=(l+r)>>1,ll=m-mid,rr=m+mid;
			if (work(p,q-ll,q,a,b)-work(p,q-rr-1,q,a,b)) {ans=min(ans,mid); r=mid-1;}else l=mid+1;
		}
		printf("%lld\n",min(solve(p,q,m-ans),solve(p,q,m+ans)));
	}
return 0;
}

 

syh0313
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
Codeforces Round #563 (Div. 2)
On the way
06-05 146
Codeforces Round #563 (Div. 2) D: 题意是给定 nnn 和 xxx ，构建这样最长的数组 aaa 满足：1≤ai≤2n1\leq a_i \leq 2^n1≤ai​≤2n 且任意一个连续区间的异或和不等于 000 或 xxx 。 一个区间[l, r]的异或和可以由l-1的异或前缀和和r的异或前缀和异或起来表示。 连续区间异或和不为0则说明构造出的数组中，没有任何一...
Codeforces Round #566 (Div. 2) A-E
过路人丶的博客
06-13 339
题目链接:http://codeforces.com/contest/1182 A. Filling Shapes 每两列有2种方法，所以是2的n/2次方 #include <bits/stdc++.h> using namespace std; typedef long long ll; const int mod = 998244353; const int mx =...
Codeforces Round #566 (Div. 2) C
mzip的博客
06-17 185
链接: https://codeforces.com/problemset/problem/1182/C 题意: A lyric is beautiful if and only if it satisfies all conditions below. The number of vowels in the first word of the first line is the same as ...
Codeforces Round #566 (Div. 2) E Product Oriented Recurrence（矩阵快速幂 + 欧拉降幂）
Wind_bow的博客
06-12 394
E. Product Oriented Recurrence time limit per test 1 second memory limit per test 256 megabytes input standard input output standard output Let fx=c2x−6⋅fx−1⋅fx−2⋅fx−3fx=c2x−6⋅fx−1⋅fx−2⋅fx−3 ...
Codeforces #566 div2
weixin_44203780的博客
06-12 157
https://codeforces.com/contest/1182 A.找规律题，偶数输出2^(n/2)，奇数0； #include<bits/stdc++.h> using namespace std; #define ll long long int main() { int n; scanf("%I64d",&n); if(n%2==0) { ...
Codeforces Round #723 (Div. 2).md
10-13
Codeforces Round #723 (Div. 2).md
Codeforces Round #627 (Div. 3) B. Yet Another Palindrome Problem
01-03
就是把所有相等的数放到一个vector里，如果他出现大于2次，看最远的间距是否大于2即可，找到一个就可以 代码： #include #include #include #include #include #include #include #include #include #include #...
Codeforces Round #627 (Div. 3) C. Frog Jumps（思维）
01-20
传送门 题意： 开始位置在0，问能否跳到n+1位置 每步只能跳d 在1——n每个位置有方向，L,R，求d的最小值 思路： 只用找相邻两个R之间的最大值即可 代码： #include #include ...typedef long long l
Codeforces Round #479 (Div. 3) E. Cyclic Components
01-03
E. Cyclic Components 题目链接-E. Cyclic Components 题目大意 给你nnn个点和mmm条边，求所构成图中单圈环的个数 ...并查集并查集并查集 很明显单圈环每个点的度都为222，所以我们可以用数组cnt[]记录每个点的度，...
Codeforces Round #628 (Div. 2)
01-03
给两两节点放一个数字（0~n-2 唯一） 给你一棵树，求所有任意两节点相连的路以外的路上的数字的最小值最小 思路 构造 若一个点连了三条边及以上，则这个点的边从最小值开始赋值。其他边从最大点开始赋值。 证明：一...
C. Beautiful Lyrics （模拟）构造
caowenbo2311694605的博客
06-13 331
C. Beautiful Lyrics time limit per test 1 second memory limit per test 256 megabytes input standard input output standard output You are givennnwords, each of which consists of lowercase al...
Codeforces Round #545 (Div. 2)B. Circus（思维）
sugarbliss
03-10 457
题目链接：http://codeforces.com/contest/1138/problem/B 题意：给你两列01串，从第一列中挑出一些0和1，并且下标只能用一次，也就是在第一列挑下标为1的数第二列就不能用了，最终挑选n/2个数，并且使第一列1的个数和第二列相等，输出第一列挑选数的下标，如果不能输出-1。 思路：挑选第一列为1第二列为0的数相当于第一列+1，挑选第一列为0第二列为1的数相当...
【题解】A. Filling Shapes⭐ 【水题】
Suprit's blog
06-12 506
A. Filling Shapes You have a given integer n. Find the number of ways to fill all 3×n tiles with the shape described in the picture below. Upon filling, no empty spaces are allowed. Shapes cannot over...
Product Oriented Recurrence(Codeforces Round #566 (Div. 2)E+矩阵快速幂+欧拉降幂)
Dillonh的博客
06-12 440
传送门 题目 fn=c2∗n−6fn−1fn−2fn−3 \begin{aligned} &amp;f_n=c^{2*n-6}f_{n-1}f_{n-2}f_{n-3}&amp;\\ \end{aligned} ​fn​=c2∗n−6fn−1​fn−2​fn−3​​​ 思路 我们通过迭代发现fnf_nfn​其实就是由cx1,f1x2,f2x3,f3x4c^{x_1},f_1^{x_...
Codeforces Round #566 (Div. 2)B. Plus from Picture
轻舟载君不载愁
06-15 260
You have a given picture with size w×h. Determine if the given picture has a single “+” shape or not. A “+” shape is described below: A “+” shape has one center nonempty cell. There should be some (at...
Arcmap天地图历史影像插件全国版-免费！免费！
最新发布
07-22
2024版“天地图”首次开放多时相影像专题，公众可查看32个省级天地图节点近半个世纪的历史影像；本arcmap插件集成了全国天地图历史影像接口，累计上千期。 免费使用。不收费！！不收费！！不收费！！
codeforces round #772 (div. 2)
03-16
好的，Codeforces Round #772 (Div. 2) 是一场由Codeforces主办的在线编程竞赛。 该比赛的比赛链接为：https://codeforces.com/contest/772。 该比赛分为两个级别：Div.1和Div.2。 Div. 1级别通常比Div. 2级别更...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值