2021暑假集训 Day 1

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_{ijk}$ 表示在此区间内能不能通过某种合并方式来获得这个数，那么怎么转移，我们可以先枚举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\ast b_j。$
定义矩阵的权值为矩阵内所有权值的和
给出限制[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*$(m-1{)}^{(n-1)}$，即第一个人任选一种，后面每个人都相对于前一个人而言有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;
}