哈尔滨理工大学第12届程序设计竞赛（牛客）题解

已于 2022-04-03 20:41:17 修改
阅读量368 收藏
点赞数
分类专栏： # 哈尔滨理工大学第12届程序设计竞赛（牛客）题解 ACM 文章标签： c++ 算法
于 2022-04-03 17:14:59 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/zbbxndm/article/details/123940143
版权

 涉及知识点： 逆元，组合数优化，欧拉降幂，分式化简，DFS

B题：

思路：

        显然概率分析有P=/;因为q次的概率其实是相同的（有放回），对其概率求逆元有:                                                qpow(C(n,k)，q,mod)*qpow(C(n+m,k),mod-2,mod); 

这里我们使用到了快速幂；其次我们要注意的是求组合数的时候，我发现一开始老是T了；

错误组合数代码：

long long FA(long long a) //定义阶乘函数FA
{
    long long b=1;//定义变量b
    for(int i=1;i<=a;i++)//计算阶乘
        b=(q_mul(b,i,mod))%mod;
    return b;//返回值得到b=a！
}
int F(int n,int k)
{
    return (FA(n)%mod)/(FA(k)*FA(n-k))%mod;
}

这里我使用该效率高的组合数代码：

const int maxn = 1e6 + 5;
 
int f(int n,int m){
    static int M=0,inv[maxn],mul[maxn],invMul[maxn];
    while(M<=n){
        if(M){
            inv[M]=M==1?1:(mod-mod/M)*inv[mod%M]%mod;
            mul[M]=mul[M-1]*M%mod;
            invMul[M]=invMul[M-1]*inv[M]%mod;
        }
        else mul[M]=1,invMul[M]=1;
        M++;
    }
    return mul[n]*invMul[m]%mod*invMul[n-m]%mod;
}

这里给出AC代码：

#include<queue>
#include<iostream>
#include<algorithm>
#include<string.h>
#include<cstdio>
#include <map>
#include<string>
#include<cstring>
#include<iomanip>
#include<set>
#include<map>
#include<numeric>
#include<math.h>
using namespace std;
#define int  long long
typedef unsigned long long ull;
typedef double db;
typedef long double ldb;
#define heap priority_queue
#define pll pair<long long, long long>
#define fir first
#define sec second
#define SPO(n) fixed <<setprecision(n)
#define FOR(i, l, r) for (long long i = l; i <= r; ++i)
#define ROF(i, r, l) for (long long i = r; i >= l; --i)
#define endl '\n'
const double pi=3.1415926535;
const double EPS = 1e-9;
const long long MAX = 1e9+5;
int  n,sum,m,q,k;
const int mod=1e9+7;


unsigned int q_mul(int a, int b, int mod)//快速乘法
{
    unsigned int ans = 0;
    while (b)
    {
        if (b & 1)
        {
            b--;
            ans = (ans + a) % mod;
        }
        b >>= 1;
        a = (a + a) % mod;
    }
    return ans;
}
int qpow(int a, int n, int mod)//快速幂
{
    int ans = 1;
    while (n)
    {
        if (n & 01)ans = ans * a % mod;
        a = a * a % mod;
        n >>= 1;
    }
    return ans;
}
//求组合数！！！效率高！！！！
const int maxn = 1e6 + 5;

int f(int n,int m){
    static int M=0,inv[maxn],mul[maxn],invMul[maxn];
    while(M<=n){
        if(M){
            inv[M]=M==1?1:(mod-mod/M)*inv[mod%M]%mod;
            mul[M]=mul[M-1]*M%mod;
            invMul[M]=invMul[M-1]*inv[M]%mod;
        }
        else mul[M]=1,invMul[M]=1;
        M++;
    }
    return mul[n]*invMul[m]%mod*invMul[n-m]%mod;
}

