求结构可靠性分析的分支限界法的程序

helenyuan

请问各位，有没有人有结构可靠性分析中寻找失效模式的分支限界法的程序啊，急用！
多谢大家的帮忙！

风花雪月

我有几个分支限界法的应用程序，不过不是用于可靠性的，属于应用范例

比如用分支限界法解决0/1背包问题、最小圆排列的分支限界法

helenyuan

能否传给我，也许有参考价值
多谢了
我的邮箱　king101010101010@163.com

风花雪月

分支限界法中的单源最短路径问题的实现，是用VC＋＋编写的

风花雪月

最小圆排列的分支限界法！

1. #include <queue>
   
2. #include <fstream>
   
3. #include<iostream>
   
4. #include <math.h>
   
5. #include <time.h>
   
6. using namespace std;
   
7. double **dic;//dic[i][j]存放第i个圆与第j个圆得圆心
   
8. double *r;//存放n个圆的半径
   
9. double best;
   
10. 
    
11. 
    
12. template<class type>
    
13. int partition (type a[],int low,int high){
    
14.         int pivotpos=low;
    
15.         type pivot=a[low];
    
16.     for (int i=low+1;i<=high;i++)
    
17.                 if (a[i]<pivot&&++pivotpos!=i)
    
18.                         swap(a[pivotpos],a[i]);
    
19.                 swap(a[low],a[pivotpos]);
    
20.                 return pivotpos;}
    
21. template<class type>
    
22. void quicksort (type a[],int p,int r){
    
23.         int i;
    
24.         if(p<r)
    
25.         {i=partition(a,p,r);
    
26.         quicksort(a,p,i-1);
    
27.         quicksort(a,i+1,r);}}
    
28. class circlenode
    
29. {
    
30.         friend void circleperm(double *r,int n);
    
31.        
    
32. private:
    
33.         int mink;//代排的圆中第mink个圆的半径最小
    
34.         int s;//算法完成了s步，即排好了s个圆
    
35.         int k;//镜像剪枝
    
36.         double *x;//圆心坐标
    
37.         int *rp;//所选的第s个圆得半径为r[rp[s]]
    
38.         double compute(int n);
    
39.         double center(int t);
    
40.        
    
41. };
    
42. 
    
43. 
    
44. 
    
45. 
    
46. double circlenode::compute(int n)
    
47. {
    
48.         int i;
    
49.         double low=0.0;
    
50.         double high=x[n-1]+r[rp[n-1]];
    
51.         for(i=0;i<n;i++)
    
52.         {
    
53.                 if(x[i]-r[rp[i]]<low)
    
54.                         low=(double)x[i]-r[rp[i]];
    
55.                 if(x[i]+r[rp[i]]>high)
    
56.                         high=(double)x[i]+r[rp[i]];
    
57.         }
    
58.         return (double)high-low;
    
59. };
    
60. 
    
61. 
    
62. 
    
63. double circlenode::center(int t)
    
64. {
    
65.         double temp=0.0;
    
66.         double valuex;
    
67.         for(int i=0;i<t;i++)
    
68.         {
    
69.                 valuex=x[i]+(double)2*sqrt(r[rp[i]]*r[rp[t]]);
    
70.                 if(valuex>temp)
    
71.                         temp=valuex;
    
72.         }
    
73.         return temp;
    
74. };
    
75. 
    
76. 
    
77. void circleperm(double *r,int n)
    
