数的三次方根
属于浮点数的二分法
*保留6位小数的话 r-l>1e-8 多2;
#include<iostream>
using namespace std;
int main(){
double x;
cin>>x;
double l=-10000,r=10000;
while(r-l>1e-8){
double mid=(l+r)/2;
if(mid*mid*mid>=x) r=mid;
else l=mid;
}
printf("%lf",l);
return 0;
}
前缀和
int输入的最大是108,如果超出输入的范围要采用longlong进行输入。
原数组:a1,a2,a3,……an;原数组下标从1开始
前缀和:si=a1+……ai;
问题1:如何求si
递归相加。
问题2:前缀和的作用
快速的求出某一段之和。[l,r]的和,可以表示成:sr-sl-1。
问题3:原数组从1开始
可以定义出s0,使得[1,x]可以用相同的方法进行计算。
例题:795、前缀和:输入长度为n的整数序列,输入m个询问,对于每个询问包含l,r输出l r序列和。
#include<iostream>
using namespace std;
const int N=100010;
int a[N],s[N];
int n,m;
int main(){
cin>>n>>m;
for(int i=1;i<=n;i++) cin>>a[i];
for(int i=1;i<=n;i++) s[i]=s[i-1]+a[i];
while(m--){
int l,r;
cin>>l>>r;
cout<<s[r]-s[l-1]<<endl;
}
return 0;
}
存在i-1的时候,下标一般从1开始
二维前缀和
应用:子矩阵的和
s[i,j]是i,j为下标的左上角的和。
//求s(i,j)
s[i,j]=s[i-1,j]+s[i,j-1]-s[i-1],j-1]+a[i,j];
//求(x1,y1)(x2,y2)构成矩阵的和
s[x2,y2]-s[x2,y1-1]-s[x1-1,y2]+s[x1-1,y1-1];
#include<iostream>
const int N=1010;
int n,m,q;
int a[N][N],s[N][N];
int main(){
scanf("%d%d%d",&n,&m,&q);
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
scanf("%d",&a[i][j]);
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
s[i][j]=s[i-1][j]+s[i][j-1]-s[i-1][j-1]+a[i][j];
while(q--){
int x1,y1,x2,y2;
scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
printf("%d\n",s[x2][y2]-s[x1-1][y2]-s[x2][y1-1]+s[x1-1][y1-1]);
return 0;
}
}
差分
相当于是前缀和的逆运算。
给定a[1],a[2],a[3],……a[n]
构造差分数组b[i]使得:
a[i]=b[1]+……b[i]
即a数组是b数组的前缀和。
b1=a1
b2=a2-a1
b3=a3–a2
差分的作用:
O(n)的复杂度内,由B数组可以得到A数组。将A数组的某一段加上某一个数的操作。
核心操作:将a[l~r]全部加上c等价于b[l]+=c;b[r+1]-=c;
#include<iostream>
using namespace std;
const int N=100010;
int a[N],b[N];
int n,m;
void insert(int l,int r,int c){
b[l]+=c;
b[r+1]-=c;
}
int main(){
scanf("%d%d",&a,&b);
for(int i=1;i<=n;i++) scanf("%d",&a[i]);
for(int i=1;i<=n;i++) insert(i,i,a[i]);
while(m--)
{
int l,r,c;
scanf("%d%d%d",&l,&r,&c);
insert(l,r,c);
}
for(int i=1;i<=n;i++)b[i]+=b[i-1];
for(int i=1;i<=n;i++) printf("%d\n",b[i]);
return 0;
}
差分矩阵
差分核心操作:给以(x1,y1)为左上角,(x2,y2)为右下角的子矩阵中所有的数a(i,j)加上c;
对于差分数组的影响:
s[x1,y1]+=c
s[x1,y2+1]-=c
s[x2+1,y1]-=c
s[[x2+1.y2+1]+=c
#include<iostream>
using namespace std;
const int N=1010;
int a[N][N],b[N][N];
int n,m,q;
void insert(int x1,int y1,int x2,int y2,int c){
b[x1][y1]+=c;
b[x1][y2+1]-=c;
b[x2+1][y1]-=c;
b[x2+1][y2+1]+=c;
return;
}
int main(){
cin>>n>>m>>q;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
scanf("%d",a[i][j]);
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++){
insert(i,j,i,j,a[i][j]);
}
while(q--){
int x1,y1,x2,y2,c;
scanf("%d%d%d%d%d",&x1,&y1,&x2,&y2,&c);
insert(x1,y1,x2,y2,c);
}
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
{a[i][j]=a[i-1][j]+a[i][j-1]-a[i-1][j-1]+b[i][j];
printf("%d ",a[i][j]);
}
return 0;
}