| [699] | 1 | #include<math.h> | 
|---|
| [921] | 2 | #include<lemon/list_graph.h> | 
|---|
| [711] | 3 |  | 
|---|
|  | 4 | #include"bench_tools.h" | 
|---|
| [699] | 5 |  | 
|---|
| [921] | 6 | using namespace lemon; | 
|---|
| [699] | 7 |  | 
|---|
|  | 8 | ///Makes a full graph by adding and deleting a lot of edges; | 
|---|
|  | 9 |  | 
|---|
|  | 10 | ///\param n Number of nodes. | 
|---|
|  | 11 | ///\param rat The funcion will make \f$rat\timesn^2\f$ edge addition and | 
|---|
|  | 12 | ///\f$(rat-1)\timesn^2\f$ deletion. | 
|---|
|  | 13 | ///\param p Tuning parameters. | 
|---|
|  | 14 | ///\warning \c rat, \c p, and \c n must be pairwise relative primes. | 
|---|
|  | 15 | template <class Graph> | 
|---|
|  | 16 | void makeFullGraph(int n, int rat, int p) | 
|---|
|  | 17 | { | 
|---|
|  | 18 | GRAPH_TYPEDEF_FACTORY(Graph); | 
|---|
|  | 19 |  | 
|---|
|  | 20 | Graph G; | 
|---|
|  | 21 |  | 
|---|
| [708] | 22 | //  Node nodes[n]; | 
|---|
|  | 23 | std::vector<Node> nodes(n); | 
|---|
| [699] | 24 | for(int i=0;i<n;i++) nodes[i]=G.addNode(); | 
|---|
|  | 25 |  | 
|---|
| [708] | 26 | //Edge equ[rat]; | 
|---|
|  | 27 | std::vector<Edge> equ(rat); | 
|---|
| [699] | 28 |  | 
|---|
| [718] | 29 | long long int count; | 
|---|
| [699] | 30 |  | 
|---|
|  | 31 | for(count=0;count<rat;count++) { | 
|---|
|  | 32 | equ[count%rat]=G.addEdge(nodes[(count*p)%n],nodes[(count*p/n)%n]); | 
|---|
|  | 33 | } | 
|---|
|  | 34 | for(;(count%rat)||((count*p)%n)||((count*p/n)%n);count++) { | 
|---|
|  | 35 | //    if(!(count%1000000)) fprintf(stderr,"%d\r",count); | 
|---|
|  | 36 | if(count%rat) G.erase(equ[count%rat]); | 
|---|
|  | 37 | equ[count%rat]=G.addEdge(nodes[(count*p)%n],nodes[(count*p/n)%n]); | 
|---|
|  | 38 | } | 
|---|
| [718] | 39 | //   std::cout << "Added " << count | 
|---|
|  | 40 | //          << " ( " << n << "^2 * " << rat << " ) edges\n"; | 
|---|
|  | 41 |  | 
|---|
|  | 42 |  | 
|---|
| [699] | 43 | //  for(int i=0;1;i++) ; | 
|---|
|  | 44 | } | 
|---|
|  | 45 |  | 
|---|
|  | 46 | int main() | 
|---|
|  | 47 | { | 
|---|
| [921] | 48 | lemon::Timer T; | 
|---|
| [699] | 49 | makeFullGraph<ListGraph>(nextPrim(1000),nextPrim(300),nextPrim(100)); | 
|---|
| [718] | 50 |  | 
|---|
|  | 51 | PrintTime("BIG",T); | 
|---|
| [699] | 52 | T.reset(); | 
|---|
|  | 53 | makeFullGraph<ListGraph>(nextPrim(100),nextPrim(30000),nextPrim(150)); | 
|---|
| [718] | 54 |  | 
|---|
|  | 55 | PrintTime("SMALL",T); | 
|---|
| [699] | 56 | } | 
|---|