結果

問題 No.317 辺の追加
ユーザー LayCurseLayCurse
提出日時 2015-12-10 17:04:41
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 22 ms / 2,000 ms
コード長 7,779 bytes
コンパイル時間 1,510 ms
コンパイル使用メモリ 164,520 KB
実行使用メモリ 8,064 KB
最終ジャッジ日時 2024-09-15 07:38:07
合計ジャッジ時間 4,315 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 13 ms
6,016 KB
testcase_03 AC 15 ms
6,528 KB
testcase_04 AC 11 ms
6,016 KB
testcase_05 AC 15 ms
6,528 KB
testcase_06 AC 19 ms
7,168 KB
testcase_07 AC 7 ms
5,376 KB
testcase_08 AC 13 ms
6,400 KB
testcase_09 AC 17 ms
6,784 KB
testcase_10 AC 22 ms
7,808 KB
testcase_11 AC 7 ms
5,760 KB
testcase_12 AC 21 ms
7,936 KB
testcase_13 AC 9 ms
5,888 KB
testcase_14 AC 19 ms
7,040 KB
testcase_15 AC 21 ms
7,552 KB
testcase_16 AC 14 ms
6,528 KB
testcase_17 AC 20 ms
7,424 KB
testcase_18 AC 22 ms
7,808 KB
testcase_19 AC 22 ms
8,064 KB
testcase_20 AC 14 ms
6,400 KB
testcase_21 AC 6 ms
5,376 KB
testcase_22 AC 4 ms
5,376 KB
testcase_23 AC 4 ms
5,376 KB
testcase_24 AC 9 ms
6,144 KB
testcase_25 AC 7 ms
5,504 KB
testcase_26 AC 7 ms
5,376 KB
testcase_27 AC 6 ms
5,376 KB
testcase_28 AC 4 ms
5,376 KB
testcase_29 AC 2 ms
5,376 KB
testcase_30 AC 14 ms
6,528 KB
testcase_31 AC 13 ms
6,528 KB
testcase_32 AC 13 ms
6,656 KB
testcase_33 AC 13 ms
6,656 KB
testcase_34 AC 13 ms
6,528 KB
testcase_35 AC 13 ms
6,528 KB
testcase_36 AC 13 ms
6,528 KB
testcase_37 AC 13 ms
6,528 KB
testcase_38 AC 14 ms
6,528 KB
testcase_39 AC 13 ms
6,528 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘void reader(double*)’:
main.cpp:16:29: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   16 | void reader(double *x){scanf("%lf",x);}
      |                        ~~~~~^~~~~~~~~

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;

#define REP(i,a,b) for(i=a;i<b;i++)
#define rep(i,n) REP(i,0,n)

#define mygc(c) (c)=getchar_unlocked()
#define mypc(c) putchar_unlocked(c)

#define ll long long
#define ull unsigned ll

void reader(int *x){int k,m=0;*x=0;for(;;){mygc(k);if(k=='-'){m=1;break;}if('0'<=k&&k<='9'){*x=k-'0';break;}}for(;;){mygc(k);if(k<'0'||k>'9')break;*x=(*x)*10+k-'0';}if(m)(*x)=-(*x);}
void reader(unsigned *x){int k;*x=0;for(;;){mygc(k);if('0'<=k&&k<='9'){*x=k-'0';break;}}for(;;){mygc(k);if(k<'0'||k>'9')break;*x=(*x)*10+k-'0';}}
void reader(ll *x){int k,m=0;*x=0;for(;;){mygc(k);if(k=='-'){m=1;break;}if('0'<=k&&k<='9'){*x=k-'0';break;}}for(;;){mygc(k);if(k<'0'||k>'9')break;*x=(*x)*10+k-'0';}if(m)(*x)=-(*x);}
void reader(double *x){scanf("%lf",x);}
int reader(char c[]){int i,s=0;for(;;){mygc(i);if(i!=' '&&i!='\n'&&i!='\r'&&i!='\t'&&i!=EOF) break;}c[s++]=i;for(;;){mygc(i);if(i==' '||i=='\n'||i=='\r'||i=='\t'||i==EOF) break;c[s++]=i;}c[s]='\0';return s;}
template <class T, class S> void reader(T *x, S *y){reader(x);reader(y);}
template <class T, class S, class U> void reader(T *x, S *y, U *z){reader(x);reader(y);reader(z);}
template <class T, class S, class U, class V> void reader(T *x, S *y, U *z, V *w){reader(x);reader(y);reader(z);reader(w);}

