库存规划

0 阅读 作者:fogsail 2017-10-19

Rinky Dink公司是一家制造溜冰场冰面修整设备的公司。这种设备每个月的需求量都在变化,因此公司希望设计一种策略来规划生产,需求是给定的,即它虽然是波动的,但是是可预测的。公司希望设计接下来的$n$个月的生产计划。

对第$i$个月,公司知道需求$d_i$,即该月能够销售出去的设备的总量。令$D=sum_{i=1}^{n}d_i$
为后$n$个月的总需求。公司雇佣的全职员工,可以提供一个月制造$m$台设备的劳动力。如果公司希望一个月内制造多于$m$台的设备,可以雇佣额外的兼职劳动力,雇佣的成本为每制造一台机器付出$c$美元。而且,如果在月末还有设备未售出,公司将付出库存成本。
<!--more-->

保存$j$台设备的成本可以描述为一个函数$h(j)$,$j=1,2,cdots,D$,其中对所有的$1 leq j leq D$,$h(j) geq 0$,对$1 leq j leq D-1$,$h(j) leq h(j+1)$。

设计库存规划算法。

能量守恒的观点分析

库存规划的问题,可以用物理学中的能量守恒的观点进行分析

$store cost+HR cost=total value$
其中,$total value$表示总能量,总能量在$store cost$和$HR cost$中能够相互转换。

算法分析:

inventory_dynamic.h

#include <iostream>
#include <cmath>
using namespace std;

#define m 15
#define n 10
#define c 3


void find_create(int *extra_create,int *create_result,int *d)
{
    for(int i=0;i<=n;i++)
        create_result[i]=extra_create[i]+d[i];
}

int min_val(int a ,int b)
{
    return a<b?a:b;
}


int man_cost(int man_num)
{
    if(man_num<=m)
        return 0;
    else
        return c*(man_num-m);
}

int store_cost(int store_num)
{
    if(store_num>0)
        return (int)log((double)store_num);
    else
        return 0;
}

int create_material(int d[],int create[][n+1],int start,int end)
//start从0开始
{
    if(start==end)
    {
        create[start][end]=0;
    }
    if(end==start+1)
    {
        create[start][end]=d[end];
    }
    if(end>start)
    {
        create[start][end]=create_material(d,create,start,end-1)+d[end];
        //实现结果 d[start+1]+d[start+2]+.....+d[end]
    }
    else
    {
        return 0;
    }
    return create[start][end];
}

void init_create(int create[][n+1],int d[])
{
    for(int i=0;i<n;i++)
        create_material(d,create,i,n);
}

int min_cost(int d[],int HR_extra[],int start,int end)
{
    if(start>=end)
        return 0;

    int cost=0; //最初cost[]值为0,一开始每一个位置所保存的create[]就是当前需求值

    for(int i=start;i<=end;i++)
    {

        if(d[i]<=m)
        {
            //此时不需要花费额外的人力成本,所有的create[i]==d[i]均会销售出去
            cost+=0; //不需要花费任何代价,月末就可以出售完毕,没有库存
            HR_extra[i]=m-d[i];
        }
        else
        {
            HR_extra[i]=m-d[i];
            int HR_cost=d[i];  //人力资源成本超出预算,这部分预算可以考虑转换成库存成本
            int HR_cost_copy=d[i]; //副本,用来输出

            int cost_original=man_cost(HR_cost);
            int extra_store=0; //额外付出的库存代价
            //注意判断cost_original和cost_adjust哪个大?哪个小?

            for(int l=start;l<i;l++)  //(l,i)
            {
                int cur_extra=HR_extra[l];
                if(cur_extra>0 && HR_cost!=0)  //这部分人力资源成本可以用来存放库存,转换成库存成本
                {
                    int temp_excess=min_val(cur_extra,HR_cost);
                    //势能守恒来求解
                    HR_cost-=temp_excess;
                    extra_store+=temp_excess;
                }
            }

            int cost_adjust=man_cost(HR_cost)+store_cost(extra_store);

            if(cost_adjust<cost_original)
            {
                cost+=cost_adjust;
                //HR_extra也要调整
                for(int l=start;l<i;l++)
                {
                    if(HR_extra[l]>0 && HR_cost_copy!=0)
                    {
                        int temp_excess_copy=min_val(HR_extra[l],HR_cost_copy);

                        HR_cost_copy-=temp_excess_copy;

                        //保持势能守恒
                        HR_extra[l]-=temp_excess_copy;
                        HR_extra[i]+=temp_excess_copy;
                    }
                }
            }
            else
            {
                cost+=cost_original;
            }
        }
    }
    return cost;
}

inventory_dynamic.cpp

#include "inventory_dynamic.h"


int main()
{
    int result_cost,original_create;
    int d[n+1]={0,10,11,13,14,20,25,29,9,8,7};  //146

    int HR_extra[n+1]={0};
    int result[n+1]={0};

    int create[n+1][n+1]={0};
    //int cost[n+1][n+1]={0};

    original_create=create_material(d,create,0,n);
    init_create(create,d);

    result_cost=min_cost(d,HR_extra,1,n);

    cout<<"The result is :"<<endl;
    cout<<result_cost<<endl;

    cout<<"The create result is :"<<endl;

    for(int i=1;i<=n;i++)
    {
        cout<<m-HR_extra[i]<<" ";
    }
}

算法运行结果

原文地址:https://segmentfault.com/a/1190000009261639
广告一下
热门教程
PHP7报A non well formed numeric value encountered 0
Linux系统下关闭mongodb的几种命令分享 0
mongodb删除数据、删除集合、删除数据库的命令 0
Git&Github极速入门与攻坚实战课程 0
python爬虫教程使用Django和scrapy实现 0
libnetsnmpmibs.so.31: cannot open shared object file 0
数据结构和算法视频教程 0
redis的hash结构怎么删除数据呢? 0
C++和LUA解析器的数据交互实战视频 0
mongodb errmsg" : "too many users are authenticated 0
C++基础入门视频教程 0
用30个小时精通C++视频教程可能吗? 0
C++分布式多线程游戏服务器开发视频教程socket tcp boost库 0
C++培训教程就业班教程 0
layui的util工具格式时间戳为字符串 0
C++实战教程之远程桌面远程控制实战 1
网络安全培训视频教程 0
LINUX_C++软件工程师视频教程高级项目实战 0
C++高级数据结构与算法视频教程 0
跨域问题很头疼?通过配置nginx轻松解决ajax跨域问题 0
相关文章
【译】JavaScript数据结构(3):单向链表与双向链表 16
10个JavaScript难点 16
【译】苹果拒绝支持PWA,有损Web的未来 16
iView 一周年了,同时发布了 2.0 正式版,但这只是开始... 16
nodejs+mongodb构建一个简单登录注册功能 16
【译】JavaScript数据结构(4):树 16
组件化开发与黑箱 16
TypeScript - 不止稳,而且快 16
webpack3+anujs+ReactCSSTransitionGroup 16
原生js实现图片放大镜效果 16
WEB缓存探究第二弹——实战 16
纯笔记:vfork 的一些使用场景(顺便讲一下 fork 的原理) 16
Android APP 内部捐赠实现(支付宝&amp;微信) 16
WKWebView 的一些小总结 16
模型评价(一) AUC大法 16
开始使用GraphQL 16
Webpack模块化原理简析 16
gulp使用问题记录 16
使用Angular4动画为页面添彩 16
Python27 Matplotlib (win64 python2.7) 安装及简单使用 16