POJ 1487 Single-Player Games

最新推荐文章于 2019-03-30 20:40:59 发布
最新推荐文章于 2019-03-30 20:40:59 发布
阅读量187 收藏
点赞数
原文链接：http://www.cnblogs.com/lin375691011/p/3996793.html
版权

高斯消元第二题，这次的多元一次方程组的系数是浮点数，高斯消元的模版就有了些改动，但是主要思路还是那样。


题目大意：

给出一个数n，给出前n个字母的描述，问这个字母获得得知的期望是多少。


解题思路：

根据给出的描述列方程。主要难点在字符串处理上。


下面是代码：

#include <set>
#include <map>
#include <queue>
#include <math.h>
#include <vector>
#include <string>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>

#define eps 1e-6
#define pi acos(-1.0)
#define inf 107374182
#define inf64 1152921504606846976
#define lc l,m,tr<<1
#define rc m + 1,r,tr<<1|1
#define iabs(x)  ((x) > 0 ? (x) : -(x))
#define clear1(A, X, SIZE) memset(A, X, sizeof(A[0]) * (SIZE))
#define clearall(A, X) memset(A, X, sizeof(A))
#define memcopy1(A , X, SIZE) memcpy(A , X ,sizeof(X[0])*(SIZE))
#define memcopyall(A, X) memcpy(A , X ,sizeof(X))
#define max( x, y )  ( ((x) > (y)) ? (x) : (y) )
#define min( x, y )  ( ((x) < (y)) ? (x) : (y) )

using namespace std;

struct node
{
    double num[28];
    //int cnt[28];
    node()
    {
        clearall(num,0);
        //clearall(cnt,0);
    }
    void clen()
    {
        clearall(num,0);
        //clearall(cnt,0);
    }
};

struct node matrix[30];
int n,m,len,p;
bool free_x[305];

double X[305];

inline double zero(double x)
{
    if(fabs(x)>eps)return x;
    return 0;
}

void Debug(void)
{
    puts("");
    int i, j;
    for (i = 0; i < m; i++)
    {
        for (j = 0; j < n + 1; j++)
        {
            cout << matrix[i].num[j] << " ";
        }
        cout << endl;
    }
    cout << endl;
}

