結果

問題 No.317 辺の追加
ユーザー anta
提出日時 2015-12-10 19:49:21
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
CE  
(最新)
AC  
(最初)
実行時間 -
コード長 8,875 bytes
コンパイル時間 253 ms
コンパイル使用メモリ 41,600 KB
最終ジャッジ日時 2024-11-14 19:30:36
合計ジャッジ時間 926 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
コンパイルエラー時のメッセージ・ソースコードは、提出者また管理者しか表示できないようにしております。（リジャッジ後のコンパイルエラーは公開されます）
ただし、clay言語の場合は開発者のデバッグのため、公開されます。

コンパイルメッセージ
main.cpp:46:31: error: ‘std::string’ has not been declared
   46 |         void read_string(std::string *str) {
      |                               ^~~~~~
main.cpp:52:29: error: ‘std::string’ has not been declared
   52 |         void read_line(std::string *str) {
      |                             ^~~~~~
main.cpp:79:36: error: ‘std::string’ has not been declared
   79 |         FastInput &operator()(std::string &res) { read_string(&res); return *this; }
      |                                    ^~~~~~
main.cpp:79:20: error: ‘FastInput& FastInput::operator()(int&)’ cannot be overloaded with ‘FastInput& FastInput::operator()(int&)’
   79 |         FastInput &operator()(std::string &res) { read_string(&res); return *this; }
      |                    ^~~~~~~~
main.cpp:74:20: note: previous declaration ‘FastInput& FastInput::operator()(int&)’
   74 |         FastInput &operator()(int &res) { read_signed(&res); return *this; }
      |                    ^~~~~~~~
main.cpp: In member function ‘void FastInput::read_string(int*)’:
main.cpp:47:22: error: request for member ‘clear’ in ‘* str’, which is of non-class type ‘int’
   47 |                 str->clear();
      |                      ^~~~~
main.cpp: In member function ‘void FastInput::read_line(int*)’:
main.cpp:53:22: error: request for member ‘clear’ in ‘* str’, which is of non-class type ‘int’
   53 |                 str->clear();
      |                      ^~~~~
main.cpp:56:26: error: request for member ‘empty’ in ‘* str’, which is of non-class type ‘int’
   56 |                 if(!str->empty() && (*str)[str->size() - 1] == '\r') str->resize(str->size() - 1);
      |                          ^~~~~
main.cpp:56:49: error: request for member ‘size’ in ‘* str’, which is of non-class type ‘int’
   56 |                 if(!str->empty() && (*str)[str->size() - 1] == '\r') str->resize(str->size() - 1);
      |

ソースコード

diff # 

#include <cstdio>
#include <vector>
#include <algorithm>
#include <cassert>
#define rep(i,n) for(int (i)=0;(i)<(int)(n);++(i))
#ifndef __GNUC__
#include <intrin.h>
#endif
using namespace std;


class FastInput {
	bool _end;
public:
	FastInput() : _end(false) { }

	operator void*() { return _end ? 0 : (void*)this; }

	template<typename T>
	void read_unsigned(T *res) {
		T x = 0;
		for(char c = skip(); '0' <= c && c <= '9'; c = gc())
			x = x * 10 + (c - '0');
		*res = x;
	}

	template<typename T>
	void read_signed(T *res) {
		char c = skip();
		bool sign = false;
		if(c == '-') sign = true, c = gc();
		T x = 0;
		for(; '0' <= c && c <= '9'; c = gc())
			x = x * 10 + (c - '0');
		*res = !sign ? x : -x;
	}

	void read_c_string(char *str, int *len) {
		int n = 0;
		for(char c = skip(); !is_delim(c); c = gc())
			str[n ++] = c;
		str[n] = 0;
		*len = n;
	}

	void read_string(std::string *str) {
		str->clear();
		for(char c = skip(); !is_delim(c); c = gc())
			*str += c;
	}

	void read_line(std::string *str) {
		str->clear();
		for(char c = gc(); c != '\n'; c = gc())
			*str += c;
		if(!str->empty() && (*str)[str->size() - 1] == '\r') str->resize(str->size() - 1);
	}

	void read_double(double *res) {
		std::string buf;
		read_string(&buf);
		sscanf(buf.c_str(), "%lf", res);
	}

	void read_char(char *res) {
		*res = skip();
	}

	void read_string_buf(char *res, int n) {
		*res = skip();
	}

	FastInput &operator()(char &res) { read_char(&res); return *this; }
	FastInput &operator()(int &res) { read_signed(&res); return *this; }
	FastInput &operator()(unsigned &res) { read_unsigned(&res); return *this; }
	FastInput &operator()(long long &res) { read_signed(&res); return *this; }
	FastInput &operator()(unsigned long long &res) { read_unsigned(&res); return *this; }
	FastInput &operator()(char *res) { int len; read_c_string(res, &len); return *this; }
	FastInput &operator()(std::string &res) { read_string(&res); return *this; }
	FastInput &operator()(double &res) { read_double(&res); return *this; }

	template<typename T1, typename T2>
	FastInput &operator()(T1 &res1, T2 &res2) {
		return operator()(res1)(res2);
	}
	template<typename T1, typename T2, typename T3>
	FastInput &operator()(T1 &res1, T2 &res2, T3 &res3) {
		return operator()(res1)(res2)(res3);
	}

	template<typename T>
	FastInput &a(T *a, int n) {
		for(int i = 0; i < n; ++ i)
			operator()(a[i]);
		return *this;
	}

	template<typename T>
	FastInput &operator()(vector<T> &v) {
		for(size_t i = 0; i < v.size(); ++ i)
			operator()(v[i]);
		return *this;
	}

private:
	static char gc() {
#if defined(__GNUC__) && !defined(__MINGW32__)
		return (char)getchar_unlocked();
#elif defined(_MSC_VER)
		return (char)_getchar_nolock();
#else
		return (char)getchar();
#endif
	}

	static bool is_delim(char c) {
		return c == ' ' || c == '\t' || c == '\n' || c == '\r' || c == EOF;
	}

	char skip() {
		if(_end) return EOF;
		char c;
		for(c = gc(); c != -1 && is_delim(c); c = gc()) ;
		if(c == EOF) _end = true;
		return c;
	}
} in;

class FastOutput {
public:

	template<typename T>
	void print_unsigned(T x) {
		char buf[24]; int n = 0;
		do buf[n ++] = x % 10, x /= 10; while(x != 0);
		while(n > 0) pc('0' + buf[-- n]);
	}

	template<typename T>
	void print_signed(T x) {
		char buf[24]; int n = 0; bool sign = false;
		if(x < 0) sign = true;
		//depends on the -x % y behaviour
		do buf[n ++] = x % 10, x /= 10; while(x != 0);
		if(!sign) {
			while(n > 0) pc('0' + buf[-- n]);
		} else {
			pc('-');
			while(n > 0) pc('0' - buf[-- n]);
		}
	}

	void print_c_string(const char *str) {
		for(const char *p = str; *p; ++ p)
			pc(*p);
	}

	void print_string(const std::string &str) {
		print_c_string(str.c_str());
	}

	void print_double(double x, int digits) {
		//|res| < 2^1024 ~ 10^308
		char buf[512];
#ifdef _MSC_VER
		sprintf_s(buf, "%.*f", digits, x);
#else
		sprintf(buf, "%.*f", digits, x);
#endif
		print_c_string(buf);
	}

	void print_char(char x) {
		pc(x);
	}

	FastOutput &operator()(char x) { print_char(x); return *this; }
	FastOutput &operator()(int x) { print_signed(x); return *this; }
	FastOutput &operator()(unsigned x) { print_unsigned(x); return *this; }
	FastOutput &operator()(long long x) { print_signed(x); return *this; }
	FastOutput &operator()(unsigned long long x) { print_unsigned(x); return *this; }
	FastOutput &operator()(const char *str) { print_c_string(str); return *this; }
	FastOutput &operator()(const std::string &str) { print_string(str); return *this; }
	FastOutput &operator()(double x) { print_double(x, 10); return *this; }

	template<typename T1, typename T2>
	FastOutput &operator()(T1 x1, T2 x2) { return operator()(x1)(x2); }
	template<typename T1, typename T2, typename T3>
	FastOutput &operator()(T1 x1, T2 x2, T3 x3) { return operator()(x1)(x2)(x3); }
	template<typename T1, typename T2, typename T3, typename T4>
	FastOutput &operator()(T1 x1, T2 x2, T3 x3, T4 x4) { return operator()(x1)(x2)(x3)(x4); }

private:
	static void pc(char c) {
#if defined(__GNUC__) && !defined(__MINGW32__)
		putchar_unlocked(c);
#elif defined(_MSC_VER)
		_putchar_nolock(c);
#else
		putchar(c);
#endif
	}
} out;

class FastOutputN {
public:
	template<typename T1>
	FastOutputN &operator()(T1 x1) { out(x1)('\n'); return *this; }
	template<typename T1, typename T2>
	FastOutputN &operator()(T1 x1, T2 x2) { out(x1, x2)('\n'); return *this; }
	template<typename T1, typename T2, typename T3>
	FastOutputN &operator()(T1 x1, T2 x2, T3 x3) { out(x1, x2, x3)('\n'); return *this; }
	template<typename T1, typename T2, typename T3, typename T4>
	FastOutputN &operator()(T1 x1, T2 x2, T3 x3, T4 x4) { out(x1, x2, x3, x4)('\n'); return *this; }

	template<typename T>
	FastOutputN &a(const T *a, int n) {
		for(int i = 0; i < n; ++ i) {
			if(i != 0) out(' ');
			out(a[i]);
		}
		out('\n');
		return *this;
	}

	template<typename T>
	FastOutputN &operator()(const vector<T> &v) {
		for(size_t i = 0; i < v.size(); ++ i) {
			if(i != 0) out(' ');
			out(v[i]);
		}
		out('\n');
		return *this;
	}
} outn;


struct UnionFind {
	vector<int> data;
	void init(int n) { data.assign(n, -1); }
	bool unionSet(int x, int y) {
		x = root(x); y = root(y);
		if(x != y) {
			if(data[y] < data[x]) swap(x, y);
			data[x] += data[y]; data[y] = x;
		}
		return x != y;
	}
	bool findSet(int x, int y) { return root(x) == root(y); }
	int root(int x) { return data[x] < 0 ? x : data[x] = root(data[x]); }
	int size(int x) { return -data[root(x)]; }
};

void update(unsigned short *dp, int n, int s, int t) {
	assert(s > 1);
	for(; n >= 0 && n % 8 != 7; -- n) {
		unsigned short &x = dp[n + s], y = dp[n] + t;
		if(y >= t && x > y) x = y;
	}
#ifndef __GNUC__
	const __m128i tttttttt = _mm_set1_epi16(t);
	const __m128i infs = _mm_set1_epi32(0xffffffffU);
	for(int i = n - 7; i >= 0; i -= 8) {
		const __m128i x = _mm_loadu_si128(reinterpret_cast<const __m128i*>(dp + (i + s)));
		const __m128i a = _mm_loadu_si128(reinterpret_cast<const __m128i*>(dp + i));
		const __m128i b = _mm_add_epi16(a, tttttttt);
		const __m128i c = _mm_cmpeq_epi16(a, infs);
		const __m128i y = _mm_or_si128(b, c);
		const __m128i z = _mm_min_epu16(x, y);
		_mm_storeu_si128(reinterpret_cast<__m128i*>(dp + (i + s)), z);
	}
#else
	__attribute__((aligned(16))) const unsigned char pat[16] = { 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1 };
	__asm("vmovd %0, %%xmm2\n\tvpshufb %1, %%xmm2, %%xmm2":: "r"(t), "m"(pat) : "xmm2");
	__asm("vpcmpeqd %%xmm3, %%xmm3, %%xmm3"::: "xmm3");
	for(int i = n - 7; i >= 0; i -= 8) {
		const unsigned short * const p = dp + i;
		unsigned short * const q = dp + (i + s);
		__asm("vmovdqu %1, %%xmm0\n\tvpaddw %%xmm0, %%xmm2, %%xmm1\n\tvpcmpeqw %%xmm3, %%xmm0, %%xmm0\n\tvpor %%xmm0, %%xmm1, %%xmm0\n\tvpminuw %0, %%xmm0, %%xmm0\n\tvmovups %%xmm0, %0"
			: "=m"(*q) : "m"(*p) : "xmm0", "xmm1", "xmm2", "xmm3");
	}
#endif
}


int main() {
	int N; int M;
	while(in(N, M)) {
		UnionFind uf; uf.init(N);
		rep(i, M) {
			int u; int v;
			in(u, v), -- u, -- v;
			uf.unionSet(u, v);
		}
		vector<int> counts(N + 1);
		rep(i, N) if(uf.root(i) == i)
			++ counts[uf.size(i)];
		vector<pair<int, int> > v;
		for(int i = 2; i <= N; ++ i) {
			int n = counts[i];
			while(n > 0) {
				int m = (n + 1) / 2;
				v.emplace_back(i * m, m);
				n -= m;
			}
		}
		const unsigned short INFS = 0xffff;
		vector<unsigned short> dp(N + 1, INFS);
		sort(v.begin(), v.end());
		int sum = 0;
		dp[0] = 0;
		for(auto p : v) {
			update(&dp[0], sum, p.first, p.second);
			sum += p.first;
		}
		int n = counts[1];
		vector<pair<int, int> > q(N + 1);
		int h = 0, t = 0;
		for(int i = 0; i <= N; ++ i) {
			while(h < t && q[h].first < i - n)
				++ h;
			if(dp[i] != INFS) {
				int y = dp[i] - i;
				while(h < t && q[t - 1].second >= y)
					-- t;
				q[t ++] = {i, y};
			}
			if(1 <= i) {
				if(h == t)
					outn("-1");
				else
					outn(q[h].second + i - 1);
			}
		}
	}
	return 0;
}
0