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=ai∗bj。
定义矩阵的权值为矩阵内所有权值的和
给出限制[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)} (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;
}