【dp】2014(高斯消元 DP 数论)Multi-University Training Contest 1 题解

http://blog.sina.com.cn/s/blog_6bddecdc0102uy9g.html

A.Couple doubi(hdu 4861):

 

 

#include<stdio.h> int main(){ int k,p; while(~scanf("%d %d",&k,&p)){ int num=k/(p-1); if(num%2==0) printf("NO\n"); else printf("YES\n"); } return 0; }

 

 

 

 

 

D.Task(hdu 4864)

基本思想是贪心。

价值c=500*time+2*level，level最大影响100*2<500.所以近似求time最大，可以先对任务按照时间大到小排序（时间相等时按照等级由大到小排序）。从最大时间的任务开始，找出满足任务时间要求的所有机器，从中找出满足等级最低且满足等级要求的机器与之匹配。依次对任务寻找满足要求的机器(能匹配当前的任务，那么剩下的所有的任务都是能匹配到的，那么这个机器一定会被用掉，当然选价值最大的。最大匹配数无论如何都会加1，如果不选当前的任务，那么剩下的机器中可能没有能完成该任务的了。而其他任务在剩下的机器中还可能得到匹配)

 

 

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
#define TIME 1449
#define LEVEL 109
#define MAX 100009
struct Node {
    int time;
    int level;
} task[MAX],mach[MAX];
int map[TIME][LEVEL];
bool cmp(struct Node a, struct Node b){
    if(a.time==b.time) return a.level > b.level;
    return a.time > b.time;
}

int main() {
    int n,m;
    while(~scanf("%d %d",&n,&m)) {
        memset(map,0,sizeof(map));
        for(int i=0; i<n; i++) {
            scanf("%d %d",&mach[i].time,&mach[i].level);
            map[mach[i].time][mach[i].level]++;
        }
        for(int i=0; i<m; i++) {
            scanf("%d %d",&task[i].time,&task[i].level);
        }
        sort(task,task+m,cmp);
        __int64 ans=0,cnt=0;
        bool ok=false;
        for(int i=0; i<m; i++) {
            ok=false;
            for(int t=task[i].level; t<LEVEL; t++) {
                for(int j=task[i].time; j<TIME; j++) {
                    if(map[j][t]!=0) {
                        cnt++;
                        map[j][t]--;
                        ans+=(__int64)(500*task[i].time+2*task[i].level);
                        ok=true;
                        break;
                    }
                }
                if(ok==true) break;
            }
        }
        printf("%I64d %I64d\n",cnt,ans);
    }
    return 0;
}
/***
1 2
300 3
100 2
200 1
***/

 

E.Peter's Hobby(hdu 4865)

 

类型：dp

本题的题意是：给出n天的叶子的湿度，让你求出天气最有可能出现的序列。

叶子湿度序列：{b1,b2...bn},若天气序列为： {a1,a2....an}，则该序列的概率应为：

f[a1]*wh[a1][b1]* ww[a1][a2] *wh[a2][b2]* ww[a2][a3]* wh[a3][b3]*........*ww[an-1][an]*wh[an][bn];

最后选择概率最大的序列即可。

由于概率的值非常小，同过log把乘法变成加法，减少精度的损失。

 

#include<stdio.h>
#include<string.h>
#include <math.h>
#include<algorithm>
using namespace std;
#define eps 1e-7
const int N=100005;
double f[20050][5];
int pre[20050][5];
double d[3][3]= {
    {0.5 ,0.375,0.125},
    {0.25,0.125,0.625},
    {0.25,0.375,0.375}
};
double p[3][4]= {
    {0.6 ,0.2 ,0.15,0.05},
    {0.25,0.3,0.2,0.25},
    {0.05,0.10,0.35,0.50}
};
char ss[3][8]= {"Sunny","Cloudy","Rainy"};
int check(char *s) {
    if(strcmp(s,"Dry")==0)  return 0;
    if(strcmp(s,"Dryish")==0)  return 1;
    if(strcmp(s,"Damp")==0)  return 2;
    if(strcmp(s,"Soggy")==0)  return 3;
    return 0;
}

