// ===================================================================================== // // Filename: test_dimacs.cpp // // Description: test if the instance contain loop or multi edge // // Version: 1.0 // Created: 2017年 12月 28日 星期四 06:27:25 -00 // Revision: none // Compiler: g++ // // Author: Jinkun Lin, jkunlin@gmail.com // Organization: School of EECS, Peking University // // ===================================================================================== #include #include #include #include using namespace std; vector> adjacency_list; vector vertex_weight; // read graph file into adjacency_list void fast_read(string file_name) {/*{{{*/ ifstream in_file(file_name); if (!in_file) { cout << "in_file error" << endl; exit(1); } in_file.seekg (0, in_file.end); size_t file_len = in_file.tellg(); in_file.seekg (0, in_file.beg); char *data = new char [file_len]; in_file.read(data, file_len); in_file.close(); //skip comments char *pos = data; while (*pos == '%') { while(*(pos++) != '\n'); } //read vertex_count int vertex_count = 0, edge_count = 0; while (*pos < '0' || *pos > '9') { ++pos; } while (*pos != ' ') { vertex_count = vertex_count * 10 + *pos - '0'; ++pos; } //read edge_count while (*pos < '0' || *pos > '9') { ++pos; } while (*pos >= '0' && *pos <= '9') { edge_count = edge_count * 10 + *pos - '0'; ++pos; } //read vertex_weight vertex_weight.resize(vertex_count + 1); int v,w; for(vector>::size_type i=1; i '9') { ++pos; } while (*pos != ' ') { v = v * 10 + *pos - '0'; ++pos; } //read weight while (*pos < '0' || *pos > '9') { ++pos; } while (*pos >= '0' && *pos <= '9') { w = w * 10 + *pos - '0'; ++pos; } vertex_weight[v]=w; } //read adjacency_list adjacency_list.resize(vertex_count + 1); vector vertex_degree(vertex_count + 1, 0); char *stash_pos = pos; int v1, v2; for (int i = 0; i < edge_count; ++i) { v1 = v2 = 0; //read v1 while (*pos < '0' || *pos > '9') { ++pos; } while (*pos != ' ') { v1 = v1 * 10 + *pos - '0'; ++pos; } //read v2 while (*pos < '0' || *pos > '9') { ++pos; } while (*pos >= '0' && *pos <= '9') { v2 = v2 * 10 + *pos - '0'; ++pos; } vertex_degree[v1]++; vertex_degree[v2]++; } for (size_t v = 1; v < adjacency_list.size(); ++v) { adjacency_list[v].reserve(vertex_degree[v]); } pos = stash_pos; for (int i = 0; i < edge_count; ++i) { v1 = v2 = 0; //read v1 while (*pos < '0' || *pos > '9') { ++pos; } while (*pos != ' ') { v1 = v1 * 10 + *pos - '0'; ++pos; } //read weight while (*pos < '0' || *pos > '9') { ++pos; } while (*pos >= '0' && *pos <= '9') { v2 = v2 * 10 + *pos - '0'; ++pos; } adjacency_list[v1].push_back(v2); adjacency_list[v2].push_back(v1); } delete[] data; }/*}}}*/ void test() { for (size_t v = 1; v < adjacency_list.size(); ++v) { std::sort(adjacency_list[v].begin(), adjacency_list[v].end()); for (size_t i = 0; i < adjacency_list[v].size(); ++i) { if (adjacency_list[v][i] == (int)v) { std::cout << "loop" << std::endl; abort(); } if (i > 0 && adjacency_list[v][i] == adjacency_list[v][i - 1]) { std::cout << "multi edge" << std::endl; abort(); } } } } int main(int argc, char const *argv[]) { if (argc != 2) { std::cout << "usage" << std::endl; return 1; } fast_read(argv[1]); test(); return 0; }