#pragma GCC optimize("Ofast") #pragma GCC optimize("unroll-loops") #pragma GCC optimize("inline") #include using namespace std; template struct cLtraits_identity{ using type = T; } ; template using cLtraits_try_make_signed = typename conditional< is_integral::value, make_signed, cLtraits_identity >::type; template struct cLtraits_common_type{ using tS = typename cLtraits_try_make_signed::type; using tT = typename cLtraits_try_make_signed::type; using type = typename common_type::type; } ; void*wmem; char memarr[96000000]; template inline auto min_L(S a, T b) -> typename cLtraits_common_type::type{ return (typename cLtraits_common_type::type) a <= (typename cLtraits_common_type::type) b ? a : b; } template inline auto max_L(S a, T b) -> typename cLtraits_common_type::type{ return (typename cLtraits_common_type::type) a >= (typename cLtraits_common_type::type) b ? a : b; } template inline void walloc1d(T **arr, int x, void **mem = &wmem){ static int skip[16] = {0, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1}; (*mem) = (void*)( ((char*)(*mem)) + skip[((unsigned long long)(*mem)) & 15] ); (*arr)=(T*)(*mem); (*mem)=((*arr)+x); } template inline void walloc1d(T **arr, int x1, int x2, void **mem = &wmem){ walloc1d(arr, x2-x1, mem); (*arr) -= x1; } inline int my_getchar_unlocked(){ static char buf[1048576]; static int s = 1048576; static int e = 1048576; if(s == e && e == 1048576){ e = fread_unlocked(buf, 1, 1048576, stdin); s = 0; } if(s == e){ return EOF; } return buf[s++]; } inline void rd(int &x){ int k; int m=0; x=0; for(;;){ k = my_getchar_unlocked(); if(k=='-'){ m=1; break; } if('0'<=k&&k<='9'){ x=k-'0'; break; } } for(;;){ k = my_getchar_unlocked(); if(k<'0'||k>'9'){ break; } x=x*10+k-'0'; } if(m){ x=-x; } } struct MY_WRITER{ char buf[1048576]; int s; int e; MY_WRITER(){ s = 0; e = 1048576; } ~MY_WRITER(){ if(s){ fwrite_unlocked(buf, 1, s, stdout); } } } ; MY_WRITER MY_WRITER_VAR; void my_putchar_unlocked(int a){ if(MY_WRITER_VAR.s == MY_WRITER_VAR.e){ fwrite_unlocked(MY_WRITER_VAR.buf, 1, MY_WRITER_VAR.s, stdout); MY_WRITER_VAR.s = 0; } MY_WRITER_VAR.buf[MY_WRITER_VAR.s++] = a; } inline void wt_L(char a){ my_putchar_unlocked(a); } inline void wt_L(int x){ int s=0; int 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){ my_putchar_unlocked('-'); } while(s--){ my_putchar_unlocked(f[s]+'0'); } } inline void wt_L(unsigned x){ int s=0; char f[10]; while(x){ f[s++]=x%10; x/=10; } if(!s){ f[s++]=0; } while(s--){ my_putchar_unlocked(f[s]+'0'); } } inline void wt_L(long long x){ int s=0; int 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){ my_putchar_unlocked('-'); } while(s--){ my_putchar_unlocked(f[s]+'0'); } } inline void wt_L(unsigned long long x){ int s=0; char f[21]; while(x){ f[s++]=x%10; x/=10; } if(!s){ f[s++]=0; } while(s--){ my_putchar_unlocked(f[s]+'0'); } } int WRITER_DOUBLE_DIGIT = 15; inline int writerDigit_double(){ return WRITER_DOUBLE_DIGIT; } inline void writerDigit_double(int d){ WRITER_DOUBLE_DIGIT = d; } inline void wt_L(double x){ const int d = WRITER_DOUBLE_DIGIT; int k; int r; double v; if(x!=x || (x==x+1 && x==2*x)){ my_putchar_unlocked('E'); my_putchar_unlocked('r'); my_putchar_unlocked('r'); return; } if(x < 0){ my_putchar_unlocked('-'); x = -x; } x += 0.5 * pow(0.1, d); r = 0; v = 1; while(x >= 10*v){ v *= 10; r++; } while(r >= 0){ r--; k = floor(x / v); if(k >= 10){ k = 9; } if(k <= -1){ k = 0; } x -= k * v; v *= 0.1; my_putchar_unlocked(k + '0'); } if(d > 0){ my_putchar_unlocked('.'); v = 1; for(r=(0);r<(d);r++){ v *= 0.1; k = floor(x / v); if(k >= 10){ k = 9; } if(k <= -1){ k = 0; } x -= k * v; my_putchar_unlocked(k + '0'); } } } inline void wt_L(const char c[]){ int i=0; for(i=0;c[i]!='\0';i++){ my_putchar_unlocked(c[i]); } } inline void wt_L(string &x){ int i=0; for(i=0;x[i]!='\0';i++){ my_putchar_unlocked(x[i]); } } template inline S chmin(S &a, T b){ if(a>b){ a=b; } return a; } template struct LHeap{ int*hp; int*place; int size; T*val; void malloc(int N){ hp = (int*)std::malloc(N*sizeof(int)); place=(int*)std::malloc(N*sizeof(int)); val=(T*)std::malloc(N*sizeof(T)); } void malloc(int N, int ini){ hp = (int*)std::malloc(N*sizeof(int)); place=(int*)std::malloc(N*sizeof(int)); val=(T*)std::malloc(N*sizeof(T)); if(ini){ init(N); } } void walloc(int N, void **mem=&wmem){ walloc1d(&hp, N, mem); walloc1d(&place, N, mem); walloc1d(&val, N, mem); } void walloc(int N, int ini, void **mem=&wmem){ walloc1d(&hp, N, mem); walloc1d(&place, N, mem); walloc1d(&val, N, mem); if(ini){ init(N); } } void free(){ std::free(hp); std::free(place); std::free(val); } void init(int N){ int i; size=0; for(i=(0);i<(N);i++){ place[i]=-1; } } void up(int n){ int m; while(n){ m=(n-1)/2; if(val[hp[m]]<=val[hp[n]]){ break; } swap(hp[m],hp[n]); swap(place[hp[m]],place[hp[n]]); n=m; } } void down(int n){ int m; for(;;){ m=2*n+1; if(m>=size){ break; } if(m+1val[hp[m+1]]){ m++; } if(val[hp[m]]>=val[hp[n]]){ break; } swap(hp[m],hp[n]); swap(place[hp[m]],place[hp[n]]); n=m; } } void change(int n, T v){ T f = val[n]; val[n] = v; if(place[n]==-1){ place[n] = size; hp[size++] = n; up(place[n]); } else{ if(f < v){ down(place[n]); } else if(f > v){ up(place[n]); } } } int pop(void){ int res = hp[0]; place[res] = -1; size--; if(size){ hp[0]=hp[size]; place[hp[0]]=0; down(0); } return res; } } ; template struct minCostFlow{ int node; int*es; int*emem; int**edge; int**rev; FT**flow; FT f_eps; CT**cost; CT*potential; CT c_eps; LHeap hp; char*reached; FT*cur_flow; CT*cur_cost; int*back_edge; void malloc(int N){ int i; es = (int*)std::malloc(N*sizeof(int)); emem = (int*)std::malloc(N*sizeof(int)); edge = (int**)std::malloc(N*sizeof(int*)); rev = (int**)std::malloc(N*sizeof(int*)); flow = (FT**)std::malloc(N*sizeof(FT*)); cost = (CT**)std::malloc(N*sizeof(CT*)); for(i=(0);i<(N);i++){ emem[i] = 0; edge[i] = rev[i] = NULL; flow[i] = NULL; cost[i] = NULL; } hp.malloc(N); reached = (char*)std::malloc(N*sizeof(char)); cur_flow = (FT*)std::malloc(N*sizeof(FT)); cur_cost = (CT*)std::malloc(N*sizeof(CT)); potential = (CT*)std::malloc(N*sizeof(CT)); back_edge = (int*)std::malloc(N*sizeof(int)); node = N; for(i=(0);i<(N);i++){ es[i] = 0; } f_eps = (FT)1e-9; c_eps = (CT)1e-9; } void init(int N){ int i; node = N; for(i=(0);i<(N);i++){ es[i] = 0; } f_eps = (FT)1e-9; c_eps = (CT)1e-9; } void memoryExpand(int i, int sz){ if(sz <= emem[i]){ return; } sz =max_L(max_L(sz, 3), 2*emem[i]); emem[i] = sz; edge[i] = (int*)realloc(edge[i], sz*sizeof(int)); rev[i] = (int*)realloc(rev[i], sz*sizeof(int)); flow[i] = (FT*)realloc(flow[i], sz*sizeof(FT)); cost[i] = (CT*)realloc(cost[i], sz*sizeof(CT)); } void addEdge(int n1, int n2, FT f, CT c){ int s1 = es[n1]++; int s2 = es[n2]++; if(s1 >= emem[n1]){ memoryExpand(n1, es[n1]); } if(s2 >= emem[n2]){ memoryExpand(n2, es[n2]); } edge[n1][s1] = n2; edge[n2][s2] = n1; rev[n1][s1] = s2; rev[n2][s2] = s1; flow[n1][s1] = f; flow[n2][s2] = 0; cost[n1][s1] = c; cost[n2][s2] = -c; } template void solve(int st, int ed, FTS &fres, CTS &cres, FT flim = -1, CT clim = 0){ int i; int j; int k; int l; FT f; CT nc; fres = 0; cres = 0; for(i=(0);i<(node);i++){ potential[i] = 0; } for(;;){ if(flim >= -f_eps && flim <= f_eps){ break; } hp.init(node); for(i=(0);i<(node);i++){ reached[i] = 0; } reached[st] = 1; cur_cost[st] = 0; l = 0; hp.change(st, cur_cost[st]); while(hp.size){ i = hp.pop(); for(j=(0);j<(es[i]);j++){ if(flow[i][j] <= f_eps){ continue; } k = edge[i][j]; nc = cur_cost[i] + cost[i][j] + potential[i] - potential[k]; if(reached[k]==0 || cur_cost[k] > nc+c_eps){ reached[k] = 1; cur_cost[k] = nc; cur_flow[k] = flow[i][j]; if(i != st){ chmin(cur_flow[k], cur_flow[i]); } back_edge[k] = rev[i][j]; hp.change(k, cur_cost[k]); } } } if(reached[ed]==0){ break; } if(flim==-2 && cur_cost[ed] + potential[ed] >= clim){ break; } f = cur_flow[ed]; if(flim >= -f_eps){ chmin(f, flim); flim -= f; } if(f <= f_eps){ break; } for(i=(0);i<(node);i++){ if(reached[i]){ potential[i] += cur_cost[i]; } } fres += f; cres += f * potential[ed]; i = ed; while(i != st){ j = back_edge[i]; k = edge[i][j]; flow[i][j] += f; flow[k][rev[i][j]] -= f; i = k; } } } } ; int N; int K; int A[100000]; int M; int X[100000]; int Y[100000]; int main(){ wmem = memarr; int i; int j; int k; int node; int st; int ed; minCostFlow f; int flow; long long cost; rd(N); rd(K); for(i=(0);i<(N);i++){ int Q5VJL1cS; rd(A[i]); rd(k); for(Q5VJL1cS=(0);Q5VJL1cS<(k);Q5VJL1cS++){ rd(X[M]);X[M] += (-1); Y[M] = i; M++; } } node = N; st = 0; ed = N-1; f.malloc(node); f.init(node); for(i=(1);i<(N);i++){ f.addEdge(i-1,i,1073709056,1000000000); } for(i=(0);i<(M);i++){ if(A[Y[i]] > A[X[i]]){ f.addEdge(X[i],Y[i],1,-A[Y[i]]+A[X[i]]+(long long)1000000000*(Y[i]-X[i])); } } f.solve(st,ed,flow,cost,K); wt_L((long long)(N-1)*1000000000*K-cost); wt_L('\n'); return 0; } // cLay version 20210904-1 // --- original code --- // int N, K, A[1d5], M, X[], Y[]; // { // int i, j, k; // int node, st, ed; // minCostFlow f; // int flow; ll cost; // rd(N,K); // rep(i,N){ // rd(A[i],k); // rep(k){ // rd(X[M]--); // Y[M] = i; // M++; // } // } // // node = N; // st = 0; // ed = N-1; // f.malloc(node); // f.init(node); // // rep(i,1,N) f.addEdge(i-1,i,int_inf,1d9); // rep(i,M) if(A[Y[i]] > A[X[i]]) f.addEdge(X[i],Y[i],1,-A[Y[i]]+A[X[i]]+(ll)1d9*(Y[i]-X[i])); // // f.solve(st,ed,flow,cost,K); // Wt((ll)(N-1)*1d9*K-cost); // }