void writer(int x, char c){int s=0,m=0;char f[10];if(x<0)m=1,x=-x;while(x)f[s++]=x%10,x/=10;if(!s)f[s++]=0;if(m)mypc('-');while(s--)mypc(f[s]+'0');mypc(c);}
void writer(unsigned x, char c){int s=0;char f[10];while(x)f[s++]=x%10,x/=10;if(!s)f[s++]=0;while(s--)mypc(f[s]+'0');mypc(c);}
void writer(ll x, char c){int s=0,m=0;char f[20];if(x<0)m=1,x=-x;while(x)f[s++]=x%10,x/=10;if(!s)f[s++]=0;if(m)mypc('-');while(s--)mypc(f[s]+'0');mypc(c);}
void writer(double x, char c){printf("%.15f",x);mypc(c);}
void writer(const char c[]){int i;for(i=0;c[i]!='\0';i++)mypc(c[i]);}
void writer(const char x[], char c){int i;for(i=0;x[i]!='\0';i++)mypc(x[i]);mypc(c);}
template<class T> void writerLn(T x){writer(x,'\n');}
template<class T, class S> void writerLn(T x, S y){writer(x,' ');writer(y,'\n');}
template<class T, class S, class U> void writerLn(T x, S y, U z){writer(x,' ');writer(y,' ');writer(z,'\n');}
template<class T> void writerArr(T x[], int n){int i;if(!n){mypc('\n');return;}rep(i,n-1)writer(x[i],' ');writer(x[n-1],'\n');}

void reader_u(int &x, char c){int k;x=0;for(;;){mygc(k);if(k==c)break;x=x*10+k-'0';}}
void writer_u(int x, char c){int s=0;char f[10];while(x)f[s++]=x%10,x/=10;if(!s)f[s++]=0;while(s--)mypc(f[s]+'0');mypc(c);}


char memarr[17000000]; void *mem = memarr;
#define INF 1000000000

template<class T> void sort(int N, T a[], void *mem = NULL){sort(a,a+N);}
template<class T1, class T2> void sort(int N, T1 a[], T2 b[], void *mem){int i;pair<T1,T2> *r=(pair<T1, T2>*)mem;rep(i,N)r[i].first=a[i],r[i].second=b[i];sort(r,r+N);rep(i,N)a[i]=r[i].first,b[i]=r[i].second;}
template<class T1, class T2, class T3> void sort(int N, T1 a[], T2 b[], T3 c[], void *mem){int i;pair<T1,pair<T2,T3> > *r=(pair<T1,pair<T2,T3> >*)mem;rep(i,N)r[i].first=a[i],r[i].second.first=b[i],r[i].second.second=c[i];sort(r,r+N);rep(i,N)a[i]=r[i].first,b[i]=r[i].second.first,c[i]=r[i].second.second;}
template<class T1, class T2, class T3, class T4> void sort(int N, T1 a[], T2 b[], T3 c[], T4 d[], void *mem){int i;pair<pair<T1,T2>,pair<T3,T4> > *r=(pair<pair<T1,T2>,pair<T3,T4> >*)mem;rep(i,N)r[i].first.first=a[i],r[i].first.second=b[i],r[i].second.first=c[i],r[i].second.second=d[i];sort(r,r+N);rep(i,N)a[i]=r[i].first.first,b[i]=r[i].first.second,c[i]=r[i].second.first,d[i]=r[i].second.second;}

struct unionFind{
  int *d, N, M;
  inline void malloc(int n){
    d = (int*)std::malloc(n*sizeof(int));
    M = n;
  }
  inline void *malloc(int n, void *mem){
    d = (int*)mem;
    M = n;
    return d + n;
  }
  inline void init(int n){
    int i;
    N = n;
    rep(i,n) d[i] = -1;
  }
  inline void init(void){
    init(M);
  }
  inline int get(int a){
    int t = a, k;
    while(d[t]>=0) t=d[t];
    while(d[a]>=0) k=d[a], d[a]=t, a=k;
    return a;
  }
  inline int connect(int a, int b){
    if(d[a]>=0) a=get(a);
    if(d[b]>=0) b=get(b);
    if(a==b) return 0;
    if(d[a] < d[b]) d[a] += d[b], d[b] = a;
    else            d[b] += d[a], d[a] = b;
    return 1;
  }
  inline int operator()(int a){
    return get(a);
  }
  inline int operator()(int a, int b){
    return connect(a,b);
  }
  inline int& operator[](int a){
    return d[a];
  }
  inline int sizeList(int res[]){
    int i, sz=0;
    rep(i,N) if(d[i]<0) res[sz++] = -d[i];
    return sz;
  }
};

template<class S, class T> inline void chmin(S &a, T b){if(a>b)a=b;}
template<class S, class T> inline void chmax(S &a, T b){if(a<b)a=b;}

#define SSE2_M128I_ADD32(c,a,b) asm volatile("movdqu %1, %%xmm0;" "movdqu %2, %%xmm1;" "paddd %%xmm1, %%xmm0;" "movdqu %%xmm0, %0;" :"=m"(c) :"m"(a), "m"(b) :);
#define SSE4_M128U_MIN32(c,a,b) asm volatile("movdqu %1, %%xmm0;" "movdqu %2, %%xmm1;" "pminud %%xmm1, %%xmm0;" "movdqu %%xmm0, %0;" :"=m"(c) :"m"(a), "m"(b) :);

