“蔚来杯“2022牛客暑期多校训练营1 GADI

dai _ tu

已于 2022-11-10 10:45:04 修改

阅读量251

点赞数 2

分类专栏：算法文章标签： c++ 算法图论

于 2022-07-19 15:01:51 首次发布

本文链接：https://blog.csdn.net/qq_57862276/article/details/125870669

版权

算法专栏收录该内容

14 篇文章 0 订阅

订阅专栏

这篇博客涵盖了三场编程竞赛的解题策略。第一场涉及处理字符串，检查是否全由9组成；第二场是覆盖区间问题，求解覆盖所需的最小长度；第三场涉及概率计算，用到了概率dp和除数取模的知识。博主分享了每道题目的解题思路和代码实现，并提供了快速幂等算法来简化除数取模操作。

摘要由CSDN通过智能技术生成

G.Villages: Landlines

官方题解：
请添加图片描述

代码如下：

#include<bits/stdc++.h>
using namespace std;
string s;
typedef long long ll;
bool flag=0;
int main()
{
	cin>>s;
	if(s.length()==1)
	{
	    cout<<s;
	    return 0;
	}
	ll len=s.length()-1;
	for(ll i=0;i<len;i++)
	{
		if(s[i]!='9')
		{
		    flag=1;
		    break;
		}	
	}
	if(flag)
	{
		for(ll i=1;i<=len;i++)
		  cout<<9;
	}
	else
	{
		cout<<s;
		return 0;
	}
	return 0;
	
}

A.Villages: Landlines

请添加图片描述
解：覆盖区间问题
代码如下：

#include<bits/stdc++.h>
using namespace std;
int n;
const int N=2e5+11;
typedef long long ll;
struct node
{
	ll l,r;
}e[N];
bool cmp(node x,node y)
{	
	return x.l<y.l;
}
ll ans=0;
int main()
{
	cin>>n;
ll x,y;
	for(int i=1;i<=n;i++)
	{
		cin>>x>>y;
	    e[i].l=x-y;
	    e[i].r=y+x;
	}
	sort(e+1,e+1+n,cmp);
	ll r=e[1].r;
	for(int i=2;i<=n;i++)
	{
		if(e[i].l>r)
		    ans+=e[i].l-r;
		r=max(e[i].r,r);
	}
	cout<<ans;
	return 0;
}

D.Mocha and Railgun

请添加图片描述
在这里插入图片描述
代码如下：

#include<bits/stdc++.h>
using namespace std;
#define endl '\n'
double x,y,r,d,t;
double ans,m,th,l;
void solve()
{
	double dis=sqrt(x*x+y*y);
	double y1=sqrt(r*r-(dis-d)*(dis-d)),y2=sqrt(r*r-(dis+d)*(dis+d));
	d=sqrt((y1-y2)*(y1-y2)+(2*d)*(2*d));
    double th=d/(2*r);
	ans=2*r*asin(th);
}
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);
	cin>>t;
	while(t--)
	{
		cin>>r>>x>>y>>d;
	    solve();
		printf("%.10lf",ans);
	}
	return 0;
 }

I.Chiitoitsu

官方题解：
在这里插入图片描述
请添加图片描述
解：

1.知识点：概率dp，除数取模。
其中"除数取模":
根据 费马小定理 1/a % m == a^m-2%m , 其中m要为素数。
所以：(b/a)%m == (b*a^m-2)%m。(快速幂)
函数代码如下：

int work(int a, int m)//1/a取模
{
    int ans = 1;
    int t=m-2;
    while(t)
    {
        if(t & 1)
          ans=ans* a % m;
        a=a*a%mod;
        t>>=1;
    }
    return ans;
}

2.转移方程：

dp[i][j]=(f[i+1][j-1]*t1*t2   +   f[i][j-1]*t2*t3   +   1 )%mod;

其中

t1 =(13-i*2)*3
t2=j-t1
t3=work(j,mod)
代码如下：

#include<bits/stdc++.h>
using namespace std;
#define endl '\n'
typedef long long ll;
ll dp[30][211];
const int mod=1e9+7;
int T;
ll s,r;
string S;
map<string,int>e;
ll work(ll a, ll m)
{
    ll ans = 1;
   ll t=m-2;
    while(t)
    {
        if(t & 1)
          ans=ans*a%mod;
        a=a*a%mod;
        t>>=1;
    }
    return ans;
}

void solve()
{
	cin>>S;
	e.clear();
	s=0;
	for(int i=0;i<S.length();i+=2)
	{
		if(e[S.substr(i,2)])
		  s++;
		else
		  e[S.substr(i,2)]=1;
	}
}
void init()
{
	for(int i=6;i>=0;i--)
	{
		ll t1=(13-i*2)*3;
		for(int j=t1;j<=123;j++)
		{
			ll t3=j-t1;
			ll t2=work(j,mod); 
			dp[i][j]=(dp[i+1][j-1]%mod*t1%mod*t2%mod +dp[i][j-1]%mod*t2%mod*t3%mod + 1 )%mod; // 后面数据很大，每次相乘都要取模一次
		}
	}
}
int main()
{
	init();
	cin>>T;
	for(int i=1;i<=T;i++)
	{
		solve();
		cout<<"Case #"<<i<<": "<<dp[s][123]<<endl;
	}
	
}