void Solve(void) {
    cin>>n>>m>>k>>q;
    if(n==0||k>n)
    {
        cout<<0<<endl;
        return ;
    }
    int ans1=f(n,k);
    int ans2=f(n+m,k);
    ans1= qpow(ans1,q,mod);
    ans2= qpow(ans2,q,mod);
    cout<<ans1*qpow(ans2,(mod-2),mod)%mod<<endl;
    return ;
}

signed main() {
    std::ios::sync_with_stdio(false);
    cin.tie(0);cout.tie(0);
    int t;
    cin>>t;
    while(t--)
    Solve();

}

 J题：

题解：

         这里我们发现 其实就是 欧拉降幂！！！当然我们只需要特判0^0=1就可以了（前提是要学过欧拉特判），这里我们给出C语言代码：

上AC代码：

#include <bits/stdc++.h>
#define ll long long
#define mod 10000000007
using namespace std;
char a[1000006];
ll x,z;
ll quickpow(ll x,ll y,ll z)
{
    ll ans=1;
    while(y)
    {
        if(y&1)
            ans=ans*x%z;
        x=x*x%z;
        y>>=1;
    }
    return ans;
}
ll phi(ll n)
{
    ll i,rea=n;
    for(i=2;i*i<=n;i++)
    {
        if(n%i==0)
        {
            rea=rea-rea/i;
            while(n%i==0)
                n/=i;
         }
    }
    if(n>1)
        rea=rea-rea/n;
    return rea;
}
int main()
{
	int t;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%lld %s %lld",&x,a,&z);
		ll len=strlen(a);
        ll p=phi(z);
        ll ans=0;
        for(ll i=0;i<len;i++)
            ans=(ans*10+a[i]-'0')%p;
        ans+=p;
        printf("%lld\n",quickpow(x,ans,z));
        
	}
    return 0;
}

C题：

 题解;

       因为没有障碍物，我们可以直接坐标差最小

#include<queue>
#include<iostream>
#include<algorithm>
#include<string.h>
#include<cstdio>
#include <map>
#include<string>
#include<cstring>
#include<iomanip>
#include<set>
#include<map>
#include<numeric>
#include<math.h>
using namespace std;
#define ll  long long
typedef unsigned long long ull;
typedef double db;
typedef long double ldb;
#define heap priority_queue
#define pll pair<long long, long long>
#define fir first
#define sec second
#define SPO(n) fixed <<setprecision(n)
#define FOR(i, l, r) for (long long i = l; i <= r; ++i)
#define ROF(i, r, l) for (long long i = r; i >= l; --i)
#define endl '\n'
const double pi=3.1415926535;
const double EPS = 1e-9;
const long long MAX = 1e5+5;
int  n,sum,m,q,k;

int dx[7]={1,-1,0,0,0,0};
int dy[7]={0,0,1,-1,0,0,};
int dz[7]={0,0,0,0,1,-1};
void Solve()
{
    int n,m,q,op,x,h,y,z;
    int st=0,mp[4][MAX],s;
    cin>>n>>m>>h>>q;
    while(q--)
    {
        cin>>op>>x>>y>>z;
        if(op==1)
        {
            mp[0][st]=x,mp[1][st]=y,mp[2][st]=z;
            st++;continue;
        }else
        {
            s=1e5;
            for(int i=0;i<st;++i)
            {
                int k=abs(mp[0][i]-x)+ abs(mp[1][i]-y)+abs(mp[2][i]-z);
                if(s>k)s=k;
            }
            cout<<s<<endl;
        }
    }
}
signed main() {
    std::ios::sync_with_stdio(false);
    cin.tie(0);cout.tie(0);
    int t;
    cin>>t;
    while(t--)
        Solve();

}

