jzoj2018.7.6【2018提高组】模拟B组

1.二分图计数
10

#include <cstdio>
#include <cstring>

using namespace std;

const int q[18]={1,2,4,8,16,32,64,128,256,512,1024,2048,4096,8192,16384,32768,65536,131072};
const long long N=1000000007;
int n,m,g,d[17],e[17],a[17],h[17]; 
long long ans,b;
long long f[17];
long long s,ss;

void work(int tot){ 
    int bz=1;
    s=tot;
    f[tot+1]=1;
    for (int i=tot;i>=1;i--)
        f[i]=f[i+1]*(m-i+1);
    b=f[1]%N;
    for (int i=1;i<=d[0];i++) h[i]=d[i];
    for (int i=3;i<=d[0];i++){
        ss=0;
        for (int j=1;j<=d[0];j++)
            ss+=(s-h[j])*d[j],h[j]=(s-h[j])*d[j];
        if (bz) b-=(ss/i)*f[i]; else b+=(ss/i)*f[i];
        bz^=1;b%=N;s=ss;
    }
    if (d[0]==2) b-=f[2]*tot;
    if (d[0]==1) b--;
}

void dfs(int i,int k,long long c){
    if (i>n){
        if (k) work(k),ans=(ans+(b%N)*(c%N))%N;
        return;
    }
    dfs(i+1,k,c);
    for (int j=1;j<=d[0];j++)
        if (e[j]==a[i]) {g=j;break;}
    if (g) d[g]++; else d[++d[0]]++,e[d[0]]=a[i];
    dfs(i+1,k+1,c+q[i-1]);
    if (g) d[g]--,g=0; else e[d[0]]=e[0],d[d[0]--]=0;
}

int main(){
//  freopen("bipartite.in","r",stdin);
//  freopen("bipartite.out","w",stdout);
    scanf("%d%d",&n,&m);
    memset(e,0x3f,sizeof(e));
    for (int i=1;i<=n;i++)
        scanf("%d",&a[i]);
    dfs(1,0,0);
    printf("%lld",ans);
}

2.扑克游戏
100
=合并果子

=贪心，每次找最小的两个相加计入答案，删去这两个，把这两个的和加入

#include <cstdio>
#include <cstring>

using namespace std;

int n,a[10004],ans;

void down(int x,int n){
    while (x*2<n){
        int p=x*2;
        if (a[x]>a[p]||a[x]>a[p+1]){
            if (a[p+1]<a[p]&&p+1<=n) p++;
            int t=a[x];
            a[x]=a[p];a[p]=t;
        }
        x=p;
    }
}

int main(){
    scanf("%d",&n);
    memset(a,0x3f,sizeof(a));
    for (int i=1;i<=n;i++){
        scanf("%d",&a[i]);
    }
    for (int i=n;i>=1;i--) down(i,n);
    for (int i=n;i>=2;i--){
        int p;
        if (a[2]<a[3]) p=2; else p=3;       
        ans+=a[p]+a[1];
        a[p]+=a[1];
        a[1]=a[i];
        a[i]=a[0];
        down(p,i-1);
        down(1,i-1);
    }
    printf("%d",ans);
}

3.表白
100
二分+dp
二分答案，dp判断可行性
科学的题解：
那么我们二分t 之后可以得到每个人的两个参数Q1i-C1i*t 以及Q2i-C2i*t，采取最科学合理的方式安排人手即可。对于小数据我们可以采用O(N^3)的动态规划，即f[i][j][k] 表示前i 个人中挑了j 个人到第一队，挑了k个人到第二队，转移方程是显而易见的。
对于100%的大数据我们需要分析观察数据，我们发现，对N1+N2个人，我们只考虑如何安排这些人时，我们一定可以按照他们的(Q1i-C1i*t)-(Q2i-C2i*t)排序之后进行选择。同样的，我们对N个人先进行排序之后，用f[i][j] 表示前i 个人取了j 个进第一队，用g[i][j] 表示后i 个人取了j 个进第二队，最后用f[i][N1] +g[N-i][N2]来更新答案即可。另外注
意精度误差就可以了。
被精度误差卡了一会，发现读入时不要用小数读

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>

using namespace std;

const double e=0.000000001;
int q1[505],c1[505],q2[505],c2[505];
double l,r,mid;
int n,n1,n2,bz[505];
struct node{
    double a1,a2,a3;
}a[505];
double f[505][505],g[505][505];

bool comp(node a,node b){
    return a.a3>b.a3;
}

bool work(double k){
    for (int i=1;i<=n;i++){
        a[i].a1=(double)q1[i]-c1[i]*k;
        a[i].a2=(double)q2[i]-c2[i]*k;
        a[i].a3=a[i].a1-a[i].a2;
    }
    sort(a+1,a+1+n,comp);   
    memset(f,200,sizeof(f));
    memset(g,200,sizeof(g));
    g[0][0]=0;f[0][0]=0;
    for (int i=1;i<=n;i++){
        g[i][0]=0;
        for (int j=1;j<=min(n1,i);j++)
            g[i][j]=max(g[i-1][j],g[i-1][j-1]+a[i].a1);
    }
    for (int i=n;i>=1;i--)
    {
        f[n-i+1][0]=0;
        for (int j=1;j<=min(n-i+1,n2);j++)
            f[n-i+1][j]=max(f[n-i][j],f[n-i][j-1]+a[i].a2);
    }
    for (int i=n1;i<=n-n2;i++)
    if (g[i][n1]+f[n-i][n2]>=0) return true;
    return false;
}

int main(){
    scanf("%d%d%d",&n,&n1,&n2);
    for (int i=1;i<=n;i++){
        scanf("%d%d%d%d",&q1[i],&c1[i],&q2[i],&c2[i]);
    }
    l=0.01,r=2000;
    while (l+e<r){
        mid=(l+r)/2;
        if (work(mid)) l=mid;
                else r=mid-e;
    }
    printf("%.6f",l);
}