复习计划-01分数规划

最新推荐文章于 2023-05-24 20:55:34 发布
最新推荐文章于 2023-05-24 20:55:34 发布
阅读量97 收藏
点赞数
分类专栏： 二分
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Zenith_Habitant/article/details/108808695
版权

01分数规划

01分数规划的题目一般比较明显，是一个sigma(a[i] * x[i]) / sigma(b[i] * x[i]) 的形式，这类题目套路比较固定，写出数学表达式之后，二分就好了，有时候需要一些贪心，二分的Check需要根据表达式判断大于0还是小于0

  • 经典例题：P1730
  • 二分 + 最短路
//#define LOCAL
#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define mem(a, b) memset(a,b,sizeof(a))
#define sz(a) (int)a.size()
#define INF 0x3f3f3f3f
#define DNF 0x7f
#define DBG printf("this is a input\n")
#define fi first
#define se second
#define mk(a, b) make_pair(a,b)
#define pb push_back
#define LF putchar('\n')
#define SP putchar(' ')
#define p_queue priority_queue
#define CLOSE ios::sync_with_stdio(0); cin.tie(0)

template<typename T>
void read(T &x) {x = 0;char ch = getchar();ll f = 1;while(!isdigit(ch)){if(ch == '-')f *= -1;ch = getchar();}while(isdigit(ch)){x = x * 10 + ch - 48; ch = getchar();}x *= f;}
template<typename T, typename... Args>
void read(T &first, Args& ... args) {read(first);read(args...);}
template<typename T>
void write(T arg) {T x = arg;if(x < 0) {putchar('-'); x =- x;}if(x > 9) {write(x / 10);}putchar(x % 10 + '0');}
template<typename T, typename ... Ts>
void write(T arg, Ts ... args) {write(arg);if(sizeof...(args) != 0) {putchar(' ');write(args ...);}}
using namespace std;

ll gcd(ll a, ll b) {
    return b == 0 ? a : gcd(b, a % b);
}

ll lcm(ll a, ll b) {
    return a / gcd(a, b) * b;
}

const int N = 55, M = 1005;
const double eps = 1e-6;
int n , m;
int head[N] , cnt;
struct node {
    int t, next;
    double w;
}edge[M << 1];
void add (int f, int t, double w)
{
    edge[cnt].t = t;
    edge[cnt].w = w;
    edge[cnt].next = head[f];
    head[f] = cnt ++;
}
int a[N][N];
double ans[N][N];
void floyd ()
{
    for (int k = 1 ; k <= n ; k ++)
        for (int i = 1 ; i <= n ; i ++)
            for (int j = 1 ; j <= n ; j ++)
                a[i][j] |= (a[i][k] && a[k][j]);
}
double dis[N];
int vis[N], times[N];
bool spfa (int s, int t, double x)
{
    queue <int> q;
    for (int i = 1 ;i <= n ; i ++)
        dis[i] = 0x3f3f3f3f, vis[i] = 0;
    dis[s] = 0.0;
    vis[s] = 1;
    q.push(s);
    while (!q.empty())
    {
        int u = q.front();
        q.pop();
        vis[u] = 0;
        for (int i = head[u] ; i != -1 ; i = edge[i].next)
        {
            int v = edge[i].t;
            double w = edge[i].w;
            if (dis[v] > dis[u] + w - x)
            {
                dis[v] = dis[u] + w - x;
                if (!vis[v])
                    q.push(v) , vis[v] = 1;
            }
        }
    }
    return dis[t] > 0;
}
void solve ()
{
    for (int i = 1 ; i <= n ; i ++)
    {
        for (int j = 1 ; j <= n ; j ++)
        {
            if (a[i][j] && i != j)
            {
                double l = 0.0 , r = 100000.0;
                while(r-l > 1e-6){
                    double mid = (l+r)/2;
                    if(spfa(i,j,mid))
                        l=mid;
                    else
                        r=mid;
                }
                ans[i][j]  = l;
            }
        }
    }
}
int main (void)
{
    mem(head,-1);
    read (n, m);
    for (int i = 1 ; i <= m ; i ++)
    {
        int u , v;
        double w;
        read (u , v);
        scanf("%lf",&w);
        add (u , v , w);
        a[u][v] = 1;
    }
    floyd ();
    solve ();
    int Q;
    read (Q);
    while (Q --)
    {
        int u , v;
        read (u , v);
        if (a[u][v])
            printf("%.3lf\n", ans[u][v]);
        else
            printf("OMG!\n");
    }
}
  • AcWing361 01分数划分+spfa判负环
