题目要读很久才能理解它的意思和笑点(如果你也看过那个笑话的话),读懂之后就会发现是一个高斯消元法的题目,对于我来说难点不在高斯消元,而在于字符串处理。先来说说题意吧:
总共有n个人,n个人都会有一段话,先是princess说话,里面如果提到了a1,a2,a3...这几个不同的人的话,对应提到的次数是x1,x2,x3..的话,那么下一个对话是ai这个人说的概率是xi/(x1+x2+x3)....,然后下一个人的对话里也会提到别的人,然后也有一定的概率会有下一轮对话,现在要问的就是,给定了这些对话,问你期望的对话次数是多少。
我们可以设第i个人持续的对话的期望是xi,那么xi应该等于 xi=p1*x1+p2*x2+p3*x3+...+1 对应的pi即为在该人对话出现的次数所对应的频率。然后对每个人列出这样的方程就会构成一系列的方程组。然后高斯消元即可。
难度在于算在某个对话里出现了多少次,因为像如何你用KMP去做匹配 prince是会出现在princess那里的,题目说了人名之间是用五种间隔符隔开的,所以匹配的时候或者匹配前将对话里的每一个词分出来,然后和对应的字符串判相等即可。我写的不好跪了好几发呀- -0
#pragma warning(disable:4996) #include<iostream> #include<cstring> #include<cstdio> #include<algorithm> #include<cmath> #include<vector> #define ll long long #define maxn 220 #define eps 1e-7 using namespace std; char str[120][1050]; char name[120][20]; int n; double mat[120][120]; char word[1200]; int dcmp(double x) { return (x > eps) - (x < -eps); } int cnt(char*T, char *S) { int res = 0; int idx = 0; for (int i = 0; T[i]; i++){ if (!(T[i] <= 'z'&&T[i] >= 'a')) continue; int id = 0; while (T[i] <= 'z'&&T[i] >= 'a'){ word[id++] = T[i++]; } word[id] = ' '; if (strcmp(word, S) == 0) res++; i--; } return res; } bool gauss(double mat[120][120], int n) { for (int i = 0; i < n; i++){ int pivot = i; for (int j = i; j < n; j++){ if (abs(mat[j][i])>abs(mat[pivot][i])) pivot = j; } swap(mat[i], mat[pivot]); if (abs(mat[i][i]) < eps) return false; for (int j = i + 1; j <= n; j++) mat[i][j] /= mat[i][i]; for (int j = 0; j < n; j++){ if (i != j) for (int k = i + 1; k <= n; k++){ mat[j][k] -= mat[j][i] * mat[i][k]; } } } return true; } int main() { while (cin >> n) { getchar(); memset(name, 0, sizeof(name)); memset(str, 0, sizeof(str)); for (int i = 0; i < n; i++){ gets(str[i]); int j; for (j = 0; str[i][j]; j++){ if (str[i][j] == ':') break; name[i][j] = str[i][j]; } name[i][j + 1] = ' '; int k = 0; j++; for (; str[i][j]; j++,k++){ str[i][k] = str[i][j]; } str[i][k] = ' '; } memset(mat, 0, sizeof(mat)); for (int i = 0; i < n; i++){ double tot = 0; double tmp = 0; for (int j = 0; j < n; j++){ if (i == j) continue; tmp = cnt(str[i], name[j]); mat[i][j] = tmp; tot += tmp; } mat[i][i] = -tot; mat[i][n] = -tot; if (dcmp(tot) == 0){ mat[i][i] = mat[i][n] = -1; } } if (gauss(mat, n)){ printf("%.3lf ", mat[0][n]); } else puts("Infinity"); } return 0; }