day 2 数学基础 &数据结构1

本文介绍了算法中的质数判断、质数筛、并查集等基础数据结构,并详细讲解了树状数组、逆序对及线段树的实现方法。

好吧,实际上好多我也不会

质数

bool is_prime(int n)
{
    if(n<2) return flase;
    for(int i=2;i<=sqrt(n);i++)
    if(n%i==0) return flase;
    return true;
    }
    

质数筛

基本上用欧式晒

void prime(int n)
{
	memset(v,0,sizeof(v));
	for(int i=2;i<=n;i++)
	{
		if(v[i]) continue;
		cout<<i<<endl;
		for(int j=2;j<=n/i;j++) v[i*j]=1;
	}
}

剩下的估计初赛也不考

数据结构

并查集

int find(int x)
{
	if(x==fa[x])return  fa[x];
	return fa[x]=get(fa[x]);
}
void merge(int x,int y)
{
	fa[get(x)]=get(y);
}

另外还有扩展与,边带全
在我博客,不再重复

树状数组

#include<bits/stdc++.h>
using namespace std;
inline long long read(){
    long long num=0;int z=1;char c=getchar();
    if(c=='-') z=-1;
    while((c<'0'||c>'9')&&c!='-') c=getchar();
    if(c=='-') z=-1,c=getchar();
    while(c>='0'&&c<='9') num=(num<<1)+(num<<3)+(c^48),c=getchar();
    return z*num;
}
int m,n;
int a;
int tree[2000005];
int lowbit(int x)
{
    return x&-x;
}
long long add(int x,int k)
{
    while(x<=n)
    {
        tree[x]+=k;
        x+=lowbit(x);
    }
}
long long ask(int x)
{
    long long ans=0;
    while(x!=0)
    {
        ans+=tree[x];
        x-=lowbit(x);
    }
    return ans;
}
int main(){
    n=read();m=read();
    int last=0,now;
    for(int i=1;i<=n;i++)
    {
        now=read();
        add(i,now-last);
        last=now;
    }
    for(int i=1;i<=m;i++)
    {
        int o,x,y,k;
        o=read();
        if(o==1)
        {
            x=read();y=read();k=read();
            add(x,k);
            add(y+1,-k);//区间
        }
        else{
            x=read();
            cout<<ask(x)<<endl;
        }
    }
    return 0;
}

逆序对

楼兰图图,很经典

for(int i=n;i;i--)
{
	ans+=ask(a[i]-1);
	add(a[i],1);
}
#include <cstdio>
#include <cstring>
#include <iostream>
#define ll long long
using namespace std;
const int N = 200006;
int n, a[N], c[N], l[N], r[N];

int ask(int x) {
    int ans = 0;
    while (x) {
        ans += c[x];
        x -= x & -x;
    }
    return ans;
}

void add(int x, int y) {
    while (x <= n) {
        c[x] += y;
        x += x & -x;
    }
}

int main() {
    cin >> n;
    for (int i = 1; i <= n; i++) scanf("%d", &a[i]);
    memset(c, 0, sizeof(c));
    for (int i = n; i; i--) {
        r[i] = ask(a[i] - 1);
        add(a[i], 1);
    }
    memset(c, 0, sizeof(c));
    for (int i = 1; i <= n; i++) {
        l[i] = ask(a[i] - 1);
        add(a[i], 1);
    }
    ll ans = 0;
    for (int i = 1; i <= n; i++)
        ans += (ll)(n - i - r[i]) * (i - 1 - l[i]);
    cout << ans << " ";
    ans = 0;
    for (int i = 1; i <= n; i++) ans += (ll)r[i] * l[i];
    cout << ans << endl;
    return 0;
}

作者:ruanmowen
链接:https://www.acwing.com/activity/content/code/content/323834/
来源:AcWing
著作权归作者所有。商业转载请联系作者获得授权,非商业转载请注明出处。

线段树

注意4倍空间

struct segmentTree
{
	int l,r;
	int dat;
}t[SIZE*4];
void build(int p,int l,int r)
{
	t[p].l=l,t[p].r=r;
	if(l==r) {
		t[p].dat=a[l];return ;
	}
	int mid=l+r>>1;
	build(p*2,l,mid);
	build(p*2+1,mid+1,r);
	t[p].dat=t[p*2+1].dat+t[p*2].dat;
}

build(1,1,n);
void change(int p,int x,int v)
{
	if(t[p].l==t[p].r)
	{
		t[p].dat=v;return ;
	}
	int mid=t[p].l+t[p].r>>1;
	if(x<=mid)change(p*2,x,v);
	else change(p*2+1,x,v);
	t[p].dat=t[p*2].dat+t[p*2+1].dat;
}

change(1,x,v); 
【四轴飞行器】非线性三自由度四轴飞行器模拟器研究(Matlab代码实现)内容概要:本文围绕非线性三自由度四轴飞行器模拟器的研究展开,重点介绍了基于Matlab的建模与仿真方法。通过对四轴飞行器的动力学特性进行分析,构建了非线性状态空间模型,并实现了姿态与位置的动态模拟。研究涵盖了飞行器运动方程的建立、控制系统设计及数值仿真验证等环节,突出非线性系统的精确建模与仿真优势,有助于深入理解飞行器在复杂工况下的行为特征。此外,文中还提到了多种配套技术如PID控制、状态估计与路径规划等,展示了Matlab在航空航天仿真中的综合应用能力。; 适合人群:具备一定自动控制理论基础和Matlab编程能力的高校学生、科研人员及从事无人机系统开发的工程技术人员,尤其适合研究生及以上层次的研究者。; 使用场景及目标:①用于四轴飞行器控制系统的设计与验证,支持算法快速原型开发;②作为教学工具帮助理解非线性动力学系统建模与仿真过程;③支撑科研项目中对飞行器姿态控制、轨迹跟踪等问题的深入研究; 阅读建议:建议读者结合文中提供的Matlab代码进行实践操作,重点关注动力学建模与控制模块的实现细节,同时可延伸学习文档中提及的PID控制、状态估计等相关技术内容,以全面提升系统仿真与分析能力。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值