//#define LOCAL
#include <bits/stdc++.h>
 
using namespace std;
#define ll long long
#define mem(a, b) memset(a,b,sizeof(a))
#define INF 0x3f3f3f3f
#define DNF 0x7f
#define DBG printf("this is a input\n")
#define fi first
#define se second
#define mk(a, b) make_pair(a,b)
#define p_queue priority_queue
 
ll gcd(ll a, ll b) {
    return b == 0 ? a : gcd(b, a % b);
}
 
ll lcm(ll a, ll b) {
    return a / gcd(a, b) * b;
}
int n, m;
int head[1005], cnt;
struct e{
    int t, next;
    double w;
}edge[10005];
void add(int f, int t, double w)
{
    edge[cnt].t = t;
    edge[cnt].w = w;
    edge[cnt].next = head[f];
    head[f] = cnt ++;
}
double ve[1005];
double f[1005][1005];
bool vis[1005];
double dis[1005];
int times[1005];
bool check(double x)
{
    mem(dis,DNF);
    mem(times,0);
    mem(vis,0);
    queue <int> q;
    dis[1] = 0.0;
    vis[1] = 1;
    q.push(1);
    while(!q.empty())
    {
        int no = q.front();
        q.pop();
        vis[no] = 0;
        for(int i = head[no] ; i != -1 ; i = edge[i].next)
        {
            int v = edge[i].t;
            double w = x * edge[i].w - ve[no];
            if(dis[v] > dis[no] + w)
            {
                dis[v] = dis[no] + w;
                times[v] = times[no] + 1;
                if(times[v] >= n)
                    return true;
                if(!vis[v])
                {
                    q.push(v);
                    vis[v] = 1;
                }
            }
        }
    }
    return false;
}
int main(void)
{
#ifdef LOCAL
    freopen("data.in", "r", stdin);
    freopen("odata.out", "w", stdout);
#endif
    mem(head,-1);
    int u ,v;
    double w;
    scanf("%d %d",&n , &m);
    for(int i = 1 ; i <= n ; i ++)
        scanf("%lf",&ve[i]);
    for(int i = 1 ; i <= m ; i ++)
    {
        scanf("%d %d %lf",&u,&v,&w);
        add(u,v,w);
    }
    double l = 0.0, r = 1000.0;
    for(int i = 1 ; i <= 50 ; i ++)
    {
        double mid = (l + r) / 2.0;
        if(check(mid))
            l = mid;
        else
            r = mid;
    }
    printf("%.2lf\n",r);
}
  • AcWing348 最优比率生成树
#include <bits/stdc++.h>
 
using namespace std;
#define ll long long
#define mem(a, b) memset(a,b,sizeof(a))
#define INF 0x3f3f3f3f
#define DBG printf("this is a input\n")
#define fi first
#define se second
#define mk(a, b) make_pair(a,b)
#define p_queue priority_queue
 
ll gcd(ll a, ll b) {
    return b == 0 ? a : gcd(b, a % b);
}
 