但是，我们为了熟悉bfs也可以这个样子

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll MAXN=5e5+5;
ll n,m,h,qw;
struct pp
{
    ll x,y,z,dis;
};
ll change(ll x,ll y,ll z)
{
    return x-1+(y-1)*n+(z-1)*n*m;
}
ll vis[MAXN];
ll dx[10]={1,0,0,-1,0,0};
ll dy[10]={0,1,0,0,-1,0};
ll dz[10]={0,0,1,0,0,-1};
ll ans[MAXN];
ll mnx=1e18,mxx=0,mny=1e18,mxy=0,mnz=1e18,mxz=0;
bool check(ll x,ll y,ll z)
{
    return x >= 1 && x <=n  && y >= 1 && y <= m && z >= 1 && z <= h;
}
queue<pp>q;
void bfs()
{
    while(!q.empty())
    {
        pp ro=q.front();
        q.pop();
        for(int i=0;i<6;++i)
        {
            ll tx=ro.x+dx[i];
            ll ty=ro.y+dy[i];
            ll tz=ro.z+dz[i];
            if(!check(tx,ty,tz))continue;
            if(vis[change(tx,ty,tz)]>ro.dis+1) {
                q.push({tx, ty, tz, ro.dis + 1});
                vis[change(tx,ty,tz)]=ro.dis+1;
            }
        }
    }
}
int main()
{
    ll t;
    cin>>t;
    while(t--)
    {
        cin>>n>>m>>h>>qw;
        for(int i=0;i<=n*m*h;++i)vis[i]=1e18;
        for(int i=0;i<qw;++i)
        {
            ll op,x,y,z;
            scanf("%lld %lld %lld %lld",&op,&x,&y,&z);
            if(op==1)
            {
                q.push({x,y,z,1});
                vis[change(x,y,z)]=1;
            }
            else {
                bfs();
                printf("%lld\n", vis[change(x, y, z)] - 1);
            }
        }
    }
}

 

补充知识点：

 1.欧拉降幂：

求a^b(mod p)的值 ，当b很大很大很大很大很大的时候，可以使用欧拉降幂

欧拉定理：若n，a为正整数，且n，a互质，则：{​{\color{Red} }a^{\phi (n)}\equiv 1(mod n)}

2.C++逆元

                       qpow(C(n,k)，q,mod)*qpow(C(n+m,k),mod-2,mod); 

