poj2112二分+最大流

本文详细介绍了SAP最大流算法的实现原理及具体代码实现,通过实例展示了如何使用该算法解决特定问题,并针对错误进行了深入分析。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

跟上题一样,因为一开始那份模板用的少,总怀疑那份sap模板出错了,今天换了一个模板做另外一题很像的题。(先前的模板不能处理过大的数)

事实证明我过程中处理错了,好吧,先放着,校赛完在来找错了。

 ////////////////////////////////////////////

最后在错误在建图的过程中,严格建图过程,因为增加了牛与牛和机器与机器的距离,导致求最大流结果时出错。

代码:

 

#include<iostream> #include<cstring> #include<cstdio> #define MAXN 405 #define INF 1e8 #define min(a,b) (a<b?a:b) #define max(a,b) (a>b?a:b) using namespace std; struct edge {     int u,v,w,next; }E[200000]; int head[MAXN],ecnt; int gap[MAXN],cur[MAXN],pre[MAXN],dis[MAXN]; int map[300][300]; int l,r,mid; int N,M,scr,sink,vn,num; int k,c1,m; void Floyd() {     for(int k=1;k<=sink-1;k++)     for(int i=1;i<=sink-1;i++)     if(i!=k&&map[i][k]!=-1)     for(int j=1;j<=sink-1;j++)     {         if(i!=j&&j!=k&&map[k][j]!=-1)         {             if(map[i][j]==-1)                 map[i][j]=map[i][k]+map[k][j];             else             map[i][j]=map[i][j]>map[i][k]+map[k][j]?map[i][k]+map[k][j]:map[i][j];         }     } } void Insert(int u,int v,int w) {     E[ecnt].u=u;     E[ecnt].v=v;     E[ecnt].w=w;     E[ecnt].next=head[u];     head[u]=ecnt++;     E[ecnt].u=v;     E[ecnt].v=u;     E[ecnt].w=0;     E[ecnt].next=head[v];     head[v]=ecnt++; } int Sap(int s,int t,int n)//核心代码(模版) {

    int ans=0,aug=INF;//aug表示增广路的流量     int i,v,u=pre[s]=s;     for(i=0;i<=n;i++)     {         cur[i]=head[i];         dis[i]=gap[i]=0;     }     gap[s]=n;     bool flag;     while(dis[s]<n)     {         flag=false;         for(int &j=cur[u];j!=-1;j=E[j].next)//一定要定义成int &j,why         {             v=E[j].v;             if(E[j].w>0&&dis[u]==dis[v]+1)             {                 flag=true;//找到容许边                 aug=min(aug,E[j].w);                 pre[v]=u;                 u=v;                 if(u==t)                 {                     ans+=aug;                     while(u!=s)                     {                         u=pre[u];                         E[cur[u]].w-=aug;                         E[cur[u]^1].w+=aug;//注意                     }                     aug=INF;                 }                 break;//找到一条就退出             }         }         if(flag) continue;         int mindis=n;         for(i=head[u];i!=-1;i=E[i].next)         {             v=E[i].v;             if(E[i].w>0&&dis[v]<mindis)             {                 mindis=dis[v];                 cur[u]=i;             }         }         if((--gap[dis[u]])==0) break;         gap[dis[u]=mindis+1]++;         u=pre[u];     }     return ans; } bool judge(int mid) {     memset(head,-1,sizeof(head));ecnt=0;         for(int i=k+1;i<=k+c1;i++)//源点         {             Insert(scr,i,1);         }         for(int i=1;i<=k;i++)//汇点         {             Insert(i,sink,m);         }            for(int i=k+1;i<=k+c1;i++)            for(int j=1;j<=k;j++)            if(map[i][j]<=mid&&map[i][j]!=-1)            {                Insert(i,j,INF);            }            //cout<<l<<' '<<r<<endl;            int ans=Sap(scr,sink,vn);

           return ans==c1; } int main() {    while(scanf("%d%d%d",&k,&c1,&m)!=EOF)    {        memset(map,-1,sizeof(map));        for(int i=1;i<=k+c1;i++)        {            for(int j=1;j<=k+c1;j++)            {                scanf("%d",&map[i][j]);                if(map[i][j]==0&&i!=j) map[i][j]=-1;            }        }

       scr=0;sink=k+c1+1;vn=sink+1;        Floyd();        /*for(int i=1;i<=k+c1;i++)        {            cout<<endl;            for(int j=1;j<=k+c1;j++)            {                cout<<map[i][j]<<' ';            }        }*/        int tmp=-1;        for(int i=1;i<=sink-1;i++)   //找最短路的最大值        for(int j=1;j<=sink-1;j++)        {            if(map[i][j]!=-1&&map[i][j]>tmp)            tmp=map[i][j];        }        l=0;r=0;mid=0;        r=tmp;        //cout<<endl<<tmp<<endl;        while(l!=r)        {            mid=(l+r)>>1;            if(judge(mid))   //大了            {                r=mid;            }            else l=mid+1;            //cout<<l<<' '<<r<<endl;        }        printf("%d\n",r);    }    return 0; }

 

资源下载链接为: https://pan.quark.cn/s/67c535f75d4c 在机器人技术中,轨迹规划是实现机器人从一个位置平稳高效移动到另一个位置的核心环节。本资源提供了一套基于 MATLAB 的机器人轨迹规划程序,涵盖了关节空间和笛卡尔空间两种规划方式。MATLAB 是一种强大的数值计算与可视化工具,凭借其灵活易用的特点,常被用于机器人控制算法的开发与仿真。 关节空间轨迹规划主要关注机器人各关节角度的变化,生成从初始配置到目标配置的连续路径。其关键知识点包括: 关节变量:指机器人各关节的旋转角度或伸缩长度。 运动学逆解:通过数学方法从末端执行器的目标位置反推关节变量。 路径平滑:确保关节变量轨迹连续且无抖动,常用方法有 S 型曲线拟合、多项式插值等。 速度和加速度限制:考虑关节的实际物理限制,确保轨迹在允许的动态范围内。 碰撞避免:在规划过程中避免关节与其他物体发生碰撞。 笛卡尔空间轨迹规划直接处理机器人末端执行器在工作空间中的位置和姿态变化,涉及以下内容: 工作空间:机器人可到达的所有三维空间点的集合。 路径规划:在工作空间中找到一条从起点到终点的无碰撞路径。 障碍物表示:采用二维或三维网格、Voronoi 图、Octree 等数据结构表示工作空间中的障碍物。 轨迹生成:通过样条曲线、直线插值等方法生成平滑路径。 实时更新:在规划过程中实时检测并避开新出现的障碍物。 在 MATLAB 中实现上述规划方法,可以借助其内置函数和工具箱: 优化工具箱:用于解决运动学逆解和路径规划中的优化问题。 Simulink:可视化建模环境,适合构建和仿真复杂的控制系统。 ODE 求解器:如 ode45,用于求解机器人动力学方程和轨迹执行过程中的运动学问题。 在实际应用中,通常会结合关节空间和笛卡尔空间的规划方法。先在关节空间生成平滑轨迹,再通过运动学正解将关节轨迹转换为笛卡
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值