2021暑假集训 Day 1

学校刷题Day 1 -- 2021-8-9

A

Decription:

定义操作a(X)B=((a&b)+a(a|b))>>1
给出n个数ai
每次用操作合并任意相邻ai
求进行n-1次操作以后可能得到的最终结果
所有结果从小到大输出

Input

第一行一个n
第二行给出a1,a2,…,an-1,an

Output

按题目要求作答

Sample Input:

4
1 4 3 2

Sample Output:

1 2

Data Constraint

30% n<=10 0<=a[i]<=7
另外40% n<=150 0<=a[i]<=3
100% n<=150 0<=a[i]<=7

看到这道题,额,想了很久,最后想到正解的区间dp但那时候时间不够了,就没打,赛后做了做A了。
首先先讲一下思路,就是定义一个 f i j k f_{i j k} fijk 表示在此区间内能不能通过某种合并方式来获得这个数,那么怎么转移,我们可以先枚举1到n-1轮,然后因为只能一次吞并一个数,那么操作次数即就是即将枚举的区间的长度,然后我们要枚举的是区间的左驱,用L来表示,那么R就可以通过左驱和区间长度来获得,为了合并区间,我们还要枚举在区间内的一把“刀”将此区间斩成两个区间来合并,最后就是枚举数字来进行合并。

代码

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
int n,a,f[151][151][8];
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%d",&a),f[i][i][a]=1;
	for(int len=1;len<=n;len++)
	for(int l=1;l<=n-len+1;l++)
	for(int mid=l;mid<l+len-1;mid++)
	{
		int r=l+len-1;
		for(int num1=0;num1<8;num1++)
		for(int num2=0;num2<8;num2++)
		{
			f[l][r][(num1+num2)>>1]+=f[l][mid][num1]&f[mid+1][r][num2];
			if(f[l][r][(num1+num2)>>1]>1) f[l][r][(num1+num2)>>1]=1; 
		}
	}
	for(int i=0;i<8;i++) if(f[1][n][i]) printf("%d ",i);
	return 0;
} 

B

Decription:

a数组长度为n,b数组长度为m,c是n行m列的矩阵,且 c i , j = a i ∗ b j 。 c_{i,j}=a_{i}*b_{j}。 ci,j=aibj
定义矩阵的权值为矩阵内所有权值的和
给出限制[L,R]
问权值在这个限制内子矩阵有多少个

Input

第一行依次给出n,m,L,R
第二行n个数表示数组a
第三行m个数表示数组b

Output

输出满足条件的子矩阵个数

Sample Input:

3 3 3 8
3 2 3
2 3 1

Sample Output:

10

Data Constraint

30% 1<=n,m<=100
100% 1<=n,m<=1000,1<=L<=R<=1e18
1<=a[i],b[i]<=1e6

思路:首先赛时想到的就是二维前缀和,但是实际上并不用将c数组给构出来,因为我们推导一下一个公式,对于一个子矩阵来讲当左上角为[x1,y1],右下角为[x2,y2],那么他的权值可以表示为ax1*(by1+…+by2) +…ax2*(by1+…+by2)即(ax1+…+ax2)*(by1+…+by2),那么我们可以知道ax和by之间没有任何必然关系,那么我们就可以分开处理一下,可以分开处理。枚举a的可能区间[x1,x2],预处理所有的by区间和对于一个限制[L,R]若x选定为[x1,x2],则(L/suma[x1,x2],R/suma[x1,x2])为by区间和可以选择的范围(注意考虑小数取整)对于这个范围我就可以在预处理的by区间在排序后进行二分,找到符合的区间个数。----elba

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#define ll long long
#define N 1000
using namespace std;
ll n,m,l,r,a[N+1],b[N+1],sumb[N*N+1],len,l1,r1,ans;
int main()
{
    scanf("%lld%lld%lld%lld",&n,&m,&l,&r);
    for(int i=1;i<=n;i++) scanf("%d",&a[i]);
    for(int i=1;i<=m;i++) scanf("%d",&b[i]);
    for(int i=2;i<=n;i++) a[i]+=a[i-1];
    for(int i=1;i<=m;i++) b[i]+=b[i-1];
    for(int i=1;i<=m;i++)
    for(int j=1;j<=i;j++) sumb[++len]=b[i]-b[j-1];
    sort(sumb+1,sumb+len+1);
    for(int i=1;i<=n;i++)
    for(int j=1;j<=i;j++)
      {
      	l1=l/(a[i]-a[j-1]),r1=r/(a[i]-a[j-1]);
      	if(l%(a[i]-a[j-1])!=0) l1++;
      	int q=lower_bound(sumb+1,sumb+len+1,l1)-sumb;
		int p=upper_bound(sumb+1,sumb+len+1,r1)-sumb;
		ans+=p-q;
	  }
	printf("%lld",ans);
	return 0;
} 

C不会哦,后面更新。

D

Decription:

有n个人依次排队打饭,有m种饭菜可以选择,每个人可能选择其中一种,如果相邻排队的人打的菜一样,那么就会影响彼此吃饭的心情,求这个队伍中有人被影响心情的状态数,对100003取余。

Input

一行整数依次给出m,n

Output

相应的状态数

Sample Input:

2 3

Sample Output:

6

Data Constraint

30%:1<=n<=20,1<=m<=2
100%: 1<=m<=1e8,1<=n<=1e12

Hint

000 001
011 100
110 111
6种方案

思路:30%肯定是裸搜。对于特定的一种情况,我们的总方案数为 m n m^n mn,那么合法方案数显然不好想,但是我们可以反过来想一下不合法的方案数,很明显数量为m* ( m − 1 ) ( n − 1 ) (m-1)^{(n-1)} (m1)(n1),即第一个人任选一种,后面每个人都相对于前一个人而言有m-1种选择。那么ans就等于总方案数减去不合法方案数。注意要用快速幂取mod。

代码:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#define mod 100003
#define ll long long
using namespace std;
ll m,n,ans;
ll ksm(ll a,ll b)
{
	ll mul=1;
	while(b)
	{
	  if(b&1) mul=(mul*a)%mod;
	  a*=a,a%=mod,b>>=1;
	}
	return mul;
}
int main()
{
    scanf("%lld%lld",&m,&n);
    ans=(ksm(m,n)-(ksm(m-1,n-1)*m)%mod+mod)%mod;
    printf("%lld",ans);
	return 0;
} 
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值