78. {
    
79.         int i,j;
    
80.         int wk;
    
81.         int tag;
    
82.         int tt;
    
83.         double bestt;
    
84.         double *tsame;
    
85.         int tsamen;
    
86.         double rtemp;
    
87.         double mintemp;
    
88.         tsame=new double[n];
    
89.     queue<circlenode> qu;
    
90.        
    
91.         /*e.x=new double[n];
    
92.     e.s=-1;
    
93.         e.mink=0;
    
94.     e.rp=new int [n];
    
95.         for(i=0;i<n;i++)
    
96.         e.rp[i]=i;*/
    
97.         circlenode tempe;
    
98.         tempe.x=new double[n];
    
99.         tempe.rp=new int[n];
    
100.         int fi=0;
     
101.         int li=n-1;
     
102.         i=0;
     
103.         while(fi<=li)
     
104.         {
     
105.                 tempe.rp[i]=fi;
     
106.                 i++;
     
107.                 fi++;
     
108.                 if(fi<=li)
     
109.                 {
     
110.                         tempe.rp[i]=li;
     
111.                         i++;
     
112.                         li--;
     
113.                 }
     
114.         }
     
115.         for(i=0;i<n;i++)
     
116.                 tempe.x[i]=tempe.center(i);
     
117.         best=tempe.compute(n);//算初值
     
118.         /* for(i=0;i<n;i++)
     
119.         cout<<tempe.rp[i]<<" ";
     
120.          cout<<endl;
     
121.          
     
122.          
     
123.       for(i=0;i<n;i++)
     
124.         {cout<<tempe.x[i]<<endl;}//test best*/
     
125.        
     
126.         delete []tempe.x;
     
127.         delete []tempe.rp;
     
128.        
     
129.         cout<<best<<endl;
     
130.        
     
131.         // best=12345;//需修改
     
132.         circlenode e;
     
133.         for(i=0;i<n-1;i++)
     
134.         {
     
135.                
     
136.                
     
137.                
     
138.                 //cout<<r[i]<<endl;
     
139.                 if(i>0&&r[i-1]==r[i]) continue;
     
140.                 if(i==0) e.mink=1;
     
141.                 else e.mink=0;
     
142.                 e.k=n-i-1;//比i大的数字的个数
     
143.                 e.rp=new int[n];
     
144.                 e.x=new double[n];
     
145.                
     
146.                 for(j=0;j<n;j++)
     
147.                 {
     
148.                         e.rp[j]=j;
     
149.                        
     
150.                 };
     
151.                 e.x[0]=0;
     
152.                 e.rp[0]=i;
     
153.                 e.rp[i]=0;
     
154.                
     
155.                 e.s=0;
     
156.                 qu.push(e);
     
157.         }
     
158.         while(!qu.empty())
     
159.         {
     
160.                 e=qu.front();
     
161.                 qu.pop();
     
162.                 if (e.s==n-3)
     
163.                 {
     
164.                         if(e.rp[n-1]>e.rp[0])
     
165.                         {
     
166.                                 e.x[n-2]=e.center(n-2);
     
167.                                 e.x[n-1]=e.center(n-1);
     
168.                                 bestt=e.compute(n);
     
169.                                
     
170.                                
     
171.                                
     
172.                                
     
173.                                 if(bestt<best)
     
174.                                         best=bestt;
     
175.                         }
     
176.                         if(e.rp[n-2]>e.rp[0])
     
177.                         {
     
178.                                 tt=e.rp[n-2];
     
179.                                 e.rp[n-2]=e.rp[n-1];
     
180.                                 e.rp[n-1]=tt;
     
181.                                
     
182.                                
     
183.                                 e.x[n-2]=e.center(n-2);
     
184.                                 e.x[n-1]=e.center(n-1);
     
185.                                 bestt=e.compute(n);
     
186.                                
     
187.                                
     
188.                                 if(bestt<best)
     
189.                                         best=bestt;
     
190.                         }
     
191.                 }
     
192.                 else
     
193.                 {
     
194.                        
     
195.                         tsamen=0;
     
196.                        
     
197.                         for(i=e.s+1;i<n;i++)
     
198.                         {
     
199.                                 if(e.rp[i]>e.rp[0]) wk=e.k-1;
     
200.                                 else wk=e.k;
     
201.                                
     
202.                                 if(wk) //镜像剪枝
     
203.                                 {
     
204.                                        
     
205.                                         tag=0;
     
206.                                         rtemp=r[e.rp[i]];
     
207.                                         for(j=0;j<tsamen;j++)
     
208.                                         {
     
209.                                                 if(rtemp==tsame[j])
     
210.                                                         tag=1;
     
211.                                                 //        continue;
     
212.                                         }
     
213.                                         if(!tag)
     
214.                                         {
     
215.                                                
     
216.                                                 tsame[tsamen]=rtemp;
     
217.                                                 tsamen++;
     
218.                                                 circlenode w;
     
219.                                                 w.k=wk;
     
220.                                                 w.s=e.s+1;
     
221.                                                 w.x=new double[n];
     
222.                                                 w.rp=new int[n];
     
223.                                                
     
224.                                                 for(j=0;j<n;j++)
     
225.                                                 {
     
226.                                                         w.x[j]=e.x[j];
     
227.                                                         w.rp[j]=e.rp[j];
     
228.                                                 }
     
229.                                                
     
230.                                                
     
231.                                                 w.rp[w.s]=e.rp[i];
     
232.                                                 w.rp[i]=e.rp[w.s];
     
233.                                                
     
234.                                                 w.mink=e.mink;
     
235.                                                 if(w.mink==w.rp[w.s])
     
236.                                                 {
     
237.                                                         w.mink=w.rp[w.s+1];
     
238.                                                         for(j=w.s+2;j<n;j++)
     
239.                                                         {
     
240.                                                                 if(r[w.rp[j]]<r[w.mink])
     
241.                                                                         w.mink=w.rp[j];
     
242.                                                         }
     
243.                                                        
     
244.                                                 }
     
245.                                                 w.x[w.s]=w.center(w.s);
     
246.                                                 mintemp=w.x[w.s]+(2*n-2*w.s-1)*r[w.mink]+r[w.rp[0]];
     
247.                                                 if(mintemp<best)
     
248.                                                 {
     
249.                                                         qu.push(w);
     
250.                                                 }
     
251.                                                 else
     
252.                                                 {
     
253.                                                         delete []w.x;
     
254.                                                         delete []w.rp;
     
255.                                                 }
     
256.                                         }  //if(!tag)
     
257.                                 }//if(wk)
     
258.                         }//for(i=e.s+1;i<n;i++)
     
259.                        
     
260.                 }//else
     
261.                 delete []e.x;
     
262.                 delete []e.rp;       
     
263.         }//while
     
264.         delete []tsame;
     
265.        
     
266. }
     
