#pragma GCC target("avx") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #include // #include // #include // #include // using namespace __gnu_pbds; // #include // namespace multiprecisioninteger = boost::multiprecision; // using cint=multiprecisioninteger::cpp_int; using namespace std; using ll=long long; using datas=pair; using ddatas=pair; using tdata=pair; using vec=vector; using mat=vector; using pvec=vector; using pmat=vector; // using llset=tree,rb_tree_tag,tree_order_statistics_node_update>; #define For(i,a,b) for(i=a;i<(ll)b;++i) #define bFor(i,b,a) for(i=b,--i;i>=(ll)a;--i) #define rep(i,N) For(i,0,N) #define rep1(i,N) For(i,1,N) #define brep(i,N) bFor(i,N,0) #define brep1(i,N) bFor(i,N,1) #define all(v) (v).begin(),(v).end() #define allr(v) (v).rbegin(),(v).rend() #define vsort(v) sort(all(v)) #define vrsort(v) sort(allr(v)) #define uniq(v) vsort(v),(v).erase(unique(all(v)),(v).end()) #define endl "\n" #define popcount __builtin_popcountll #define eb emplace_back #define print(x) cout< ostream& operator<<(ostream& os,const pair& p){return os<<"("< ostream& operator<<(ostream& os,const vector& v){ os<<"{";bool f=false; for(auto& x:v){if(f)os<<",";os< ostream& operator<<(ostream& os,const set& v){ os<<"{";bool f=false; for(auto& x:v){if(f)os<<",";os< ostream& operator<<(ostream& os,const multiset& v){ os<<"{";bool f=false; for(auto& x:v){if(f)os<<",";os< ostream& operator<<(ostream& os,const map& v){ os<<"{";bool f=false; for(auto& x:v){if(f)os<<",";os< inline bool chmax(T& a,const T b){bool x=a inline bool chmin(T& a,const T b){bool x=a>b;if(x)a=b;return x;} #ifdef DEBUG void debugg(){cout<void debugg(const T& x,const Args&... args){cout<<" "< int main(){ startupcpp(); int N,K; cin>>N>>K; atcoder::mcf_graph fl(N); vector a(N); vector> edge; vec dp(N,0); for(int i=0;i>a[i]>>M; while(M--){ int x; cin>>x; if(a[--x]>=a[i])continue; edge.emplace_back(x,i,1,a[x]-a[i]); chmin(dp[i],dp[x]+a[x]-a[i]); // fl.add_edge(x,i,1,a[x]-a[i]); } if(i){ chmin(dp[i],dp[i-1]); edge.emplace_back(i-1,i,K,0); } // if(i)fl.add_edge(i-1,i,K,0); } debug(dp); for(auto[from,to,cap,cost]:edge){ debug(from,to,cap,dp[from]-dp[to]+cost); fl.add_edge(from,to,cap,dp[from]-dp[to]+cost); } auto ans=fl.flow(0,N-1,K); debug(ans); print(-ans.second-dp[N-1]*K); }