即 分子*qpow(分母，mod-2,mod);      (1/4的逆元为250000002；

3.分式化简

/*分数化简，关键是求出分母和分子的最大公因数。
  这里采用“辗转相除法”求两个整数的最大公因数。
  辗转相除法， 又名欧几里德算法（Euclidean algorithm），是求最大公约数的一种方法。
  它的具体做法是：用较大数除以较小数，再用出现的余数（第一余数）去除除数，再用出现的余数（第二余数）去除第一余数，
  如此反复，直到最后余数是0为止。如果是求两个数的最大公约数，那么最后的除数就是这两个数的最大公约数。
*/
 
#include<iostream>
 
using namespace std;
 
int main() {
	int n, m;	
	cout << "请分别输入分母和分子：" << endl;
	cin >> m >> n;
	int m0 = m, n0 = n;    //m0和n0保存m和n的原始值
	while (m%n != 0) {
		int temp = m ;
		m = n;
		n = temp%n;
	}
	cout << "原分式的最简分式为： " << n0 / n << "/" << m0 / n << endl;		                                                   
}

underivativity_life
关注
专栏目录
题解 | #武汉工程大学第六届ACM程序设计竞赛（同步赛）#
2301_79125431的博客
04-29 1758
各位，如果来了，就来了，别慌了，看命吧， 各位，如果来了，就来了，别慌了，看命吧，慌改变不了什么，还是要打工吃饭的。#铜五铁六真的存在吗？(51062)# #铜五铁六真的存在吗？# 我选铜五必存！#牛客帮帮团来啦！有问必答(50227) #牛客帮帮团来啦！有问必答# 25届双非java后端求助。请问各位大佬，现在双非后端很难进大厂，所以杭州亚信科技 26号做的笔试，有牛友收到面试通知了吗本菜鸟双非零实习零竞
广西大学东信杯第四届程序设计竞赛（同步赛）部分题解
m0_60283288的博客
12-13 2692
广西大学东信杯第四届
分数的加减乘除及化简c++
01-18
分数的加减乘除及化简c++分数的加减乘除及化简c++分数的加减乘除及化简c++
分数化简
热门推荐
思
04-09 1万+
2019.04.09 将给定的真分数化简为最简分数形式，如：将6/8化为3/4 c++ /*分数化简，关键是求出分母和分子的最大公因数。 这里采用“辗转相除法”求两个整数的最大公因数。 辗转相除法， 又名欧几里德算法（Euclidean algorithm），是求最大公约数的一种方法。 它的具体做法是：用较大数除以较小数，再用出现的余数（第一余数）去除除数，再用出现的余数（...
真分式化简
weixin_45395062的博客
10-13 2292
分式化简。设计算法，将一个给定的真分数化简为最简分数形式。例如，将6/8化简为3/4。 #include<iostream> using namespace std; int main() { int m,n; int a,b; cout<<"输入分数的分子与分母："; cin>>n>>m; a=m; b=n; int r=m%n; while(r!=0) { a=n;
分数的化简
u013246297的专栏
02-12 1469
分数的化简 Time Limit: 1000ms Memory limit: 65536K 有疑问？点这里^_^ 题目描述 输入一个分数，对该分数化简后输出其最简分式的形式。 输入 输入一个分数。 输出 输出其真最简分式的形式。 示例输入 6/9 示例输出 2/3   分析：大神们就是大神！我这小人物只能看看水题，等会一
哈尔滨理工大学第12届程序设计竞赛--ADK题
weixin_45772744的博客
04-05 780
A割韭菜--线段树 题目链接: A-割韭菜_哈尔滨理工大学第12届程序设计竞赛（同步赛） (nowcoder.com) #include <bits/stdc++.h> using namespace std; typedef long long ll; const int N=1e5+10; int t,n,m; int a[N]; typedef struct Node{ int l,r; ll maxh,minh;//区间里面韭菜的最大和最小高度 ll maxv,minv
（牛客2021）第一届河北工业大学程序设计竞赛（部分题解）
m0_46335115的博客
05-03 1364
A.WELCOME! 直接输出 #include<bits/stdc++.h> using namespace std; int main(){ printf("Welcome to The First Programming Competition of Hebei University of technology!"); return 0; } B.POOLING 题意： 给出n * m大小的图，求出该图中每个大小为k * k的子图中的最大的数是多少。 思路： 因为n和m很
湖南省第十届"蓝狐网络杯"大学生计算机程序设计竞赛,2018湖南省第14届大学生计算机程序设计竞赛 C: 时间旅行...
weixin_42517607的博客
06-24 632
Description 假设 Bobo 位于时间轴(数轴)上t0点,他要使用时间机器回到区间(0, h]中. 当 Bobo 位于时间轴上t点,同时时间机器有c单位燃料时,他可以选择一个满足⌈xh⌉⋅h≤c⌈xh⌉⋅h≤c的非负整数x, 那么时间机器会在[0, x]中随机整数y,使 Bobo 回到(t − y)点,同时消耗y单位燃料. (其中⌈ ⋅ ⌉表示上取整...
校园招聘信息交流平台WEB，UNIAPP小程序（仿牛客、offershow、codetop）作为毕业设计简单实现功能
最新发布
05-14
【标题】中的“校园招聘信息交流平台WEB，UNIAPP小程序”指的是一个综合性的网络应用，旨在为大学生提供一个集中的地方获取实习、就业信息，并进行相关的交流互动。它仿照了像“牛客”、“offershow”、“codetop”...
C++（分数类，进行加减乘除的运算并化简）
FFFPAG的博客
06-18 5552
要求： （1）定义整数类和分数类。其中，包括构造函数、析构函数、显示函数等。 （2）输入 / 输出：对流提取和流插入运算符进行重载。 （3）计算功能：可进行分数的加、减、乘和除法运算。 （4）化简功能：将分数化简为最简分数。 （5）异常处理功能：分数中分母不能为零。 （6）菜单功能：每种功能的操作都是在菜单中进行相应选择。 程序： #include<iostream> #include<cmath> using namespace std; //整数 class Integer{
哈尔滨理工大学第12届程序设计竞赛（同步赛）
Cauchy的博客
04-03 1488
哈尔滨理工大学第12届程序设计竞赛（同步赛） 文章目录A 割韭菜B 抓球C 迷宫D gk的爬山之旅E gk的字符串F gk的树G gk的数字游戏H 进来DPI 又AK了J 大数乘法K 蹦蹦炸弹L NP-hard A 割韭菜 B 抓球 通过Note: 概率p=a/b 答案=a∗bmod−2%mod {答案=a*b^{mod-2}\%mod} 答案=a∗bmod−2%mod 题意：n白m黑，取k个黑，问q次全黑的“概率” 思路： 1. 如下公式，可知我们需要预处理阶乘，不然会TLE 2. 对除法转换成
哈尔滨理工大学第12届程序设计竞赛 B 抓球
追殇 | Blog
04-04 634
链接：https://ac.nowcoder.com/acm/contest/30825/B 题目描述 一个箱子里有n个黑球，m个白球，小王想要连续q次从箱子里随机的取出k个球（每次取出k个后立即放回），连续q次取球k个都为黑球的概率是多少，结果对1e9+7取模。 输入描述: 第一行给出一个t，代表t组测试数据。（1<=t<=100000） 每组测试数据给出n，m，k，q。（1<=n,m,k<=1e5，k<=n+m，1<=q<=1e12） 输出描述: 对于每组测试数
分数化小数
姚光超的专栏
10-12 1万+
《编程之美》有一个题是给定一个小数，将其转化成最简分数，思路比较简单，首先将小数转化成分数，然后对分数化简。如果将问题倒过来，给定一个分数（N/D），将其转化成对应的小数，这该如何做？ 我们先分析一下分数转化成小数的可能情况：1）小数是一个有限小数（0.abc…d）；2）小数是纯循环小数（0.）；3）小数是非纯循环小数（0.ab…c）。分数不可能产生无限不循环小数。第一种情况无需特殊考虑，只需要
哈尔滨理工大学第12届程序设计竞赛C题
m0_64045085的博客
04-04 573
迷宫
「五大常用算法」一文图解分治算法和思想
好好学java
12-16 296
点击上方好好学java，选择星标公众号重磅资讯、干货，第一时间送达 今日推荐：硬刚一周，3W字总结，一年的经验告诉你如何准备校招！ 个人原创100W+访问量博客：点击前往，查...
哈尔滨理工大学第12届程序设计竞赛（同步赛）错题笔记
热爱计算机语言！
04-03 3583
目录官网链接J 大数乘法 官网链接 J 大数乘法 题目描述 给定整数x,y，你需要输出 xy mod p 的值 输入描述: 第一行一个整数t（1 ≤ t ≤ 10）,表示接下来有 t 组测试用例，每个测试用例有三行整数x,y,p。 0 ≤ x ≤ 100000，0 ≤ y ≤ 10100000,100000 ≤ p ≤ 1000000007 输出描述: 输出t行，每行一个整数代表xy mod p 的值。 示例1 输入 5 3 4 998244353 2 10 998244353 0 100 99824435
ICPC 山东省省赛刷题 哈尔滨理工大学第12届程序设计竞赛（同步赛）牛客 EGIL题解
weixin_51273276的博客
05-18 259
E gk的字符串 题意：给你一个带有’?'的字符串，要求使用小写字母替代，相邻的字母不可以相同，要求字典序最小，输出这个字符串。 分析：字典序最小使用’a’替换’?‘，如果前后有’a’就使用’b’替换，如果前后有’a’也有’b’，就用’c’替换，这三种情况。 题解： #include<bits/stdc++.h> using namespace std; typedef long long ll; int main() { int t; cin >> t; while(t--)
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值