// d = min(a,b+c)
#define SSE4_M128U_MINADD32(d,a,b,c) asm volatile("movdqu %2, %%xmm0;" "movdqu %3, %%xmm1;" "paddd %%xmm1, %%xmm0;" "movdqu %1, %%xmm1;" "pminud %%xmm1, %%xmm0;" "movdqu %%xmm0, %0;" :"=m"(d) :"m"(a), "m"(b), "m"(c) :);

#define AVX2_M256I_ADD32(c,a,b) asm volatile("vmovdqu %1, %%ymm0;" "vmovdqu %2, %%ymm1;" "vpaddd %%ymm1, %%ymm0, %%ymm0;" "vmovdqu %%ymm0, %0;" :"=m"(c) :"m"(a), "m"(b) :);
#define AVX2_M256I_MIN32(c,a,b) asm volatile("vmovdqu %1, %%ymm0;" "vmovdqu %2, %%ymm1;" "vpminud %%ymm1, %%ymm0, %%ymm0;" "vmovdqu %%ymm0, %0;" :"=m"(c) :"m"(a), "m"(b) :);

// d = min(a,b+c)
#define AVX2_M256U_MINADD32(d,a,b,c) asm volatile("vmovdqu %2, %%ymm0;" "vmovdqu %3, %%ymm1;" "vpaddd %%ymm1, %%ymm0, %%ymm0;" "vmovdqu %1, %%ymm1;" "vpminud %%ymm1, %%ymm0, %%ymm0;" "vmovdqu %%ymm0, %0;" :"=m"(d) :"m"(a), "m"(b), "m"(c) :);

void* readUndirectedGraph(int N, int M, int **es, int ***edge, void *mem, int origin=0){int i,*A;A=((int*)(((int**)mem)+N)+N+2*M);rep(i,2*M)reader(A+i),A[i]-=origin;*es=(int*)mem;*edge=(int**)((*es)+N);memset(*es,0,sizeof(int)*N);rep(i,2*M)(*es)[A[i]]++;(*edge)[0]=(int*)((*edge)+N);REP(i,1,N)(*edge)[i]=(*edge)[i-1]+(*es)[i-1];memset(*es,0,sizeof(int)*N);for(i=0;i<2*M;i+=2)(*edge)[A[i]][(*es)[A[i]]++]=A[i+1],(*edge)[A[i+1]][(*es)[A[i+1]]++]=A[i];return A;}
int connectedComponentSizeList(int N, int *es, int **edge, int res[], void *mem){
  int i, j, k, sz=0;
  int *st, sts, *vis;

  st = (int*)mem;
  vis = st + N;
  rep(i,N) vis[i] = 0;

  rep(i,N) if(vis[i]==0){
    res[sz] = 0;
    vis[i] = 1;
    st[0] = i; sts=1;
    while(sts){
      j = st[--sts];
      res[sz]++;
      rep(k,es[j]) if(vis[edge[j][k]]==0) vis[edge[j][k]]=1, st[sts++]=edge[j][k];
    }
    sz++;
  }

  return sz;
}

int N, M;
int dp[100001];
int cnt[100001], gain[100001];
int *es, **edge;

inline void update(unsigned g, unsigned c, unsigned st){
  int i;
  if(g < 8){
    for(i=st;i>=g;i--) chmin(dp[i], dp[i-g]+c);
  } else {
    unsigned cc[4] = {c,c,c,c};
    for(i=st;i-3>=g;i-=4){
      SSE4_M128U_MINADD32(dp[i-3], dp[i-3], dp[i-g-3], cc[0]);
    }
    for(;i>=g;i--) chmin(dp[i], dp[i-g]+c);
  }
}

int qs, qcost[100001], qgain[100001];

int main(){
  int i, j, k, st;

  reader_u(N,' ');
  reader_u(M,'\n');

  mem = readUndirectedGraph(N, M, &es, &edge, mem, 1);
  k = connectedComponentSizeList(N, es, edge, cnt, mem);
  rep(i,k) gain[cnt[i]]++;

  st = 0;
  REP(i,1,N+1) dp[i] = INF;
  REP(k,1,N+1) while(gain[k]){
    i = (gain[k]+1) / 2;
    qcost[qs] = i;
    qgain[qs++] = k*i;
    gain[k] -= i;
  }

  sort(qs, qgain, qcost, mem);
  rep(i,qs) update(qgain[i], qcost[i], st+=qgain[i]);

  REP(i,1,N+1){
    if(dp[i]==INF) mypc('-'), mypc('1'), mypc('\n');
    else writer_u(dp[i]-1,'\n');
  }

  return 0;
}
0