int Guass()
{
    int i,j,k,col;
    clearall(X,0);
    clearall(free_x,1);//把解集清空，所有变量都标为自由变量

    //Debug();

    for (k = 0,col = 0; k < m && col < n; ++k, ++col) //枚举行列
    {
        int max_r = k;//找到该col列元素绝对值最大的那行与第k行交换.(为了在除法时减小误差)
        for (i = k + 1; i < m; ++i)
        {
            if (iabs(matrix[i].num[col]) > iabs(matrix[max_r].num[col])) max_r = i;
        }
        if (max_r != k) //交换
        {
            for (i = k; i < n + 1; ++i) swap(matrix[k].num[i],matrix[max_r].num[i]);
        }
        if (iabs(matrix[k].num[col])<eps ) //如果对应该列都为0，枚举该行的下一列
        {
            k--;
            continue;
        }
        for (i = k + 1; i < m; ++i) //将k后边的col进行初等变换成行阶梯矩阵
        {
            if (iabs(matrix[i].num[col])>eps)
            {
                double d=matrix[i].num[col]/matrix[k].num[col];
                for (j = col; j < n + 1; ++j)
                {
                    matrix[i].num[j] =matrix[i].num[j] - d*matrix[k].num[j];
                }
            }
        }
    }

    //Debug();
    // 1. 无解的情况: 化简的增广阵中存在(0, 0, ..., a)这样的行(a != 0). 即R(A) != R(A')无解
    for (i = k; i < m; ++i)
    {
        if (iabs(matrix[i].num[col]) >eps) return -1;
    }
    // 2. 无穷解的情况: 在n * (n + 1)的增广阵中出现(0, 0, ..., 0)这样的行，即说明没有形成严格的上三角阵.
    // 且出现的行数即为自由变元的个数.   即R(A) = R(A') < n
    //printf("%d %d\n",k,n);
    if (k < n)
    {
        //注释处为求多解的自由变量
        // 首先，自由变元有n - k个，即不确定的变元至少有n - k个.
        int num = 0,freeidx;
        for (i = k - 1; i >= 0; --i)
        {
            num = 0;// 用于判断该行中的不确定的变元的个数，如果超过1个，则无法求解，它们仍然为不确定的变元.
            double tmp = matrix[i].num[n];
            // 第i行一定不会是(0, 0, ..., 0)的情况，因为这样的行是在第k行到第m行.
            // 同样，第i行一定不会是(0, 0, ..., a), a != 0的情况，这样的无解的.
            for (j = 0; j < n; ++j)
            {
                if (iabs(matrix[i].num[j]) > eps && free_x[j])
                {
                    num++;
                    freeidx = j;
                }
            }
            if (num > 1) continue; // 无法求解出确定的变元.
            // 说明就只有一个不确定的变元free_index，那么可以求解出该变元，且该变元是确定的.
            tmp = matrix[i].num[n];
            for (j = 0; j < n; ++j)
            {
                if (iabs(matrix[i].num[j])>eps && j != freeidx) tmp -= matrix[i].num[j]*X[j];
            }
            X[freeidx] = tmp/matrix[i].num[freeidx];
            free_x[freeidx] = 0;
        }
        return n - k;
    }
    // 3. 唯一解的情况: 在n * (n + 1)的增广阵中形成严格的上三角阵.
    // 计算出Xn-1, Xn-2 ... X0.
    for (i = k - 1; i >= 0; --i)
    {
        double tmp = matrix[i].num[n];
        for (j = i + 1; j < n; ++j)
        {
            tmp =tmp- matrix[i].num[j]*X[j];
        }
        X[i] = tmp/matrix[i].num[i];
    }
    return 0;
}
char s[250];
void dfs(node &a)
{
    node temp;
    int num,cnt=0;
    while(s[p]!=')')
    {
        if(s[p]>='a'&&s[p]<='z')
        {
            a.num[s[p]-'a']+=1;
            p++;
            cnt++;
        }
        else if(s[p]==' ')p++;
        else if(s[p]=='-')
        {
            p++;
            num=0;
            while(s[p]>='0'&&s[p]<='9')
            {
                num*=10;
                num+=s[p]-'0';
                p++;
            }
            a.num[n]-=num;
            cnt++;
        }
        else if(s[p]>='0'&&s[p]<='9')
        {
            num=0;
            while(s[p]>='0'&&s[p]<='9')
            {
                num*=10;
                num+=s[p]-'0';
                p++;
            }
            a.num[n]+=num;
            cnt++;
        }
        else if(s[p]=='(')
        {
            p++;
            temp.clen();
            dfs(temp);
            cnt++;
            for(int i=0; i<=n; i++)
            {
                a.num[i]+=temp.num[i];
            }
        }
    }
    if(cnt)
    {
        for(int i=0; i<27; i++)
        {
            a.num[i]/=cnt;
        }
    }
    p++;
    /*printf("%d\n",cnt);
    for(int i=0; i<n; i++)
    {
        for(int j=0; j<27; j++)
        {
            if(j!=0)printf(" ");
            printf("%.2lf",a.num[j]);
        }
        puts("");
    }*/
    return ;
}
int main()
{
   // freopen("data.out","w",stdout);
    int case1=1;
    while(scanf("%d",&n),n)
    {
        m=n;
        getchar();
        clearall(matrix,0);
        for(int i=0; i<n; i++)
        {
            gets(s);
            len=strlen(s);
            p=0;
            while(s[p]!='(')
            {
                p++;
            }
            p++;
            dfs(matrix[i]);
            for(int j=0; j<n; j++)
            {
                matrix[i].num[j]=-matrix[i].num[j];
            }
            matrix[i].num[i]++;
        }

        int sta=Guass();
        /*for(int i=0; i<n; i++)
        {
            for(int j=0; j<=n; j++)
            {
                if(j!=0)printf(" ");
                printf("%.2lf",matrix[i].num[j]);
            }
            puts("");
        }
        puts("");*/
        printf("Game %d\n",case1++);

        for(int i=0; i<n; i++)
        {
            printf("Expected score for %c ",i+'a');
            if(sta==-1)
            {
                puts("undefined");
            }
            else
            {

                    if(sta==0) printf("= %.3lf\n",zero(X[i]));
                    else if(free_x[i])puts("undefined");
                    else printf("= %.3lf\n",zero(X[i]));
            }
        }
        puts("");
    }
    return 0;
}


转载于:https://www.cnblogs.com/lin375691011/p/3996793.html

