SLMCLP.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "def rastojanje(A,B):\n",
    "    x1,y1 = A\n",
    "    x2,y2 = B\n",
    "    \n",
    "    return np.sqrt((x1-x2)**2 + (y1-y2)**2)\n",
    "\n",
    "def mat_rastojanja(niz):\n",
    "   \n",
    "    duzina = len(niz)\n",
    "    M = np.eye(duzina)\n",
    "    #print(duzina)\n",
    "    for i in range(duzina):\n",
    "        for j in range(duzina):\n",
    "            \n",
    "            M[i,j] = rastojanje(niz[i],niz[j])\n",
    "            \n",
    "    return np.array(M)\n",
    "\n",
    "\n",
    "def generisi(broj_tacaka):\n",
    "    \n",
    "    x = np.random.randint(0,100,broj_tacaka)\n",
    "    y = np.random.randint(0,100,broj_tacaka)\n",
    "    \n",
    "   \n",
    "   # x = np.random.random_sample(size = broj_tacaka)*100\n",
    "   #y = np.random.random_sample(size = broj_tacaka)*100\n",
    "    \n",
    "    return np.array([(a,b) for a,b in zip(x,y)])\n",
    "\n",
    "\n",
    "def get_locations(m,t):\n",
    "    \n",
    "    sve = list(range(t))\n",
    "    np.random.shuffle(sve)\n",
    "    \n",
    "    \n",
    "    facilities = np.sort(sve[:m]).tolist()\n",
    "    customers = np.sort(sve[m:]).tolist()\n",
    "    \n",
    "    return (facilities, customers)\n",
    "    \n",
    "\n",
    "def separate(facilities,deo):\n",
    "    \n",
    "    np.random.shuffle(facilities)\n",
    "    k = int(np.round(len(facilities)/deo))\n",
    "    \n",
    "    I2 = np.sort(facilities[:k]).tolist()\n",
    "    I1 = np.sort(facilities[k:]).tolist()\n",
    "   \n",
    "    return (I1,I2)\n",
    "\n",
    "\n",
    "\n",
    "def pref(customers,facilities, D):\n",
    "    x= np.eye(len(customers), len(facilities))\n",
    "   # print('matrica rastijanja', D)\n",
    "    D = np.array(D)\n",
    "    u = np.random.uniform(size = 1)[0]\n",
    "    for i in range(len(customers)):\n",
    "        \n",
    "        \n",
    "       # print(i+1,\"red\",D[customers[i]][facilities])\n",
    "        m = min(D[customers[i]][facilities])\n",
    "        M = max(D[customers[i]][facilities])\n",
    "        \n",
    "       # print(\"min {0:.2f}\".format(m))\n",
    "       # print(\"max{0:.2f}\".format(M))\n",
    "        \n",
    "        for j in range(len(facilities)):\n",
    "            k = (D[i,j] - m)/(M-m)    \n",
    "            \n",
    "            if 0<=u<=k:\n",
    "                x[i,j] = m + np.sqrt(u*(M-m)*(D[customers[i]][facilities[j]] - m))\n",
    "            else:\n",
    "                x[i,j] = M - np.sqrt((M-D[customers[i]][facilities[j]])*(M-m)*(1-u))\n",
    "                \n",
    "    return x\n",
    "\n",
    "def podeli_kupce(tacke, radius,customers, I2):\n",
    "    \n",
    "    J1 = []\n",
    "    J2 = []\n",
    "    \n",
    "    for i in customers:\n",
    "        for j in I2:\n",
    "            ind = False\n",
    "            \n",
    "            if rastojanje(tacke[i], tacke[j]) < radius:\n",
    "                ind = True\n",
    "                break\n",
    "            \n",
    "        if ind:\n",
    "            J2.append(i)\n",
    "        else:\n",
    "            J1.append(i)\n",
    "            \n",
    "    return (J1,J2)\n",
    "            \n",
    "               "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "def foo(radius,p,t,m,dat,deo = 10):\n",
    "\n",
    "\n",
    "\n",
    "    tacke = generisi(t)\n",
    "    #print(tacke)\n",
    "    D = mat_rastojanja(tacke)\n",
    "\n",
    "\n",
    "\n",
    "    facilities, customers = get_locations(m,t)\n",
    "    #print(facilities)\n",
    "    #print(customers)\n",
    "    I1,I2 = separate(facilities,deo)\n",
    "    I = I1 + I2 # Nadovezem I2 na I1 \n",
    "    m1 = len(I1) # new facilities\n",
    "    m2 = len(I2) #existing facilities\n",
    "\n",
    "\n",
    "    print(\"=======  I   ======\")\n",
    "    \n",
    "    print(\"Ustanove\")\n",
    "    print(\"nove,stare\")\n",
    "    print(I1,I2)\n",
    "    print(\"nove + stare\")\n",
    "    print(I)\n",
    "\n",
    "    print(\"==============\")\n",
    "\n",
    "    \n",
    "\n",
    "\n",
    "    a,b = zip(*tacke[customers])\n",
    "    c,d = zip(*tacke[I1])\n",
    "    e,f  = zip(*tacke[I2])\n",
    "\n",
    "    \n",
    "    ax = plt.gca()\n",
    "    ax.axis('equal')\n",
    "\n",
    "    plt.scatter(a,b,c='red', label = 'Customers')\n",
    "    plt.scatter(e,f,label = 'I2')\n",
    "    plt.scatter(c,d,c='g', label = 'I1')\n",
    "\n",
    "\n",
    "    for i, txt  in enumerate(tacke):\n",
    "        ax.annotate(i,tacke[i])\n",
    "\n",
    "\n",
    "    for i in tacke[facilities]:\n",
    "        circle1 = plt.Circle(tuple(i),radius = 35, color = 'gray', fill = False)\n",
    "\n",
    "        ax.add_artist(circle1)\n",
    "\n",
    "    #plt.legend(loc=1)\n",
    "    #plt.scatter()\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "    #formiram skupove  J1, J2, J = J1 U J2\n",
    "\n",
    "    J1,J2 = podeli_kupce(tacke, radius,customers, I2)\n",
    "    J = J1 + J2\n",
    "\n",
    "    print(\"========  J  =======\")\n",
    "\n",
    "    print(\"Kupci\")\n",
    "    print(customers)\n",
    "    print(\"nisu_oplsluzeni,opsluzeni\")\n",
    "    print(J1, J2)\n",
    "    print(\"nisu_opsluzeni+opsluzeni\")\n",
    "    print(J)\n",
    "\n",
    "\n",
    "    print(\"==================\")\n",
    "\n",
    "\n",
    "    k1 = len(J1) # new facilities\n",
    "    k2 = len(J2) #existing facilities\n",
    "    \n",
    "    \n",
    "    #print(pref(customers,facilities, D))\n",
    "    preference = np.argsort(pref(J,I, D)).T\n",
    "\n",
    "    print('preference \\n',preference)\n",
    "    \n",
    "    print(\"==================\")\n",
    "    \n",
    "    demand = np.random.randint(0,100,len(J))\n",
    "    print(demand)\n",
    "    \n",
    "    \n",
    "    \n",
    "    with open(dat,'w') as f:\n",
    "    \n",
    "        f.write(str(radius))\n",
    "        f.write('\\n')\n",
    "        f.write(str(p))\n",
    "        f.write('\\n')\n",
    "        #f.write(str(t-m))\n",
    "        #f.write('\\n')\n",
    "\n",
    "        f.write(str(k1))\n",
    "        f.write('\\n')\n",
    "        f.write(str(k2))\n",
    "        f.write('\\n')\n",
    "\n",
    "\n",
    "        f.write(str(m1))\n",
    "        f.write('\\n')\n",
    "        f.write(str(m2))\n",
    "        f.write('\\n')\n",
    "\n",
    "\n",
    "        f.write('\\n')\n",
    "        for line in tacke:\n",
    "            f.write(str(line).strip('[]'))\n",
    "            f.write('\\n')\n",
    "\n",
    "\n",
    "        f.write('\\n')\n",
    "        for line in preference:\n",
    "            for item in line:\n",
    "                f.write(str(item)+\" \")\n",
    "        f.write('\\n')\n",
    "\n",
    "        for i in range(len(J)):\n",
    "            f.write(str(demand[i]))\n",
    "            f.write(' ')\n",
    "        f.write('\\n')\n",
    "        f.write('\\n')\n",
    "        ''' f.write(str(customers).replace(', ',\" \").strip('[]'))\n",
    "        f.write('\\n')\n",
    "\n",
    "        '''\n",
    "        \n",
    "        f.write(str(J).replace(', ',\" \").strip('[]'))\n",
    "        f.write('\\n')\n",
    "\n",
    "        #f.write(str(J1).replace(', ',\" \").strip('[]'))\n",
    "        #f.write('\\n')\n",
    "        #f.write(str(J2).replace(', ',\" \").strip('[]').strip(','))\n",
    "        #f.write('\\n')\n",
    "\n",
    "\n",
    "        f.write(str(I).replace(', ',\" \").strip('[]'))\n",
    "        f.write('\\n')\n",
    "\n",
    "        #f.write(str(I1).replace(', ',\" \").strip('[]'))\n",
    "        #f.write('\\n')\n",
    "        #f.write(str(I2).replace(', ',\" \").strip('[]').strip(','))\n",
    "        #f.write('\\n')\n",
    "        \n",
    "\n",
    "    mat = []\n",
    "\n",
    "    print(\"I(j)\")\n",
    "\n",
    "    for j in range(len(customers)):\n",
    "       pom = []\n",
    "       for i in range(len(facilities)):\n",
    "\n",
    "           if rastojanje(tacke[J[j]],tacke[I[i]])<=radius:\n",
    "               pom.append(i)\n",
    "       mat.append(pom)\n",
    "\n",
    "    print(mat)\n",
    "\n",
    "    I_j = mat\n",
    "    mat = []\n",
    "\n",
    "\n",
    "    print(\"J(i)\") \n",
    "    for i in range(len(facilities)):\n",
    "       pom = []\n",
    "       for j in range(len(customers)):\n",
    "\n",
    "           if rastojanje(tacke[J[j]],tacke[I[i]])<=radius:\n",
    "               pom.append(j)\n",
    "       mat.append(pom)\n",
    "\n",
    "    print(mat)\n",
    "    J_i = mat\n",
    "            \n",
    "    return (preference, I_j,J_i)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "'''\n",
    "#SMALL\n",
    "\n",
    "\n",
    "t_niz = [8,8,10,10,12,12,14,14,20,20,25,35,40,40,45,50,80,80,90,100]\n",
    "m_niz = [4,4,4,4,6,6,6,6,11,11,11,11,17,17,17,17,30,30,30,30]\n",
    "r_niz = [50,50,45,45,40,40,38,38,35,35,32,32,30,30,25,25,15,15,12,12]\n",
    "p_niz = [2,3,2,3,4,5,4,5,5,6,5,6,8,9,8,9,20,22,20,22]\n",
    "\n",
    "print(len(t_niz),len(m_niz),len(r_niz),len(p_niz))\n",
    "\n",
    "\n",
    "for i in range(20):\n",
    "    t= t_niz[i]\n",
    "    m = m_niz[i]\n",
    "    p = p_niz[i]\n",
    "    radius = r_niz[i]\n",
    "    datoteka  = 'instance\\mclp_S'+str(i+1)+'.txt'\n",
    "    preference, I_j,J_i = foo(radius,p,t,m,datoteka,deo = 3)\n",
    "'''\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "'''\n",
    "#Medium\n",
    "\n",
    "###PARAMETRI ZA INSTANCE:\n",
    "t_niz = [450,450,450,450,675,675,675,675,900,900,900,900,1350,1350,1350,1350,1800,1800,1800,1800]\n",
    "m_niz = [50,50,50,50,75,75,75,75,100,100,100,100,150,150,150,150,200,200,200,200]\n",
    "r_niz = [80,80,75,75,60,60,55,55,40,40,35,35,30,30,25,25,22,22,20,20]\n",
    "p_niz = [6,7,6,7,10,11,10,11,13,15,13,15,20,22,20,22,27,30,27,30]\n",
    "\n",
    "podaci = zip(t_niz,m_niz,r_niz, p_niz)\n",
    "print(len(t_niz),len(m_niz),len(r_niz),len(p_niz))\n",
    "\n",
    "for i in list(podaci):\n",
    "    print(i)\n",
    "'''\n",
    "##Medium\n",
    "'''\n",
    "for i in range(20):\n",
    "    t= t_niz[i]\n",
    "    m = m_niz[i]\n",
    "    p = p_niz[i]\n",
    "    radius = r_niz[i]\n",
    "    datoteka  = 'instance\\mclp_M'+str(i+1)+'.txt'\n",
    "    preference, I_j,J_i = foo(radius,p,t,m,datoteka)\n",
    "  \n",
    "\n",
    "\n",
    "#pref = np.array(preference).T\n",
    "#pref\n",
    "'''\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "'''\n",
    "#Large\n",
    "\n",
    "###PARAMETRI ZA INSTANCE:\n",
    "t_niz = [1900,1900,1900,1900,2000,2000,2000,2000,2200,2200,2300,2300,2300,2400,2400,2400,2500,2500,3000,3000]\n",
    "m_niz = [100,100,100,100,150,150,150,150,200,200,210,210,200,200,200,220,220,250,250,300]\n",
    "r_niz = [80,80,75,75,60,60,55,55,40,40,35,35,30,30,25,25,22,22,20,20]\n",
    "p_niz = [10,15,10,15,20,25,20,25,22,30,27,33,40,45,40,45,55,60,55,60]\n",
    "\n",
    "podaci = zip(t_niz,m_niz,r_niz, p_niz)\n",
    "print(len(t_niz),len(m_niz),len(r_niz),len(p_niz))\n",
    "\n",
    "for i in list(podaci):\n",
    "    print(i)\n",
    "\n",
    "for i in range(20):\n",
    "    t= t_niz[i]\n",
    "    m = m_niz[i]\n",
    "    p = p_niz[i]\n",
    "    radius = r_niz[i]\n",
    "    datoteka  = 'instance\\mclp_L'+str(i+1)+'.txt'\n",
    "    preference, I_j,J_i = foo(radius,p,t,m,datoteka)\n",
    "  \n",
    "'''\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "'''\n",
    "A = []\n",
    "for j,queue in enumerate(pref):\n",
    "    pom = []\n",
    "    for i,loc in enumerate(queue):\n",
    "        if loc in I_j[j]:\n",
    "            pom.append(loc)\n",
    "    A.append(pom)\n",
    "    \n",
    "            \n",
    "print(A)\n",
    "'''"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "'''#Large\n",
    "\n",
    "###PARAMETRI ZA INSTANCE:\n",
    "t_niz = [4000]\n",
    "m_niz = [300]\n",
    "r_niz = [60]\n",
    "p_niz = [10]\n",
    "\n",
    "podaci = zip(t_niz,m_niz,r_niz, p_niz)\n",
    "print(len(t_niz),len(m_niz),len(r_niz),len(p_niz))\n",
    "\n",
    "for i in list(podaci):\n",
    "    print(i)\n",
    "\n",
    "for i in range(1):\n",
    "    print(i)\n",
    "    t= t_niz[i]\n",
    "    m = m_niz[i]\n",
    "    p = p_niz[i]\n",
    "    radius = r_niz[i]\n",
    "    \n",
    "    datoteka  = 'mclp_XL'+str(i+1)+'.txt'\n",
    "    print(datoteka)\n",
    "    preference, I_j,J_i = foo(radius,p,t,m,datoteka)\n",
    "    print(I_j)\n",
    "'''"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Large\n",
    "\n",
    "###PARAMETRI ZA INSTANCE:\n",
    "t_niz = [2000]\n",
    "m_niz = [200]\n",
    "r_niz = [20]\n",
    "p_niz = [27]\n",
    "\n",
    "podaci = zip(t_niz,m_niz,r_niz, p_niz)\n",
    "print(len(t_niz),len(m_niz),len(r_niz),len(p_niz))\n",
    "\n",
    "for i in list(podaci):\n",
    "    print(i)\n",
    "\n",
    "for i in range(10):\n",
    "    print(i)\n",
    "    t= t_niz[0]\n",
    "    m = m_niz[0]\n",
    "    p = p_niz[0]\n",
    "    radius = r_niz[0]\n",
    "    \n",
    "    datoteka  = 'inst_new/mclp_6.'+str(i+1)+'.txt'\n",
    "    print(datoteka)\n",
    "    preference, I_j,J_i = foo(radius,p,t,m,datoteka)\n",
    "    print(I_j)\n",
    "\n",
    "\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}