void output(int time,int k){
    if(time==0) return;
    output(time-1,pre[time][k]);
    printf("%s\n",ss[k]);
}
int main() {
    int T,n;
    char str[100];
    while(~scanf("%d",&T)) {
        for(int cas=1; cas<=T; cas++) {
            printf("Case #%d:\n",cas);
            scanf("%d",&n);
            f[1][0]=0.63;
            f[1][1]=0.17;
            f[1][2]=0.20;
            for(int i=1; i<=n; i++ ) {
                scanf("%s",str);
                int k=check(str);
                for(int j=0; j<3; j++) {
                    if(i==1) {
                        f[i][j]=log(f[i][j])+log(p[j][k]);
                    } else {
                        f[i][j]=-100000;
                        for(int t=0; t<3; t++) {
                            if(f[i][j]<f[i-1][t]+log(d[t][j])+log(p[j][k])) {
                                f[i][j]=f[i-1][t]+log(d[t][j])+log(p[j][k]);
                                pre[i][j]=t;
                            }
                        }
                    }
                }
                memset(str,0,sizeof(str));
            }
            int k=0;
            if(f[n][1]>f[n][k]) k=1;
            if(f[n][2]>f[n][k]) k=2;
            output(n,k);
        }
    }
    return 0;
}

 

 

I. Turn the pokers（hdu 4869）

 

 

http://www.cnblogs.com/sayeter/p/3864184.html

由于mod=10^9这么大，f[MAX]我竟然用了int,调了很久才知道是这个地方错了，用long long int(10^18)刚刚好

 

#include<stdio.h>
#define MAX 100009
#define mod 1000000009
int n,m,x;
long long int f[MAX];
void deal(){
    f[0]=1;
    for(int i=1;i<MAX;i++){
        f[i]=(f[i-1]*i)%mod;
    }
}
long long quick_mod(long long a,long long b){
    long long ans=1;
    while(b){
        if(b&1){
            ans=(ans*a)%mod;
            b--;
        }
        b=b/2;
        a=a*a%mod;
    }
    return ans;
}
int main(){
    deal();
    while(~scanf("%d %d",&n,&m)){
        int s=0,e=0,S,E;
        for(int i=0;i<n;i++){
            scanf("%d",&x);
            if(x<=s) S=s-x;
            else if(x<=e) S=((s%2)==(x%2)?0:1);
            else S=x-e;
            if(x+e<=m) E=e+x;
            else if(x+s<=m) E=((x+s)%2==(m%2)?m:m-1);
            else E=2*m-(x+s);
            s=S;
            e=E;
        }
        long long ans=0;
        for(int i=S;i<=E;i+=2){
            ans=ans+(f[m]*quick_mod((f[m-i]*f[i])%mod,mod-2))%mod;
            ans=ans%mod;
        }
        printf("%lld\n",ans);
    }
return 0;
}
/***
2 43
35 43
s=8 e=8 S=8 E=8
145008513
5 46
46 10 45 43 10
s=22 e=46 S=22 E=46
556929028
2 32
10 32
s=22 e=22 S=22 E=22
64512240
5 47
47 33 27 46 45
s=10 e=44 S=10 E=44
418169790
13 50
50 1 4 5 7 8 2 3 4 5 1 2 2
s=6 e=50 S=6 E=50
948123245
****

Rating (hdu 4870)

 

//http://www.cnblogs.com/luyingfeng/p/3864693.html
/**
题意：有两个号，初始分数都是0，每次选一个分数较小的打比赛，如果分数一样任选一个，有p的概率涨50分，
最高为1000分，有1-p的概率跌100分，最低为0分。问有一个号涨到1000需要打比赛的次数的期望。

令(x, y)表示高分为x，低分为y的状态（x >= y），E(x, y)表示从(x, y)到达(1000, ?)的比赛场数期望。
容易得到E(x, y) = P * E(x1, y1) + (1 - P) * E(x2, y2) + 1，其中，(x1, y1)表示rating上升后的状态，
（x2, y2）表示rating下降后的状态。每50分一个状态，共有210个状态。
移项后得E(x, y) -P * E(x1, y1) - (1 - P) * E(x2, y2) = 1，共有210个这样的方程组。高斯消元求解，
x[0]代表E(0,0)这个状态到目标状态的期望。 
注意精度。
*/
#include <stdio.h>
#include <string.h>
#include <iostream>
#include <math.h>

