Продолжение приложения а

147. if((pu=(int*)calloc(m,sizeof(int)))==NULL) abort();

148. if((pv=(int*)calloc(n,sizeof(int)))==NULL) abort();

149. // ПРИСВОЕНИЕ ВСЕМ ПОТЕНЦИАЛАМ ЗНАЧЕНИЯ MIN

150. for(i=0;i<m;i++)

151. *(pu+i) = MIN;

152. for(j=0;j<n;j++)

153. *(pv+j) = MIN;

154. // Выделение памяти под структуру и заполнение её значениями

155. for(i=0;i<m;i++)

156. {

157. for(j=0;j<n;j++)

158. {

159. if((*(matr+i*n+j) > 0) || (*(matr+i*n+j)==-2))

160. {

161. if((top=(struct poten*)malloc(sizeof(struct poten)))==NULL)

162. abort();

163. fprintf(fil,"top = %d",top);

164. if(!topnast){

165. topnast = top;

166. fprintf(fil,"topnast = top = %d",top);

167. }

168. else top1 -> next=top;

169. top1=top;

170. top -> next=NULL;

171. top -> b = *(st+i*n+j); //Стоимости

172. top -> v = j;

173. top -> u = i;

174. top -> zn = -1;

175. }

176. }

177. }

178. *pu = 0;

179. i=0; j = -1;

180. for(top = topnast;top!=NULL;top = top -> next)

181. {

182. if((top -> u) == i && (top -> v)!=j)

183. {

184. *(pv+(top -> v)) = (top -> b) - *(pu+i);

185. j = (top->v);

186. top -> zn = 0;

187. }

188. else{

189. for(top1 = topnast;top1 != NULL;top1 = top1->next)

190. {

191. if((top1->v) == j && (top1->u)!=i)

192. {

193. *(pu+(top1->u))=(top1->b) - *(pv+j);

194. i = (top1->u);

195. top1 ->zn = 0;

196. break;

Продолжение приложения а

197. }}}}

198. // ********** Продолжение функции подсчета потенциалов *****************

199. for(;;){

200. fl = 0;

201. for(top = topnast;top!=NULL;top =top -> next)

202. {

203. if((top -> zn) == -1)

204. {

205. if(*(pu+(top ->u)) !=MIN)

206. {

207. *(pv+(top->v))=(top->b) - *(pu+(top ->u));

208. fl = 1;

209. top -> zn = 0;

210. }

211. if(*(pv+(top->v)) !=MIN)

212. {

213. *(pu+(top->u))=(top->b) - *(pv+(top->v));

214. fl=1;

215. top->zn = 0;

216. }

217. }

218. }

219. if(!fl) break;

220. }

221. printf("\n ПОТЕНЦИАЛЫ ПО v:");

222. fprintf(fil,"\n **** ПОТЕНЦИАЛЫ ПО v:");

223. for(i = 0;i<n;i++)

224. {

225. printf("%d",*(pv+i));

226. fprintf(fil,"%5d",*(pv+i));

227. }

228. getch();

229. printf("\n ПОТЕНЦИАЛЫ ПО u: ");

230. fprintf(fil,"\n **** ПОТЕНЦИАЛЫ ПО u: ");

231. for(i=0;i<m;i++)

232. {

233. printf("%d",*(pu+i));

234. fprintf(fil,"%5d",*(pu+i));

235. }

236. fprintf(fil,"\n");

237. for(top = topnast;top!=NULL;top = top->next)

238. free(top);

239. return;

240. } // potenzial

241. // ****** PROVERKA PLANA NA OPTIMALNOST' ************************

242. void abcikl(int ik,int jk);

243. int cikl(int ik,int jk);

244. void pr(char pr[],void *st); // Предварительно

Продолжение приложения а

245. int prpoisk(int i1,int j1); // Объявлены

246. int levpoisk(int i1,int j1); //ЭТИ

247. int verpoisk(int i1,int j1); //Функции

248. int nizpoisk(int i1,int j1);

249. int optim(void)

250. {

251. int i,j;

252. for(i=0;i<m;i++)

253. {

254. for(j=0;j<n;j++)

255. {

256. // ИЩЕМ ОПТИМАЛЬНОЕ РЕШЕНИЕ В НАШЕЙ МАТРИЦЕ И ЕСЛИ ЕГО НЕ НАЙДЕМ

257. // ТО ПО СЛЕДУЮЩЕМУ УСЛОВИЮ ПРИСВОИМ ГРАФОКЛЕТКЕ С КООРДИНАТАМИ

258. // ik,jk ЗНАЧЕНИЕ -1

259. if(( *(pu+i)+ *(pv+j))>(*(st+i*n+j))&&((*(matr+i*n+j)) == 0))

260. {

261. abcikl(i,j);

262. fprintf(fil,"optim(): План не оптимален, функции main() возвращаем -1,\n а abcikl() параметры i,j ");

263. return(-1);

264. }

265. }

266. }

267. fprintf(fil,"Plan optimalen");

268. return(0);

269. } // ******* optim() ***************

270. // ************** UPGRADE PLAN **************************

271. void abcikl(int ik,int jk)

272. {

273. int i,j;

274. fprintf(fil,"Мы в abcikl()");

275. if((matr2=(int*)calloc(n*m,sizeof(int))) == NULL) abort();

276. for(i=0;i<m;i++)

277. {

278. for(j=0;j<n;j++)

279. {

280. *(matr2+i*n+j) = *(matr+i*n+j); // Создаем копию рабочей матрицы

281. }

282. }

283. *(matr2+ik*n+jk) = -1;

284. // значению матрицы с параметрами ik,jk присваеваем -1

285. printf("\n");

286. printf("Поиск Цепи: \n\n");

287. fprintf(fil,"Поиск Цепи: \n\n");

288. for(i=0;i<m;i++)

289. {

290. for(j=0;j<n;j++)

291. {