zoj 2974 Just Pour the Water (矩阵快速幂,简单)

题目

对于案例的解释请见下图:

zoj 2974 Just Pour the Water (矩阵快速幂,简单)

这道要变动提取一下矩阵,之后就简单了

具体解释可看代码:

#include <string.h>
#include <stdio.h>
#include<algorithm>
using namespace std; int num;
struct matrix
{
double a[][];//把每次倒水的比率提取出来放在这里面,例如i倒给j几分之几,以便进行计算
}origin,answ;//answ保存提取出来比率计算后的答案 matrix multiply(matrix x,matrix y)//x与y的矩阵乘法
{
matrix temp;
for(int i=;i<=num;i++)
{
for(int j=;j<=num;j++)
{
double ans=;
for(int k=;k<=num;k++)
{
ans+=x.a[i][k]*y.a[k][j];
}
temp.a[i][j]=ans;
}
} return temp;
} matrix calc(int n)//矩阵快速幂——把提取出来的比率进行n次重复的‘倒’,即n次幂
{
while(n)//以下运用二分
{
if(n%==)
answ=multiply(origin,answ);
origin=multiply(origin,origin);
n/=;
}
return answ;
} int main()
{
int t,n,k,m;
double a[];
scanf("%d",&t);
for(int id=;id<t;id++)
{
memset(answ.a,,sizeof(answ.a));
memset(origin.a,,sizeof(origin.a)); scanf("%d",&n);
for(int i=;i<=n;i++)
scanf("%lf",&a[i]); num=n;
for(int i=;i<=n;i++)
answ.a[i][i]=; for(int i=;i<=n;i++)
{
scanf("%d",&k);//i要平均的倒给之后的k个人
if(k==)
origin.a[i][i]=1.0;//不向任何地方倒,相当于都倒给自己,即1/1
for(int j=;j<=k;j++)
{
int p;
scanf("%d",&p);
origin.a[p][i]=1.0/k;//i倒给p的比率
}
} scanf("%d",&m);
calc(m);//对数据进行处理 for(int i=;i<=n;i++)//逐个输出答案
{
double ans=;//初始化
for(int j=;j<=n;j++)//一列一列处理
ans+=answ.a[i][j]*a[j];//处理后的比率*原始数据
if(i==)
printf("%.2lf",ans);
else
printf(" %.2lf",ans);
}
puts("");
}
return ;
}
上一篇:vue阻止事件冒泡和默认事件


下一篇:STL,ATL,WTL的联系与区别