using namespace std ;

const double eps = 1e-9 ;

double p,mp[250][250];// mp数组中mp[i]表示第i个方程组 
int vis[25][25] ;//表示分数为E(i,j)的这个变量对应的编号 

double Gauss() {//消元完成后，只有对角线上和最后一列会有系数，其他的地方全为零。 
    for(int i = 0 ; i < 210 ; i++) {
        int k ;
        for(k = i ; k < 210 ; k++)
            if(fabs(mp[k][i]) > eps) break ;//用第i列不为零的行(>=i的行)，进行消元 
        for(int j = 0 ; j <= 210 ; j++)
            swap(mp[i][j],mp[k][j]) ;
        for(int j = 0 ; j < 210 ; j++) {
            if(i != j) {
                if(fabs(mp[j][i]) > eps) {//防止除数为零 
                    double s = mp[j][i] / mp[i][i] ;
                    for(int k = i ; k <= 210 ; k++)
                        mp[j][k] -= mp[i][k] * s ;
                }
            }
        }
    }
    return mp[0][210]/mp[0][0] ;//mp[0][210]为:a1*x1+a2*x2+....an*xn=b1;中的b1，那么mp[0][210]/mp[0][0]即为E[0][0]的值了 
}
/***
E(x,y)中(x>=y&&x>=0&&x<=20&&y>=0&&y<=19)那么这样的组合刚好有(21*20)/2这么多个变量，则需要建210个方程来解 
同时：状态中不用计入状态(1000,?)是因为除了(1000,950)外其他状态不会出现。而(1000,950)状态只与状态(950,950)有关
并且(950,950)不会从(1000,950)得到，210个状态中没有状态是与(1000,?)状态相关，所以不需要处理。(E(1000,?)==0)

*/
void pro() {//该函数用于构造非齐次线性方程组的增广矩阵(系数阵+常数阵) 
    for(int i = 0 ; i < 20 ; i++) {
        for(int j = 0 ; j < i ; j ++) {
            int u = vis[i][j] ;
            mp[u][u] = 1 ;
            mp[u][210] = 1 ;
            int v = vis[i][max(0,j-2)] ;
            mp[u][v] -= (1-p) ;//找到对应的变量E(u,v)在mp[u]这个方程中的位置 
            v = vis[i][j+1] ;
            mp[u][v] -= p ;
        }
        int u = vis[i][i] ;
        mp[u][u] = 1 ;
        mp[u][210] = 1 ;
        int v = vis[i][max(0,i-2)] ;
        mp[u][v] -= (1-p) ;
        v = vis[i+1][i] ;
        mp[u][v] -= p ;
    }
}

int main() {
    while(scanf("%lf",&p) != EOF) {
        int cnt = 0 ;
        memset(vis,-1,sizeof(vis)) ;
        for(int i = 0 ; i < 20  ; i++)
            for(int j = 0 ; j <= i ; j++)
                vis[i][j] = cnt ++ ;
        for(int i=0; i<20; i++) {
            for(int j=0; j<=i; j++) {
                printf("%d ",vis[i][j]);
            }
            puts("");
        }
        memset(mp,0,sizeof(mp)) ;
        pro() ;
        printf("%.6lf\n",Gauss()) ;
    }
    return 0 ;
}