【ZOJ 2974】Just Pour the Water(矩阵快速幂)

传送门:http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=2974

 

题意 给出n个杯子与初始水量同时进行操作 将其中的水同时平均分入所指定的杯子 进行x次后 输出杯子剩余水量

刚拿到这个题,第一反应是递推找规律,但是因为每个杯子的初始水量是未知的,所以能找的只是每个杯子水量与其余杯子水量的关系。

但是看到了操作次数巨大,而且最多只有20个杯子,感觉可以用快速幂去做。

我们假设矩阵a[i][j]代表第i个杯子的水有a[i][j]来自第j个杯子,这样用init矩阵表示初始水量,unit矩阵表示每次操作后的变化,最后init*pow(unit,k)即为所求。

注意当某一个杯子不用向其他杯子操作时,要将相应位置填为1,因为这个原因wa了8次。

 

总结:这是道思维大于算法的题,算法只是简单的矩阵快速幂,但是思维与经验很重要,比如对于数据的理解。由于操作数巨大,所以可以考虑一发快速幂。

代码如下:

///BY: Torrance_ZHANG
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
struct JZ {
    double a[25][25];
} init, unit;
int n;
JZ multi(JZ x, JZ y) {
    JZ z;
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < n; j++) {
            z.a[i][j] = 0.0;
            for (int k = 0; k < n; k++) {
                z.a[i][j] += (x.a[i][k] * y.a[k][j]);
            }
        }
    }
    return z;
}
JZ Pow(JZ x, JZ y, int k) {
    while (k) {
        if (k % 2 != 0)
            y = multi(y, x);
        x = multi(x, x);
        k /= 2;
    }
    return y;
}
void init_JZ(int n) {
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < n; j++) {
            init.a[i][j] = 0.0;
            unit.a[i][j] = 0.0;
        }
    }
}
double aa[25];
int main() {
    int t;
    scanf("%d", &t);
    while (t--) {
        scanf("%d", &n);
        init_JZ(n);
        for (int i = 0; i < n; i++)
            scanf("%lf", &init.a[0][i]);
        int tmp = 0, xh;
        for (int i = 0; i < n; i++) {
            scanf("%d", &tmp);
            if (tmp == 0)                ///最容易被忽略的一点,若输入为0,则全部的水都来自本身
                unit.a[i][i] = 1;
            else {
                for (int j = 0; j < tmp; j++) {
                    scanf("%d", &xh);
                    unit.a[i][xh - 1] += 1.0 / tmp;
                }
            }
        }
        int time;
        scanf("%d", &time);
        JZ res = Pow(unit, init, time);
        for (int i = 0; i < n; i++) {
            printf("%.2f", res.a[0][i]);
            printf(i == n - 1 ? "\n" : " ");
        }
    }
}

 

转载于:https://www.cnblogs.com/Torrance/p/5410755.html

Python网络爬虫与推荐算法新闻推荐平台:网络爬虫:通过Python实现新浪新闻的爬取,可爬取新闻页面上的标题、文本、图片、视频链接(保留排版) 推荐算法:权重衰减+标签推荐+区域推荐+热点推荐.zip项目工程资源经过严格测试可直接运行成功且功能正常的情况才上传,可轻松复刻,拿到资料包后可轻松复现出一样的项目,本人系统开发经验充足(全领域),有任何使用问题欢迎随时与我联系,我会及时为您解惑,提供帮助。 【资源内容】:包含完整源码+工程文件+说明(如有)等。答辩评审平均分达到96分,放心下载使用!可轻松复现,设计报告也可借鉴此项目,该资源内项目代码都经过测试运行成功,功能ok的情况下才上传的。 【提供帮助】:有任何使用问题欢迎随时与我联系,我会及时解答解惑,提供帮助 【附带帮助】:若还需要相关开发工具、学习资料等,我会提供帮助,提供资料,鼓励学习进步 【项目价值】:可用在相关项目设计中,皆可应用在项目、毕业设计、课程设计、期末/期中/大作业、工程实训、大创等学科竞赛比赛、初期项目立项、学习/练手等方面,可借鉴此优质项目实现复刻,设计报告也可借鉴此项目,也可基于此项目来扩展开发出更多功能 下载后请首先打开README文件(如有),项目工程可直接复现复刻,如果基础还行,也可在此程序基础上进行修改,以实现其它功能。供开源学习/技术交流/学习参考,勿用于商业用途。质量优质,放心下载使用。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值