FZU 2168 防守阵地 I(公式推导)(经典)(中等)

Problem 2168 防守阵地 I

Accept: 377    Submit: 1280
Time Limit: 3000 mSec    Memory Limit : 32768 KB

 Problem Description

部队中共有N个士兵,每个士兵有各自的能力指数Xi,在一次演练中,指挥部确定了M个需要防守的地点,按重要程度从低到高排序,依次以数字1到M标注每个地点的重要程度,指挥部将选择M个士兵依次进入指定地点进行防守任务,能力指数为X的士兵防守重要程度为Y的地点将得到X*Y的参考指数。现在士兵们排成一排,请你选择出连续的M个士兵依次参加防守,使得总的参考指数值最大。

 Input

输入包含多组数据。

输入第一行有两个整数N,M(1<=N<=1000000,1<=M<=1000),第二行N个整数表示每个士兵对应的能力指数Xi(1<=Xi<=1000)。

对于30%的数据1<=M<=N<=1000。

 Output

输出一个整数,为最大的参考指数总和。

 Sample Input

5 32 1 3 1 4

 Sample Output

17

思路:

这题肯定是不能暴力求解的,听说可以用什么线段树,我没试过。下面说说代码最少最快的公式推导

例:
5 3 
2 1 3 1 4
5个数,假设用a,b,c,d,e表示这5个数。有两个数组,sum[],s[](sum[i]表示前i 个数的和,s[i]表示前i 个sum[]的和)。
现在求1*a+2*b+3*c+4*d+5*e的和:
=a+b+c+d+e     sum[5]
+b+c+d+e         sum[5]-sum[1]
+c+d+e             sum[5]-sum[2]
+d+e                 sum[5]-sum[3]
+e                     sum[5]-sum[4]
即:
1*a+2*b+3*c+4*d+5*e=5*sum[5]-(sum[1]+sum[2]+sum[3]+sum[4])。
现在题目求连续三个数*1*2*3的最大和。接着看:
当i=3:
3*sum[3]=3*(a+b+c)
1*a+2*b+3*c=3*sum[3]-(2*a+b)=3*sum[3]-(a+b+a)=3*sum[3]-(sum[2]+sum[1])=3*sum[3]-(s[2]);
i=4:
3*sum[4]=3(a+b+c+d)
b+2*c+3*d=3*sum[4]-(3*a+2*b+c)=3*sum[4]-(sum[1]+sum[2]+sum[3])=3*sum[4]-(s[3])
这时候如果认为公式就是:

m*sum[i]-s[i-1]
那就错了;

你看,当i=5时:
c+2d+3e=3*sum[5]-(3*a+3*b+2*c+d)=3*sum[5]-(sum[4]+sum[3]+sum[2]+sum[1]-sum[1]);
这明显和上面推导的公式不符。
这里我就不接着推了,公式是:
m*sum[i]-(s[i-1]-s[i-1-m])

代码:
#include <cstdio>
#include <cstring>
#include<iostream>
#include <algorithm>
using namespace std;
const int L = 1000010;
int sum[L],s[L];

int main()
{
    int n,m,i,j,k,num;
    while(~scanf("%d%d",&n,&m))
    {
        sum[0]=s[0]=0;
        for(i=1;i<=n;i++)
        {
            scanf("%d",&num);
            sum[i]=sum[i-1]+num;
            s[i]=s[i-1]+sum[i];
        }
        if(n<=m)
        {
            printf("%d\n",sum[n]);
            continue;
        }
        int cnt,maxn=0;
        for(int i=m;i<=n;i++)
        {
            cnt=m*sum[i]-(s[i-1]-s[i-1-m]);
            maxn=max(maxn,cnt);
        }
        printf("%d\n",maxn);
    }
    return 0;
}