dhn37379
关注
poj1995 Raising Modulo Numbers （整数快速幂）
.
02-24 2万+
Raising Modulo Numbers Time Limit:1000MS Memory Limit:30000K Total Submissions:7280 Accepted:4397 Description People are different. Some secretly read magazines
poj 1995 快速幂二进制取模算法
张小V
08-24 1539
本题大意：求a1的b1次方加a2的b2次方一直加到an的bn次方，用他们幂的和对一个数x取余，把结果输出！ 参考出处：http://www.cnblogs.com/yan-boy/archive/2012/11/29/2795294.html（感谢原作者！） 矩阵的快速幂是用来高效地计算矩阵的高次方的。将朴素的o（n）的时间复杂度，降到log（n）。 这里先对原理（主要运
【POJ 1487】Single-Player Games（建方程高消）
ACMer'
08-22 1552
Single-Player Games Time Limit: 1000MSMemory Limit: 10000KTotal Submissions: 1387Accepted: 377DescriptionPlaying games is the most fun if other people take part. But other players are not always availa
poj 2065 SETI（高斯消元）
街角的花街糖
05-19 647
SETI Time Limit: 1000MS   Memory Limit: 30000K Total Submissions: 1469   Accepted: 899 Description For some years, quite a lot of work has been put into listening to
【POJ】2065 SETI
weixin_30753873的博客
09-04 79
题意：直接拿样例，37 abc。 a～z表示1～26，*表示0。 x0*1^0+x1*1^1+x2*1^2=1(mod 37) x0*2^0+x1*2^1+x2*2^2=2(mod 37) x0*3^0+x1*3^1+x2*3^2=3(mod 37) 高斯消元，除法x等于乘以x对p的逆元。 1 #include<cstdio> 2 #include<cstr...
[POJ2065]SETI（高斯消元）
zyf2000
01-12 833
题目描述传送门 题意： 良心博主才会给出题意… 直接说方程是什么吧… 首先a…z=1..26，*=0 读入p（模数且为质数），s（下标从0开始），s长度为n 那么求方程组 10x1+11x2...1n−1xn≡s1(modp)1^0x_1+1^1x_2...1^{n-1}x_n\equiv s_1\pmod p 20x1+21x2...2n−1xn≡s2(modp)2^0x_1+2^
poj 2065 SETI
codeStone
08-24 2726
SETI Time Limit: 1000MS   Memory Limit: 30000KB   64bit IO Format: %I64d & %I64u [Submit]   [Go Back]
POJ前面的题目算法思路【转】
梁山伯的专栏
01-08 5479
1000 A+B Problem 送分题 49％ 2005-5-7 1001 Exponentiation 高精度 85％ 2005-5-7 1002 487-3279 n/a 90％ 2005-5-7 1003 Hangover 送分题 62％ 2005-5-
poj 1027（开心消消乐 模拟 + 搜索）
看起来喜欢做枯燥事情的人或许只是更清楚自己想要什么罢了
09-13 1175
Total Submissions: 6288   Accepted: 2317 Description The game named "Same" is a single person game played on a 10 \Theta 15 board. Each square contains a ball colored red (R), green (G), o
高斯消元法
热门推荐
蜗牛
04-23 3万+
矩阵，逆矩阵，秩，矩阵方程，参考代码。
SETI（POJ-2065）
Alex_McAvoy的博客
03-30 442
Problem Description For some years, quite a lot of work has been put into listening to electromagnetic radio signals received from space, in order to understand what civilizations in distant galaxie...
poj1487
weixin_30730151的博客
11-13 119
题目大意： 给一棵递归树，看链接图片，从根节点开始对于每个节点往它的子节点移动，直到叶子节点停止。每个节点选哪一个孩子节点继续往下走是随机的(等概率)。然后叶子节点都会标记一个数值，记为走到该节点的得分。 输入条件：先输入整数n(n=0时结束),接下来有n行(n <= 26)，每一行会为前n个小写字母(每个字母作为一个变量)的描述。如a = (1 b)表示 f(a) = 1/2 * (1 ...
poj_1487 Single-Player Games（高斯消元+自由元）
凳子
12-23 1149
Single-Player Games Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 1479   Accepted: 407 Description Playing games is the most fun if other people take part. B
poj2065 SETI（飘忽的英文题面）
Coco_T的博客
10-08 411
译文： 描述 多年来，为了解存在于遥远星系中的其他文明可能要告诉我们的信息，我们投入了大量的人力物力倾听来自太空的电磁无线电信号。其中一种令UTS的科学家们特别感兴趣的信号源便是星云Stupidicus。最近，我们发现，如果把每条消息都被看做是一个整数序列，a0，a1，……； 如果找到了一个正确的p值，那么函数 的函数值总是在0~26（0<=f(k)<=26）之间。其中n是传递的消息的长
uva 664 single player game
chengouxuan的专栏
09-27 457
Single-Player Games  Playing games is the most fun if other people take part. But other players are not always available if you
【POJ1487】——Single-Player Games（高斯消元）
绝尘
08-22 567
【POJ1487】——Single-Player Games（高斯消元）
POJ2065 SETI 高斯消元
AC_Gibson的专栏
10-16 975
题目链接：http://poj.org/problem?id=2065 题目大意：给你一个素数P(P 定义: f (k) = ∑0 求a0、a1.....an-1。题目保证肯定有唯一解。 分析：高斯消元。根据上面的公式显然可以列出有n个未知数的n个方程式： a0*1^0 + a1*1^1+a2*1^2+........+an-1*1^(n-
Single-Player Games hdoj
生活总是要微笑着面对~
05-24 700
很久没有写题解了，因为我发现之前写的题没有写题解的现在再让我读题就忘的差不多了，比如置换群，当初看了好久，现在忘的差不多了，也许题解靠后写的好处是可以巩固啊！最近没有多少时间刷题了，快要考试了，最后两周还要数据库课程设计，还要学html和javascrpt，而我对着是完全没有接触过的，老师说以前都是几个人一组做一个大作业的，今年特殊，一个人做一个，所以，我是考前复习加上数据库设计，稍带还要备考六级
POJ 1010题解 - 巧妙设计思路与陷阱分析
资源摘要信息:"POJ 1010 邮票问题" 该压缩文件包含了关于POJ（北京大学在线评测系统）上的典型题目“邮票问题”（题目编号1010）的详尽解题信息。解题内容不仅包括源代码，而且还涵盖了在解题过程中需要注意的问题...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值