洛谷P1667/[10.22 模拟赛] 数列 (思维+模拟)

洛谷P1667 数列

题目描述

给定一个长度是n的数列A,我们称一个数列是完美的,当且仅当对于其任意连续子序列的和都是正的。现在你有一个操作可以改变数列,选择一个区间[X,Y]满足\(A_X +A_{X+1} +…+ A_Y<0,1<X<=Y<n,\)\(S=A_X +A_{X+1} +…+ A_Y\),对于\(A_{X-1}\)\(A_{Y+1}\)分别加上S,\(A_X\)\(A_Y\)分别减去S(如果X=Y就减两次)。问最少几次这样的操作使得最终数列是完美的。

输入输出格式

输入格式:

第一行一个数n,以下n个数。

【数据规模】

对于20%的数据,满足1≤N≤5;

对于100%的数据,满足\(1≤N≤10^5; 1≤|A[i]|≤2^31-1.\)

输出格式:

一个数表示最少的操作次数,如果无解输出-1。

输入输出样例

输入样例#1:

5
13
-3
-4
-5
62

输出样例#1:

2

说明

【样例解释】

首先选择区间[2,4],之后数列变成1,9,-4,7,50,然后选择[3,3],数列变成1,5,4,3,50

Solution

按照题目意思,我们令\(T=sum[r]-sum[l-1]\),其中sum为a的前缀和

那么会有a[l-1]+=T,a[r+1]+=T,a[l]-=T,a[r]-=T,实际上对于sum[l]和sum[r+1]是没有变化的,而sum[l-1]会增加T,sum[r]会减少T,实际上就是sum[l-1]和sum[r]交换了位置

由于题目要求任意\(a_i\)均为正数,所以前缀和必须严格上升,那么很容易看出\(sum_i<=0\)或者是\(i<j\)并且\(sum_i=sum_j\)无解

正常情况下,我们要求交换次数,把前缀和离散后,它是第几小就该到哪去,所以就是模拟交换并统计次数就可以了

Code

#include<bits/stdc++.h>
#define rg register
#define il inline
#define Min(a,b) ((a)<(b)?(a):(b))
#define Max(a,b) ((a)>(b)?(a):(b))
#define lol long long
#define in(i) (i=read())
using namespace std;

const lol N=2e5+10;

lol read() {
    lol ans=0,f=1; char i=getchar();
    while(i<'0' || i>'9') {if(i=='-') f=-1; i=getchar();}
    while(i>='0' && i<='9') ans=(ans<<1)+(ans<<3)+i-'0',i=getchar();
    return ans*=f;
}

lol n,sum[N],id[N],AQ[N];

bool cmp(lol a,lol b) {return sum[a]<sum[b];}

int main()
{
    //freopen("bsum.in","r",stdin);
    //freopen("bsum.out","w",stdout);
    in(n);
    for(lol i=1;i<=n;i++) {
        in(sum[i]),id[i]=i;
        sum[i]+=sum[i-1],AQ[i]=sum[i];
    }
    sort(AQ+1,AQ+1+n);
    for(lol i=1;i<n;i++) {
        if(AQ[1]<=0 || AQ[i]==AQ[i+1])
            cout<<-1<<endl,exit(0);
    }
    sort(id+1,id+1+n,cmp);
    for(lol i=1;i<=n;i++) sum[id[i]]=i;
    lol ans=n;
    for(lol i=1;i<=n;i++) {
        if(sum[i]==i) ans--;
        else {
            swap(id[i],id[sum[i]]);
            swap(sum[i],sum[id[sum[i]]]);
        }
    }cout<<ans<<endl;
}

转载于:https://www.cnblogs.com/real-l/p/9832337.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
题目描述:一张灰度图像目标的边缘有两条,分别为A和B,边缘A上的点对应的法向量为NA,边缘B上的点对应的法向量为NB。试写出一种方法,用于估计这两条边缘的夹角。 解题思路:由于边缘A和边缘B都可以看做是一组点的集合,因此我们可以先计算出这两条边缘的法向量,然后再计算它们之间的夹角。 具体地,对于边缘A上的每个点p,可以计算出它的法向量NA(p),然后将所有的法向量取平均值,得到边缘A的平均法向量NA_avg。同理,对于边缘B上的每个点q,可以计算出它的法向量NB(q),然后将所有的法向量取平均值,得到边缘B的平均法向量NB_avg。最后,边缘A和边缘B的夹角θ可以通过计算它们的平均法向量之间的夹角得到,即: cosθ = dot(NB_avg, NA_avg) / (||NA_avg|| * ||NB_avg||) 其,dot(NB_avg, NA_avg)表示两个向量的点积,||NA_avg||和||NB_avg||表示两个向量的模长。 代码实现如下: ```python import numpy as np # 计算向量的模长 def norm(x): return np.sqrt(np.sum(x**2)) # 计算向量的夹角 def angle_between(u, v): cos_theta = np.dot(u, v) / (norm(u) * norm(v)) return np.arccos(cos_theta) # 计算边缘的平均法向量 def average_normal(points, normals): sum_normal = np.zeros_like(normals[0]) for i in range(len(points)): sum_normal += normals[i] avg_normal = sum_normal / len(points) return avg_normal # 边缘A的点集和法向量集合 points_A = [[10, 10], [20, 20], [30, 30]] normals_A = [[-1, 1], [-1, 1], [-1, 1]] # 边缘B的点集和法向量集合 points_B = [[10, 30], [20, 20], [30, 10]] normals_B = [[-1, -1], [-1, -1], [-1, -1]] # 计算边缘A和边缘B的平均法向量 NA_avg = average_normal(points_A, normals_A) NB_avg = average_normal(points_B, normals_B) # 计算边缘A和边缘B的夹角 theta = angle_between(NA_avg, NB_avg) print("边缘A和边缘B的夹角为:", np.degrees(theta)) ``` 上述代码,points_A和points_B分别表示边缘A和边缘B上的点集,normals_A和normals_B分别表示边缘A和边缘B上的法向量集合,NA_avg和NB_avg分别表示边缘A和边缘B的平均法向量,theta表示边缘A和边缘B的夹角。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值