ll lcm(ll a, ll b) {
    return a / gcd(a, b) * b;
}
int n, tot;
struct po{
    double x , y, h;
}point[1010];
struct e
{
    int u, v;
    double cost, val;
    double w;
    bool operator < (const e& no) const{
        return  w < no.w;
    }
}edge[1010];
double dp[1010];
int vis[1010];
double dis(double x1, double y1 , double x2, double y2) {
    return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
bool prime(double R)
{
    for(int i = 1 ; i <= n ; i ++)
        dp[i] = 10000001.0;
    mem(vis,0);
    double sum = 0.0;
    dp[1] = 0.0;
    for(int i = 1 ; i <= n ; i ++)
    {
        int index = 1;
        double minn = INF;
        for(int j = 1 ; j <= n ; j ++)
        {
            if(dp[j] < minn && !vis[j])
            {
                minn = dp[j];
                index = j;
            }
        }
        vis[index] = 1;
        sum += minn;
        for(int j = 1 ; j <= n ; j ++)
        {
            double cost = dis(point[index].x, point[index].y, point[j].x, point[j].y);
            double h = fabs(point[index].h - point[j].h);
            if(!vis[j] && dp[j] > h - cost*R + 1e-6) {
                dp[j] = h - cost * R;
            }
        }
    }
    if(sum < 0.0) //最小 反之最大
        return true;
    return false;
}
int main(void)
{
    while(~scanf("%d",&n),n) {
        tot = 1;
        for (int i = 1; i <= n; i++)
            scanf("%lf %lf %lf",&point[i].x,&point[i].y,&point[i].h);
        double l = 0.0, r = 10000001.0;
        for (int i = 1; i <= 40; i++) {
            double mid = (l + r) / 2.0;
            if ((prime (mid)))
                r = mid;
            else
                l = mid;
        }
        printf("%.3lf\n", r);
    }
}
  • P1768
  • 最优比率环：SPFA判断负环+二分（推式子）
//#define LOCAL
#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define mem(a, b) memset(a,b,sizeof(a))
#define sz(a) (int)a.size()
#define INF 0x3f3f3f3f
#define DNF 0x7f
#define DBG printf("this is a input\n")
#define fi first
#define se second
#define mk(a, b) make_pair(a,b)
#define pb push_back
#define LF putchar('\n')
#define SP putchar(' ')
#define p_queue priority_queue
#define CLOSE ios::sync_with_stdio(0); cin.tie(0)

template<typename T>
void read(T &x) {x = 0;char ch = getchar();ll f = 1;while(!isdigit(ch)){if(ch == '-')f *= -1;ch = getchar();}while(isdigit(ch)){x = x * 10 + ch - 48; ch = getchar();}x *= f;}
template<typename T, typename... Args>
void read(T &first, Args& ... args) {read(first);read(args...);}
template<typename T>
void write(T arg) {T x = arg;if(x < 0) {putchar('-'); x =- x;}if(x > 9) {write(x / 10);}putchar(x % 10 + '0');}
template<typename T, typename ... Ts>
void write(T arg, Ts ... args) {write(arg);if(sizeof...(args) != 0) {putchar(' ');write(args ...);}}
using namespace std;

ll gcd(ll a, ll b) {
    return b == 0 ? a : gcd(b, a % b);
}

ll lcm(ll a, ll b) {
    return a / gcd(a, b) * b;
}
const int N = 5e4 + 50;

int n , m;
int head[N] , cnt;
struct node {
    int t, next;
    double v , p;
}edge[N];
void add (int f, int t, double v, double p)
{
    edge[cnt].t = t;
    edge[cnt].v = v;
    edge[cnt].p = p;
    edge[cnt].next = head[f];
    head[f] = cnt ++;
}
double dis[N];
int vis[N], times[N], in[N];
bool spfa (double x)
{

    queue<int>q;
    for (int i = 0 ; i <= n ; i ++)
        dis[i] = 1000000000000.0 , vis[i] = times[i] = 0;
    dis[0] = 0.0;
    vis[0] = 1;
    q.push(0);
    while (!q.empty())
    {
        int u = q.front();
        q.pop();
        vis[u] = 0;
        for (int i = head[u] ; i != -1 ; i = edge[i].next)
        {
            int t = edge[i].t;
            double v = edge[i].v , p = edge[i].p;
            double w = 1.0 * x * p - 1.0 * v;
            if (dis[t] > dis[u] + w)
            {
                dis[t] = dis[u] + w;
                times[t] = times[u] + 1;
                if (times[t] >= n+1)
                    return true;
                if (!vis[t])
                    q.push(t), vis[t] = 1;
            }
        }
    }
    return false;
}
bool topo()
{
    int tot = 0;
    queue <int> q;
    for (int i = 1 ; i <= n ; i ++)
        if (!in[i])
            q.push(i), tot ++;
    while (!q.empty())
    {
        int u = q.front();
        q.pop ();
        for (int i = head[u] ; i != -1 ; i = edge[i].next)
        {
            int t = edge[i].t;
            in[t] --;
            if (!in[t])
                q.push(t), tot ++;
        }
    }

    if (tot == n)
        return true;
    return false;

}
int main (void)
{
    mem (head,-1);
    read (n , m);
    for (int i = 1 ; i <= m ; i ++)
    {
        int u , v;
        double a, b;
        scanf("%d %d %lf %lf",&u ,&v , &a , &b);
        add (u , v , a, b);
        in[v] ++;
    }
    for (int i = 1; i <= n ; i ++)
        add (0 , i , 0 , 1000000);
    if (topo())
    {
        printf("-1\n");
        return 0;
    }
    double l = 0.0 , r = 200.0;
    while (r - l > 1e-6)
    {
        double mid = (l+r)/2.0;
        if(spfa(mid))
            l=mid;
        else
            r=mid;
    }
    printf("%.1lf\n",l);
}
卡常带师
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
2021考研——复习规划（408篇）
Mr.鹏
04-29 2万+
计算机考研408是全国统考的科目，统考科目包括四门计算机专业课：数据结构、计算机组成原理、操作系统和计算机网络，这四门课程合在一起称为计算机科学专业基础综合，共150分。四门专业课在试卷中所占的分数分别为：数据结构45分，计算机组成原理45分，操作系统35分，计算机网络25分。
期末复习----习题
weixin_51676760的博客
01-23 2169
期末复习一 1.(10分) 从键盘输入任意一个3位整数，编程计算并输出它的逆序数（忽略整数前的正负号）。例如，输入-123，则忽略负号，由123分离出百位1、十位2、个位3，然后计算3100+210+1=321，输出321. 题目要求： 输入提示为Input x: 输出格式为y=%d\n #include <stdio.h> #include <math.h> int main(void) { int x = 0, y = 0; printf("Input x:")
机器学习R复习提纲-TYUT
weixin_51603782的博客
05-24 3142
复习备考
考研数据结构复习---线性表
qq_43499647的博客
05-24 2820
考纲内容 1.线性表的定义和基本操作 2.线性表的实现 顺序存储;链式存储;线性表的应用 本章知识框架 复习提示 线性表是考研命题的重点.这类算法题实现起来比较容易,代码量少,但是要求具有最优的性能(时间,空间复杂度),才能获得好的分数. 线性表目录1.线性表的定义和基本操作1.1线性表的定义1.2线性表的基本操作2.线性表的顺序表示2.1顺序表的定义2.2顺序表上基本操作的实现3.线性表的链式表示3.1单链表的定义3.2单链表上基本操作的实现3.3双链表3.4循环链表3.5静态链表3.6顺序表.
考研复习计划
qq_44396242的博客
06-01 1846
从6月初开始，以计算机专硕22408（202英语二、302数学二、408计算机学科专业基础综合）为例。
数据库作业复习 -- 子查询（自用）
SweetSky_Lee的博客
04-04 2300
1. 查询每门被选修的课程信息，输出课程号，课程名 select distinct sc.cno,cname from sc,course where sc.cno = course.cno order by sc.cno desc 2. 查询每门课程被选修的情况，输出课程号，课程名，被选修次数 select sc.cno,cname,count(*) from sc,course where sc.cno = course.cno group by ...
Linux期末复习-ch6-Shell
AlphaCat的博客
12-05 1360
! 与一个逻辑值相反的逻辑值 -a 两个逻辑值为“是”返回值才为“是”，反之为“否” -o 两个逻辑值有一个为“是”，返回值就为“是” 【例6.11】 使用if-then 语句创建简单的Shell程序。 使用vi编辑器创建Shell程序，文件名为bbbb，文件内容如下所示。 #!/bin/bash #filename:bbbb echo -n "Do you want to continue: Y or N" read ANSWER if [ $ANSWER=N -o $ANSWER=n ] then ex
CSP-S前的复习计划 +记录+总结
lahlah的博客
10-02 1623
一只菜鸡CSP-S前的挣扎
好分数阅卷3.0_好分数app下载-好分数网查成绩安卓版app v3.6.4.1-清风安卓软件网...
weixin_34188805的博客
12-23 4367
好分数网查成绩安卓版app：是一个提供快速成绩查询服务的平台，让用户可以随时随地查看成绩，以及分析离自己理想的大学相差的分数。在这里，可以查看自己的试卷，分析错题，快速提高学习。好分数app，查成绩也是学习的一大动力哦，真心希望每位同学考上理想的大学！欢迎喜欢的朋友下载体验！好分数app特色： 1、老师发布成绩后，学生第一时间快速查分；2、可以进行多维度的排名比较；3、用平时考试成绩来预测高考成绩...
中学《教育知识与能力》复习笔记-教资复习全靠它~背诵资料
weixin_45797533的博客
08-23 8606
（wps没会员凑乎看）
六年数学上册总复习题-分数乘除法四则计算精选.doc
09-26
六年数学上册总复习题-分数乘除法四则计算精选.doc
《分数乘法整理复习》教学反思.docx
01-07
《分数乘法整理复习》教学反思
六年级下册数学总复习试题-分数的基本性质专项练（通用版 含答案）.docx
06-06
六年级下册数学总复习试题-分数的基本性质专项练（通用版 含答案）.docx
六年级下册数学总复习试题-分数大小的比较专项练（通用版 含答案）.docx
06-06
六年级下册数学总复习试题-分数大小的比较专项练（通用版 含答案）.docx
六年级下册数学总复习试题-分数加减法的计算专项练通用版（含答案）.docx
06-06
六年级下册数学总复习试题-分数加减法的计算专项练通用版（含答案）.docx
“推荐系统”相关资源推荐
06-14
推荐了国内外对推荐系统的讲解相关资源
全渠道电商平台业务中台解决方案.pptx
06-14
全渠道电商平台业务中台解决方案.pptx
云计算企业私有云平台建设方案.pptx
06-14
云计算企业私有云平台建设方案.pptx
通过CNN卷积神经网络对盆栽识别-含图片数据集.zip
最新发布
06-14
本代码是基于python pytorch环境安装的。 下载本代码后，有个requirement.txt文本，里面介绍了如何安装环境，环境需要自行配置。 或可直接参考下面博文进行环境安装。 https://blog.csdn.net/no_work/article/details/139246467 如果实在不会安装的，可以直接下载免安装环境包，有偿的哦 https://download.csdn.net/download/qq_34904125/89365780 安装好环境之后， 代码需要依次运行 01数据集文本生成制作.py 02深度学习模型训练.py 和03pyqt_ui界面.py 数据集文件夹存放了本次识别的各个类别图片。 本代码对数据集进行了预处理，包括通过在较短边增加灰边，使得图片变为正方形(如果图片原本就是正方形则不会增加灰边)，和旋转角度，来扩增增强数据集， 运行01数据集文本制作.py文件，会就读取数据集下每个类别文件中的图片路径和对应的标签 运行02深度学习模型训练.py就会将txt文本中记录的训练集和验证集进行读取训练，训练好后会保存模型在本地
英语四级合理复习方式规划
04-20
针对英语四级的合理复习方式规划，以下几点可能对你有所帮助： 1. 定期模拟考试：模拟考试可以让你更好地了解自己的考试水平，以及哪些方面需要更加用心学习。建议每周至少模拟一次考试。 2. 多练习听力和阅读理解...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值