### 回答1: 你可以使用matlab中的hilbert函数来计算每一列的包络。具体步骤如下: 1. 加载数据矩阵 假设你的矩阵名字为data,包含n个样本,m个特征(列)。你可以使用load函数从文件中加载数据,或者手动创建一个矩阵。 2. 计算每一列的包络 使用hilbert函数计算每一列的复包络信号,再取模得到幅值,即可得到每一列的包络。代码如下: ```matlab env = abs(hilbert(data)); ``` 3. 绘制包络谱 使用matlab中的pwelch函数可以对每一列的包络信号进行功率谱估计,绘制包络谱。代码如下: ```matlab fs = 1000; %采样率 nfft = 2^nextpow2(size(env,1)); %FFT点数 window = hann(size(env,1)); %窗函数 noverlap = size(env,1)/2; %重叠点数 for i = 1:size(env,2) [Pxx, f] = pwelch(env(:,i), window, noverlap, nfft, fs); plot(f,Pxx); xlabel('Frequency (Hz)'); ylabel('PSD'); title(['Envelope Spectrum of Column ', num2str(i)]); pause(1); end ``` 上述代码中,先定义了采样率、FFT点数、窗函数和重叠点数等参数。然后对每一列的包络信号使用pwelch函数进行功率谱估计,并绘制包络谱。使用pause函数可以使得每一列的包络谱逐个显示出来。 完整代码如下: ```matlab %加载数据矩阵 load('data.mat'); %计算每一列的包络 env = abs(hilbert(data)); %绘制包络谱 fs = 1000; %采样率 nfft = 2^nextpow2(size(env,1)); %FFT点数 window = hann(size(env,1)); %窗函数 noverlap = size(env,1)/2; %重叠点数 for i = 1:size(env,2) [Pxx, f] = pwelch(env(:,i), window, noverlap, nfft, fs); plot(f,Pxx); xlabel('Frequency (Hz)'); ylabel('PSD'); title(['Envelope Spectrum of Column ', num2str(i)]); pause(1); end ``` ### 回答2: 要利用Matlab对一个矩阵中的每一列画出包络谱,可以按照以下步骤进行操作: 1. 首先,假设矩阵名为A,包含m行和n列。你可以使用Matlab中的size函数获取矩阵的行数m和列数n。 2. 使用一个循环从第一列开始遍历到最后一列,对每一列进行操作。 3. 在循环内,使用Matlab中的hilbert函数对当前列进行Hilbert变换。Hilbert变换用于将信号从时域转换为频域,并计算出信号的包络谱。 4. 对Hilbert变换后的信号应用Matlab中的abs函数,取其绝对值。 5. 对取得的包络谱进行可视化,使用Matlab中的plot函数绘制频率轴和包络谱的幅度。 6. 在绘制包络谱之前,可以使用Matlab中的fftshift函数将频谱进行中心化,以使低频部分位于中心位置。 7. 重复步骤2到步骤6直到所有列都被处理完毕。 8. 最后,可以使用Matlab中的标题、标签、图例等函数添加适当的文字描述和说明,以使绘图更加清晰和易读。 通过以上步骤,你可以通过Matlab对一个矩阵中的每一列画出包络谱。 ### 回答3: 利用Matlab对一个矩阵中的每一列画包络谱,我们可以按照以下步骤进行操作: 1. 首先,我们需要将矩阵中的每一列提取出来作为一个向量。 2. 接下来,我们可以利用matlab中的hilbert函数,将每一列向量转换为解析信号。 3. 然后,我们可以利用matlab中的envelope函数,求解每一列解析信号的包络谱。 4. 最后,我们可以利用matlab中的plot函数,将每一列包络谱绘制出来。 具体的Matlab代码如下: ``` % 假设我们的矩阵存储在变量A中,每一列代表一个信号 % 提取每一列作为向量 for i=1:size(A,2) x = A(:,i); % 转换为解析信号 xi = hilbert(x); % 求解包络谱 [up,lo] = envelope(xi); % 绘制包络谱 plot(abs(up)); hold on; end % 设置图例和标签 legend('Column 1', 'Column 2', ...); % 根据列的个数设置图例 xlabel('Frequency'); ylabel('Amplitude'); title('Envelope Spectrum of Matrix Columns'); ``` 以上代码可以对矩阵A中的每一列绘制包络谱,通过设置图例和标签,我们可以清晰地看到每一列信号的包络谱。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值