267. void main()
     
268. {
     
269.         clock_t start, finish;
     
270.         start=clock();
     
271.         fstream infile,outfile;
     
272.         int n;
     
273.         int i;
     
274.         int j;
     
275.        
     
276.         infile.open("input.txt",ios::in);
     
277.         outfile.open("output.txt",ios.out);
     
278.         infile>>n;
     
279.         r=new double[n];
     
280.         dic=new double*[n];
     
281.         for(i=0;i<n;i++)
     
282.                 dic[i]=new double[n];
     
283.         for(i=0;i<n;i++)
     
284.                 infile>>r[i];
     
285.        
     
286.     if(n==1)
     
287.         {
     
288.                 best=2.0*r[0];
     
289.                 outfile<<best<<endl;
     
290.                 cout<<best<<endl;
     
291.         }
     
292.         else if(n==2)
     
293.         {
     
294.                 double temp2;
     
295.                 temp2=2*sqrt(r[0]*r[1]);
     
296.                 double low2,high2;
     
297.                 low2=0-r[0];
     
298.                 if(temp2-r[1]<low2) low2=temp2-r[1];
     
299.                 high2=temp2+r[1];
     
300.                 if(high2<r[0]) high2=r[0];
     
301.                 best=high2-low2;
     
302.                 outfile<<best<<endl;
     
303.                 cout<<best<<endl;
     
304.         }
     
305.         else
     
306.         {
     
307.                 quicksort(r,0,n-1);
     
308.                 /*        for(i=0;i<n;i++)
     
309.                 cout<<r[i]<<endl;*/
     
310.                 for(i=0;i<n-1;i++)
     
311.                 {
     
312.                         for(j=i+1;j<n;j++)
     
313.                                 dic[i][j]=dic[j][i]=sqrt(r[i]*r[j]);
     
314.                 }
     
315.                
     
316.                 circleperm(r,n);
     
317.                 outfile<<best<<endl;
     
318.                 cout<<best<<endl;
     
319.         }
     
320.         infile.close();
     
321.         outfile.close();
     
322.         delete []r;
     
323.         for(i=0;i<n;i++)
     
324.                 delete []dic[i];
     
325.         delete []dic;
     
326.        
     
327.         finish=clock();
     
328.         cout<<endl<<"Elapsed Time: "<<(double)(finish-start)/CLOCKS_PER_SEC<<" secouds"<<endl;
     
329.        
     
330. }       
     

复制代码

helenyuan

多些，会好好学习一下的

cat

很好!:@)

jinbin

感谢分享！

luther_5015133

谢谢了，正在学习中

lq12131010

不要学了！30年前的东西了。。。