主要是参考别人的代码和想法做的
这题是网络流建模的题
奇偶性建模，主要是用二分图的形式，奇数二分图为一边，偶数为另外一边，主要是要注意1。 我将我的dinic模板抽象为下面这样，仅供参考。
主要是用前向星的方式存边 struct E{ int to,next,w; }edge[maxn*200]; int head[maxn]; void init() { memset(head,-1,sizeof(head)); } void addedge(int u,int v,int w) { edge[++tot].to=v; edge[tot].w=w; edge[tot].next=head[u]; head[u]=tot; edge[++tot].to=u; edge[tot].w=0; edge[tot].next=head[v]; head[v]=tot; } void build_graph() { init();//初始化 //这一部分主要是用来建图 } int level[maxn],s,t; bool bfs() { queue<int>q; memset(level,-1,sizeof(level)); level[s]=0; q.push(s); while(!q.empty()){ int u=q.front();q.pop(); for(int i=head[u];i!=-1;i=edge[i].next){ int v=edge[i].to; if(level[v]==-1&&edge[i].w>0){ level[v]=level[u]+1; q.push(v); } } } if(level[t]==-1) return false; else return true; } int dfs(int u,int flow) { if(u==t) return flow; int tmp, res = 0; for(int i=head[u];i!=-1;i=edge[i].next){ int v=edge[i].to; if(level[v]==level[u]+1&&edge[i].w>0){ tmp=dfs(v,min(edge[i].w,flow)); edge[i].w-=tmp; edge[i^1].w+=tmp; res += tmp; flow -= tmp; } } return res; } int dinic() { int flow=0; while (bfs()) { flow += dfs(s,inf); } return flow; }下面是这题的代码
#include <bits/stdc++.h> using namespace std; typedef long long LL; const LL INF=1e17L; const int inf=0x3f3f3f3f; const int maxn=1e2+5; const int mod=1e9+7; int n,k,cnt=0; struct node{ int p,c,l; }a[maxn]; bool isprime(int x) { if(x==2) return true; int t=sqrt(x); for(int i=2;i<=t;i++) if(x%i==0) return false; return true; } struct E{ int to,next,w; }edge[maxn*200]; int head[maxn],tot,s,t,total; void init() { memset(head,-1,sizeof(head)); tot=-1; } void addedge(int u,int v,int w) { edge[++tot].to=v; edge[tot].w=w; edge[tot].next=head[u]; head[u]=tot; edge[++tot].to=u; edge[tot].w=0; edge[tot].next=head[v]; head[v]=tot; } void build_graph(int mid) { s=0;t=n+1;total=0; init(); int p1=-1; for(int i=1;i<=n;i++){ if(a[i].c==1 && a[i].l<=mid){ if(p1==-1 || a[i].p>a[p1].p) p1=i; } } if(p1!=-1){ addedge(s,p1,a[p1].p); total+=a[p1].p; for(int i=1;i<=n;i++) if(a[i].c!=1 && isprime(1+a[i].c) && a[i].l<=mid) addedge(p1,i,inf); } for(int i=1;i<=n;i++){ if(a[i].c==1 || a[i].l>mid) continue; if(a[i].c&1) addedge(s,i,a[i].p); else addedge(i,t,a[i].p); total+=a[i].p; for(int j=i+1;j<=n;j++){ if(a[j].c==1 || a[j].l>mid || !isprime(a[i].c+a[j].c)) continue; if(a[i].c&1) addedge(i,j,inf); else addedge(j,i,inf); } } } int level[maxn]; bool bfs() { queue<int>q; memset(level,-1,sizeof(level)); level[s]=0; q.push(s); while(!q.empty()){ int u=q.front();q.pop(); for(int i=head[u];i!=-1;i=edge[i].next){ int v=edge[i].to; if(level[v]==-1&&edge[i].w>0){ level[v]=level[u]+1; q.push(v); } } } if(level[t]==-1) return false; else return true; } int dfs(int u,int flow) { if(u==t) return flow; int tmp, res = 0; for(int i=head[u];i!=-1;i=edge[i].next){ int v=edge[i].to; if(level[v]==level[u]+1&&edge[i].w>0){ tmp=dfs(v,min(edge[i].w,flow)); edge[i].w-=tmp; edge[i^1].w+=tmp; res += tmp; flow -= tmp; } } return res; } int dinic() { int flow=0; while (bfs()) { flow += dfs(s,inf); } return flow; } int main() { cin>>n>>k; for(int i=1;i<=n;i++) cin>>a[i].p>>a[i].c>>a[i].l; int l=1,r=100,ans=-1; while(l<=r){ int mid=(l+r)>>1; build_graph(mid); int cut=dinic(); if(total-cut>=k){ ans=mid; r=mid-1; } else l=mid+1; } cout<<ans<<endl; return 0; }