hdu 4406 GPA

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4406

 

题意：

       距离期末考试还有n天，小A同学现在有m门课要复习，他计划着每天的复习时间是k分钟，现在告诉你他的每门课现在的分数和每门课的学分，他每在一门课上投入一分钟，这门课的最终成绩就会多加1，同时因为课程教学的问题，他每天都只能复习固定的几门课，现在在保证所有课都不挂科（不低于60分）的基础上，请你计算它的最大绩点，假设第i门课的绩点为w[i]，分数为so[i]，那么他最后的绩点会是.

做法：

       emmmm，还是很明显是最大费用流的，就是怎么建图的问题了，每一天都是一个点，每一门课也是一个点，这样的话对应的日子向对应的课程连一条容量为inf，费用为0的边，表示我每天可以在这门课上投入无限大的时间，源点向每一天连接流量为k的边，表示一天只能复习k分钟，问题来了，如何优先让自己不挂科呢。

        我们将离及格还缺少的分数向汇点连边，并附上一个无限小的值（我是用最小费用的板子跑的，无非就是建负边），这样的话就会优先跑这个流量，最后我们看费用有没有超过我们的sum流量*无限小费用，如果没有超过的话就是还没跑完那些优先要跑的边，自然也就不会有其他的边了，所以我们就做完啦~

---

#include <bits/stdc++.h>
#define rep(i,a,b) for(int i=(int)a;i<=(int)b;i++)
using namespace std;
typedef long long ll;
const int maxn=105;
const int maxm=3500;
const int MAX= 100000;
const double inf=1e6;
const double eps=1e-6;
double dis[maxn];
int vis[maxn],pre[maxn];
int head[maxn],cnt;
int n,m,sp,tp,k,neadd;
int flows[maxn];
struct node{
    int to,cap,next;
    double cost;
}e[maxm];
int dcmp(double x){
    if(fabs(x)<eps) return 0;
    if(x>0) return 1;
    return -1;
}
void add(int from,int to,int cap,double cost){
    //printf("from = %d  to = %d  cap = %d  cost = %.3f\n",from,to,cap,cost);
    e[cnt].to=to; e[cnt].cap=cap;
    e[cnt].cost=cost; e[cnt].next=head[from];
    head[from]=cnt++;

    e[cnt].to=from; e[cnt].cap=0;
    e[cnt].cost=-cost; e[cnt].next=head[to];
    head[to]=cnt++;
}
bool spfa(int s,int t,int &flow,double &cost){
    queue<int> q;
    for(int i=0;i<=65;i++){
        dis[i]=inf;vis[i]=0;pre[i]=-1; flows[i]=inf;
    }
    dis[s]=0;  q.push(s);
    vis[s]=1; flows[s]=MAX;
    while(!q.empty()){
        int u=q.front();
        q.pop();
        vis[u]=0;
        for(int i=head[u];~i;i=e[i].next){
            int v=e[i].to;
            if(e[i].cap>0&&dcmp(dis[v]-(dis[u]+e[i].cost))>0){
                dis[v]=dis[u]+e[i].cost;
                pre[v]=i; flows[v]=min(flows[u],e[i].cap);
                if(!vis[v]){
                    vis[v]=1;
                    q.push(v);
                }
            }
        }
    }
    if(dcmp(dis[t]-inf)==0){
        return false;
    }
    int d=flows[t];
    for(int i=pre[t];~i;i=pre[e[i^1].to]){
        e[i].cap-=d;
        e[i^1].cap+=d;
        cost+=(double)e[i].cost*d;
    }
    flow+=d;
    return true;
}
double mcmf(int s,int t){
    int flow=0;
    double cost=0;
    while(spfa(s,t,flow,cost)){

    }
    return cost;
}
int w[maxn];
double soc[105],ans;
int bas[maxn];
double Cal(double qu,int x){
    return (double)qu*soc[x];
}

int main(){
    for(int i=60;i<=100;i++) soc[i]=(double)(4.0-3.0*(100-i)*(100-i)/1600.0);
    while(~scanf("%d%d%d",&n,&k,&m)){
        if(!n&&!m&&!k) break;
        memset(head,-1,sizeof(head));
        cnt=0; neadd=0;
        sp=0;tp=61; ans=0;
        double sum=0;
        for(int i=1;i<=n;i++) add(sp,i,k,0);
        rep(i,1,m) scanf("%d",&w[i]),sum+=w[i];
        rep(i,1,m){
            scanf("%d",&bas[i]),neadd+=max(0,60-bas[i]);
            if(bas[i]<60){
                add(i+40,tp,60-bas[i],-inf);
                bas[i]=60;
            }
            ans=ans+Cal(w[i],bas[i]);
            for(int j=bas[i]+1;j<=100;j++){
                //printf("j=%d Cal j = %.5f , Cal j-1 = %.5f\n",j,Cal(w[i],j),Cal(w[i],j-1));
                add(i+40,tp,1,-(Cal(w[i],j)-Cal(w[i],j-1)));
            }
        }

        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                int x;
                scanf("%d",&x);
                if(x) add(i,j+40,MAX,0);
            }
        }
        double ret=-mcmf(sp,tp);
        if(dcmp(ret-neadd*(inf)<0)) {
            printf("0.000000\n");
        }
        else{
            ret=ret-inf*neadd;
            printf("%.6f\n",(ans+ret)/sum);
        }
    }
    return 0;
}