图论测试2 t2 幻想乡的符卡 二分答案+最大流

考试的时候花了好久来调这道题，结果忘了1的特判就gg了。但自己二分也确实有问题，建边也不是特别优。但基本的想法是没啥问题的。
就是二分一个等级，然后对符合条件的符卡建一个最大权闭合图，判断是否合法，再一次缩小答案的区间。
每次也不用想太多，直接重新拆了暴力建边就可以了23333之前以为这样会T掉，然而23333只要不是完全暴力（每次再重新判断两个的和是不是质数）（这个也只会t掉10%的点而已）都能过。
关于建图。
这是一个二分图，奇数放一边，偶数放一边。每一边的任意两个数的和都是偶数，肯定相加不是质数。每次二分等级level 后，选出符合这个等级的点来建图。（加入这个图的边一定已经满足等级<=level了）s向选出的每个偶数连一条流量为它的火力的边，每个偶数u向等级<=level且它们时间和为质数的奇数连一条inf的边，表示这两个必须有取舍，每个奇数再向t连一条流量为火力的边。
其实就是最大权闭合子图。
看起来似乎特别完美。
但：
1的特判！1的特判！1的特判！
无论这里面有多少个T等于1的点，我们都只能选其中的一个。所以每次找出符合level限制的火力最大的一个1（有点贪心的思想吧）加入图中。

#include<cstdio>#include<cstring>#include<ctime>#include<algorithm>#define ms(x,y) memset(x,y,sizeof(x))using namespace std;const int N = 400 + 10;const int M = 1e6 + 10;const int INF = 0x73f3f3f * 4;int n,k;int p[N],T[N],l[N];int mmax=-1,mmin=0x73f3f3f;int zzz=-1,lev;struct node{    int pre,v;    int f;}edge[M];int num=1;int head[N],cur[N];void addedge(int from,int to,int f){    num++;    edge[num].pre=head[from];    edge[num].v=to;    edge[num].f=f;    head[from]=num;    num++;    edge[num].pre=head[to];    edge[num].v=from;    edge[num].f=0;    head[to]=num;}int cnt=0,primes[M];bool isnot[M];void Prime(){    memset(isnot,0,sizeof(isnot));    isnot[1]=true;    for(int i=2;i<=M - 10;i++){        if(!isnot[i]){            cnt++;            primes[cnt]=i;        }        for(int j=1;j<=cnt;j++){            int u=i*primes[j];            if(u>M-10) break;            isnot[u]=true;            if(i%primes[j]==0) break;        }    }}int number[N][2];int s,t;int dis[N],state[M];bool vis[N];bool bfs(){    int h=0,tail=1;    dis[s]=0,vis[s]=true;    state[1]=s;    do{        h++;        int u=state[h];        for(int i=head[u];i;i=edge[i].pre){            int v=edge[i].v;            if(!vis[v]&&edge[i].f){                dis[v]=dis[u]+1;                vis[v]=true;                tail++;                state[tail]=v;            }        }    }while(h<tail);    if(vis[t]==true) return true;    return false;}inline int Min(int a,int b){    return a<b?a:b;}inline int Max(int a,int b){    return a>b?a:b;}int dfs(int u,int delta){    if(u==t||delta==0) return delta;    int ans=0;    for(int i=head[u];i&&delta;i=edge[i].pre){        int v=edge[i].v;        if(edge[i].f&&dis[v]==dis[u]+1){            int dd=dfs(v,Min(delta,edge[i].f));            edge[i].f-=dd;            edge[i^1].f+=dd;            ans+=dd;            delta-=dd;        }    }    if(!ans) dis[u]=-1;    return ans;}#define ms(x,y) memset(x,y,sizeof(x))void zero(){    ms(dis,0);ms(vis,0);ms(state,0);}int maxflow(){    int ans=0;    while(1){        zero();        if(!bfs()) break;        ans+=dfs(s,INF);    }    return ans;}bool check1(int x){    for(int i=1;i<=cnt;i++){        if(x%primes[i]==0) return false;        if(primes[i]*primes[i]>x) break;    }    return true;}int numb[N][N];void noname(){    ms(numb,0);    for(int i=1;i<=number[0][0];i++){        for(int j=1;j<=number[0][1];j++){            int u=number[i][0],v=number[j][1];            if(T[u]+T[v]<M-10){                 if(isnot[T[u]+T[v]]==false){                    numb[u][v]=1;                }            }            else if(check1(T[u]+T[v])) numb[u][v]=1;        }    }}int zx=0;bool check(int mid){    num=1,ms(head,0),zx=0,zzz=-1;    for(int i=1;i<=number[0][0];i++){        for(int j=1;j<=number[0][1];j++){            int u=number[i][0],v=number[j][1];            if(l[u]>mid||l[v]>mid||T[v]==1) continue;            if(numb[u][v]==1) addedge(u,v,INF);        }    }    for(int i=1;i<=number[0][0];i++){        if(l[number[i][0]]<=mid){          addedge(s,number[i][0],p[number[i][0]]);          zx+=p[number[i][0]];        }    }    for(int i=1;i<=number[0][1];i++){        if(T[number[i][1]]==1){            if(l[number[i][1]]>mid) continue;            if(zzz<p[number[i][1]]){                zzz=p[number[i][1]];                lev=number[i][1];            }            continue;        }        if(l[number[i][1]]<=mid){            addedge(number[i][1],t,p[number[i][1]]);            zx+=p[number[i][1]];        }    }    if(zzz!=-1){        for(int i=1;i<=number[0][0];i++){            int u=number[i][0];            if(numb[u][lev]==1) addedge(u,lev,INF);        }        zx+=zzz;        addedge(lev,t,zzz);    }    int xx=maxflow();    xx=zx-xx;    if(xx<k) return false;    return true;}int main(){    freopen("card.in","r",stdin);    freopen("card.out","w",stdout);    scanf("%d%d",&n,&k);    s=0,t=n+1;    Prime();    for(int i=1;i<=n;i++){        scanf("%d%d%d",&p[i],&T[i],&l[i]);        number[0][(T[i]&1)]++;        number[number[0][T[i]&1]][T[i]&1]=i;        mmax=Max(mmax,l[i]);        mmin=Min(mmin,l[i]);    }    noname();    if(!check(mmax)){       printf("-1\n");return 0;}    int L=mmin,R=mmax;    while(L<R){        num=1,ms(head,0);        int mid=(L+R)>>1;        if(check(mid)) R=mid;        else L=mid+1;    }    printf("%d",R);    return 0;}