2020.icpc 南京

1. Evil Coordinate 构造+模拟

大意：给定长度为 n 的由 “LRUD” 组成的字符串，以及一个障碍点坐标，问任意排列字符串是否能避开障碍点，若可以避开，输出排列后的字符串。

思路：分类讨论的写法：大致构造角度，对于终点在障碍点的没办法，对于终点不在障碍点的

若L,R 均存在，若不能走到与障碍点x的坐标相等，那么显然可以先把L,R 输出完，然后对于U,D随便输出就行，对于L,R 恰好走到与 x 的坐标相等，我们可以随便拿出一个 L或 R ,然后类比于上面，最后再把拿出来的那个字符输出就行了，因为排出了终点在障碍点处的情况，最后一步一定不会走到障碍点，其他情况与之类似，细节非常多。

通过上面的分析，其实我们就是把L,R,U,D 相同的字符全部放一块，我们可以通过全排列确定L,R,U,D 的执行顺序，然后模拟就行了，只要有一种情况可行就是可行。

代码如下：

代码鸽了

2. Fireworks 期望、数学

大意：制作一个烟花花费 n 分钟，点燃之前所有做的烟花花费m 分钟，每只烟花是完美烟花的概率是 p, 当点燃的烟花至少有一只是完美烟花时就能停止制作烟花。求期望时间的最小值。

思路：之前没做过这类的题，思路完全错了。

我们先假假设每做 k 只烟花去点燃是最优。接下来实际上就是轮次重复试验。每轮成功的概率，即k只烟花中至少有一个完美的概率 $1-(1-p{)}^k$ ，符合几何分布。根据几何分布期望公式，即 期望=1/成功概率 = $\frac{1}{1-(1-p{)}^k}$ 。再乘上时间 kn+m。所以最终表达式为 $\frac{kn+m}{1-(1-p{)}^k}$。

通过求二阶导，发现是凹函数，三分求最小值就行了。

代码如下：

#include <bits/stdc++.h>
#define int long long 
#define rep(i,bbb,eee) for(int i=bbb;i<=eee;i++)
#define frep(i,bbb,eee) for(int i=bbb;i>=eee;i--)
#define mem(a,b) memset(a,b,sizeof(a))
#define cout(xxx) cout<<fixed<<setprecision(xxx)
#define inf 0x3f3f3f3f
#define pb push_back
#define AC signed
#define fi first
#define se second
using namespace std;
typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;
const int N=110,M=998244353;

int qmi(int a,int b)
{
	int res=1;
	while(b)
	{
		if(b&1)res=res*a%M;
		a=a*a%M;
		b>>=1;
	}
	return res%M;
}
inline int read()
{
    int x=0,f=1;
    char ch=getchar();
    while(!isdigit(ch)&&ch!='-')ch=getchar();
    if(ch=='-')f=-1,ch=getchar();
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return f*x;
}
inline void print(int x)
{   
   if(x<0){putchar('-');x=-x;}
   if(x>9) print(x/10);
   putchar(x%10+'0');
}
//____________________________________________//
int n,m;
double p;
double cale(int x)
{
	return (1.0*x*n+m)/(1.0-pow(1.0-p,x));
}
void solve()
{
	cin>>n>>m>>p;
	p*=0.0001;
	int l=1,r=1e6;
	while(l<r)
	{
		int lmid=l+(r-l)/3,rmid=r-(r-l)/3;
		if(cale(lmid)<cale(rmid))r=rmid-1;
		else l=lmid+1;
	}
	cout(10)<<cale(l)<<"\n";
}
AC main()
{
	ios::sync_with_stdio(false);cin.tie(0);
	int _=1;
	cin>>_;
	while(_--)solve();
	return 0;
}

3. Harmonious Rectangle 思维+打表

大意：给定一个 n*m 的矩阵，每个格子可以三种颜色，问使得存在最少一个 >=2乘2的子矩阵，该矩阵的四个角，相邻两个角的颜色相同。求方案数。

思路：显然对于 n=1||m=1 的必然是0.

对于其他的情况，我们假设任取两行，

a1 a2 a3 a4 …

b1 b2 b3 b4

a1,b1 的颜色组合方式 最多 3乘3，9种，那么一旦超过9列，必然会有一列颜色完全一样，所以颜色可以随便取。对于列同理。

所以 当 n>9||m>9 时，方案数为 $3^{nm}$

对于 n<=9&&m<=9 的，直接暴力打表就行了。

代码如下：

#include <bits/stdc++.h>
#define int long long 
#define rep(i,bbb,eee) for(int i=bbb;i<=eee;i++)
#define frep(i,bbb,eee) for(int i=bbb;i>=eee;i--)
#define mem(a,b) memset(a,b,sizeof(a))
#define cout(xxx) cout<<fixed<<setprecision(xxx)
#define inf 0x3f3f3f3f
#define pb push_back
#define AC signed
#define fi first
#define se second
using namespace std;
typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;
const int N=200010,M=1e9+7;
int qmi(int a,int b)
{
	int res=1;
	while(b)
	{
		if(b&1)res=res*a%M;
		a=a*a%M;
		b>>=1;
	}
	return res%M;
}
inline int read()
{
    int x=0,f=1;
    char ch=getchar();
    while(!isdigit(ch)&&ch!='-')ch=getchar();
    if(ch=='-')f=-1,ch=getchar();
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return f*x;
}
inline void print(int x)
{   
   if(x<0){putchar('-');x=-x;}
   if(x>9) print(x/10);
   putchar(x%10+'0');
}
//____________________________________________//
int n,m;
int del[9][9]=
{
    0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 15, 339, 4761, 52929, 517761, 4767849, 43046721, 387420489, 
    0, 339, 16485, 518265, 14321907, 387406809, 460338013, 429534507, 597431612, 
    0, 4761, 518265, 43022385, 486780060, 429534507, 792294829, 175880701, 246336683, 
    0, 52929, 14321907, 486780060, 288599194, 130653412, 748778899, 953271190, 644897553, 
    0, 517761, 387406809, 429534507, 130653412, 246336683, 579440654, 412233812, 518446848, 
    0, 4767849, 460338013, 792294829, 748778899, 579440654, 236701429, 666021604, 589237756, 
    0, 43046721, 429534507, 175880701, 953271190, 412233812, 666021604, 767713261, 966670169, 
    0, 387420489, 597431612, 246336683, 644897553, 518446848, 589237756, 966670169, 968803245
};
void solve()
{
	cin>>n>>m;
	if(n<=1||m<=1)cout<<"0\n";
	else if(n>9||m>9)cout<<qmi(3,n*m)<<"\n";
	else cout<<del[n-1][m-1]%M<<"\n";
}
AC main()
{
	ios::sync_with_stdio(false);cin.tie(0);
	int _=1;
	cin>>_;
	while(_--)solve();
	return 0;
}