10.29 T2 magic 逆序对+前缀和

在这里插入图片描述
在这里插入图片描述

思路

  1. 先特判掉-1,-2
  2. 如何?? 对a和b分别求一下前缀和,如果sum[i]==sum[i+k]则说明这个区间(i-k)和为0
  3. 考虑修改其实就是把a的前缀和相邻两个元素交换顺序
  4. 那么问题转化为把1-n的排列a通过最少的交换次数得到给定的排列b,求最少方案数
  5. 这个问题用树状数组求(或归并排序)
  6. https://mp.csdn.net/mdeditor/83586170#
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
ll read()
{
	char ch=' ';
	ll f=1;ll x=0;
	while(ch<'0'||ch>'9')
	{
		if(ch=='-') f=-1;ch=getchar();
	}
	while(ch>='0'&&ch<='9')
	{
	    x=x*10+ch-'0';ch=getchar();
	}
	return x*f;
}
const ll N=1e5+100;
ll a[N],b[N],s1[N],s2[N],s3[N],s4[N];
ll x[N],y[N];
ll c[N];
ll lb(ll x)
{
	return x&(-x);
}
void modify(ll x)
{
	for(ll i=x;i<=N;i+=lb(i))
	{
		c[i]++;
	}
}
ll query(ll x)
{
	ll ret=0;
	for(ll i=x;i>=1;i-=lb(i))
	{
		ret+=c[i];
	}
	return ret;
}
int main()
{
	freopen("magic.in","r",stdin);
	freopen("magic.out","w",stdout);
	ll n;
	n=read();
	ll i,j;
	s1[0]=0;
	for(i=1;i<=n;i++)
	{
		a[i]=read();
		s1[i]=s1[i-1]+a[i];
	}
	for(i=1;i<=n;i++)
	{
		b[i]=read();
		s2[i]=s2[i-1]+b[i];
	}
	sort(s1+1,s1+1+n);
	sort(s2+1,s2+1+n);
	bool flag1=0;
	bool flag2=0;
	for(i=1;i<=n;i++)
	{
		if(s1[i-1]==s1[i]||s2[i-1]==s2[i])
		{
			flag1=true;
		}
		if(s1[i]!=s2[i])
		{
			flag2=true;
		}
		if(flag1&&flag2)
		{
			break;
		}
	}
	if(flag1)
	{
		cout<<"-1"<<endl;
		return 0;	
	}
	if(flag2)
	{
		cout<<"-2"<<endl;
		return 0;
	}
	for(i=1;i<=n;i++)
	{
		s3[i]=s3[i-1]+a[i];
	}
	for(i=1;i<=n;i++)
	{
		s4[i]=s4[i-1]+b[i];
	} 
	for(i=1;i<=n;i++)
	{
		ll pos=lower_bound(s1+1,s1+1+n,s3[i])-s1;
		x[pos]=i;
	}
	for(i=1;i<=n;i++)
	{
		ll pos=lower_bound(s2+1,s2+1+n,s4[i])-s2;
		y[i]=x[pos];
	}
	ll cnt=0;
	for(i=1;i<=n;i++)
	{
		ll tmp=n-y[i]+1;
		modify(tmp);
		cnt+=query(tmp-1);
	}
	cout<<cnt<<endl;
	return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
根据题意,可以画出一环形音乐厅的示意图,如下所示: ``` O | | | r2 = 13.72 m | | r1 = 5.145 m | ------------O------------ | | | r2 = 13.72 m | | O ``` 其中,演奏台位于两同心圆的中心,表示为 O 点。一侧墙的直径距离为 10.29 米,对应内同心圆的半径为 r1 = 5.145 米;另一侧墙的直径距离为 27.44 米,对应外同心圆的半径为 r2 = 13.72 米。 根据题目所给的公式 y(n)=x(n)+0.76x(n-N1)+0.5x(n-N2) 可以得到系统的单位冲激响应为: h(n) = δ(n) + 0.76δ(n-N1) + 0.5δ(n-N2) 其中,δ(n) 表示单位冲激信号,N1 和 N2 分别是内外两侧墙的回声延迟,可以通过计算距离和声音传输速度来得到: N1 = round(r1 / 343 * fs) N2 = round(r2 / 343 * fs) 其中,fs 是采样率。需要注意的是,由于 N1 和 N2 是整数,所以需要使用 round 函数进行四舍五入。 接下来,可以使用 Matlab 来进行仿真。首先,定义参数: ```matlab % 系统参数 r1 = 5.145; % 内同心圆半径,单位:米 r2 = 13.72; % 外同心圆半径,单位:米 v = 343; % 声音传输速度,单位:米/秒 fs = 44100; % 采样率,单位:赫兹 % 计算回声延迟 N1 = round(r1 / v * fs); N2 = round(r2 / v * fs); ``` 然后,可以生成单位冲激信号,并计算系统的单位冲激响应: ```matlab % 生成单位冲激信号 x = [1, zeros(1, fs - 1)]; % 计算系统的单位冲激响应 h = [1, zeros(1, N1 - 1), 0.76, zeros(1, N2 - N1 - 1), 0.5, zeros(1, fs - N2)]; ``` 最后,可以通过卷积来模拟音乐在一环形音乐厅中的效果: ```matlab % 读取音乐文件 [x, fs] = audioread('music.wav'); % 计算系统的输出 y = conv(x, h); % 播放系统的输出 sound(y, fs); ``` 在上面的代码中,使用了 Matlab 自带的读取音乐文件的函数 audioread,可以将音乐文件读入到 x 变量中。然后,使用 conv 函数对 x 和 h 进行卷积,得到系统的输出信号 y。最后,使用 Matlab 自带的播放音频的函数 sound,可以听到模拟音乐在一环形音乐厅中的效果。 需要注意的是,由于一环形音乐厅的形状比较特殊,这里只是简单地模拟了系统的单位冲激响应,实际上还需要考虑更多因素,比如音乐的频率响应、演奏者和听众的位置等等。

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值