Wannafly挑战赛18

Wannafly挑战赛18


A. 序列

先考虑暴力,相邻两个树之间乘上给定的三种数,递推出下一个位置填什么,然后再check一下,最后一位是否为1即可。这样时间显然不行,但是给我们一种思路,就是中间的转换关系,确定唯一一个序列。现在的目标是让最后一位出现1,可以如果不管1,由-2和0.5取凑出1需要两个-2和两个0.5。那所有的转换中,就只要保证有若干组(-2,-2,0.5,0.5)存在,其他地方为1即可。具体公式见代码

#include <bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<=b;++i)
#define pb push_back
#define mem(W) memset(W,0,sizeof(W))
typedef long long ll;
inline int read() {
    char c=getchar();int x=0,f=1;
    while(!isdigit(c)){if(f==‘-‘)f=-1;c=getchar();}
    while(isdigit(c)){x=x*10+c-‘0‘;c=getchar();}
    return x*f;
}
using namespace std;
double a[1111];
int n;
const int mod = 1e9 + 7;
ll ans,c[1111][1111];
ll C(int n,int m){
    return c[n][m];
}
int main() {
    scanf("%d",&n);
    rep(i,0,n)c[i][0]=c[i][i]=1;
    rep(i,1,n)rep(j,1,i){
        c[i][j]=(c[i-1][j-1]+c[i-1][j])%mod;
    }
    rep(i,0,(n-1)/4){
        ans = (ans + (C(n-1,i*4)*C(4*i,2*i))%mod)%mod;
    }
    printf("%lld\n",ans);
    return 0;
}

B. 随机数

设\(P_n\)是前n次出现奇数次的概率,单次出现1的概率\(X = \frac{a}{10000}\), 那么有

\[P_n=(1-X)P_{n-1}+X(1-P_{n-1})=(1-2X)P_{n-1}+X, P_0=0\]

由中学数学可知,通项为:\(P_n=\frac{1-(1-2X)^n}{2}\)

之后发现n很大,所以我胆子很大的用java写了快速幂......单次运算复杂度太高TLE,实际上碰见这种指数项远大于模数的情况第一个就该想到循环节,对于mod为素数(mod-1)就是循环节,然后就可以把n降低,直接套公式即可。

#include <bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<=b;++i)
#define pb push_back
#define mem(W) memset(W,0,sizeof(W))
typedef long long ll;
inline int read() {
    char c=getchar();int x=0,f=1;
    while(!isdigit(c)){if(f==‘-‘)f=-1;c=getchar();}
    while(isdigit(c)){x=x*10+c-‘0‘;c=getchar();}
    return x*f;
}
using namespace std;
ll n,a;
const int mod = 1e9 + 7;
ll q_pow(ll a,ll b){
    ll ans=1;
    while(b){
        if(b&1)ans=(ans*a)%mod;
        a=(a*a)%mod;
        b>>=1;
    }
    return ans;
}
char s[1000100];
int main() {
    scanf(" %lld %s",&a,s);
    int len=strlen(s);
    rep(i,0,len-1)n=(n*10+s[i]-‘0‘)%(mod-1);
    ll X = (a*q_pow(10000LL,mod-2))%mod;
    ll ans = (1-q_pow(1-2*X,n)+mod)%mod*q_pow(2,mod-2)%mod;
    ans=(ans+mod)%mod;
    printf("%lld\n",ans);
    return 0;
}

C. 异或和

思考如何快速求出一个位置的期望,曼哈顿距离常见套路,把x,y拆开考虑,那么把x,y都投影到坐标轴上,分别预处理每个位置左右两边的1到它的距离和,o(n)预处理,然后枚举一下就解决了。

#include <bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<=b;++i)
#define frep(i,a,b) for(int i=a;i>=b;--i)
#define pb push_back
#define mem(W) memset(W,0,sizeof(W))
typedef long long ll;
inline int read() {
    char c=getchar();int x=0,f=1;
    while(!isdigit(c)){if(f==‘-‘)f=-1;c=getchar();}
    while(isdigit(c)){x=x*10+c-‘0‘;c=getchar();}
    return x*f;
}
using namespace std;
const int N = 2000+50;
int n,m;
const int mod = 1e9 + 7;
ll q_pow(ll a,ll b){
    ll ans=1;
    while(b){
        if(b&1)ans=(ans*a)%mod;
        a=(a*a)%mod;
        b>>=1;
    }
    return ans;
}
char s[2002];
ll num,X[N],Y[N],LX[N],RX[N],LY[N],RY[N];
int main() {
    scanf("%d%d",&n,&m);
    rep(i,1,n){
        scanf(" %s",s+1);
        rep(j,1,m)if(s[j]==‘1‘){
            ++X[i];
            ++Y[j];
            ++num;
        }
    }
    int t=X[1];
    rep(i,2,n){
        LX[i] = (LX[i-1] + t)%mod;
        t=(t+X[i])%mod;
    }
    t=Y[1];
    rep(i,2,m){
        LY[i]=(LY[i-1]+t)%mod;
        t=(t+Y[i])%mod;
    }
    t=X[n];
    frep(i,n-1,1){
        RX[i] = (RX[i+1] + t)%mod;
        t=(t+X[i])%mod;
    }
    t=Y[m];
    frep(i,m-1,1){
        RY[i] = RY[i+1] + t;
        t=(t+Y[i])%mod;
    }
    ll ans=0;
    rep(i,1,n)rep(j,1,m){
        ll tmp = (((RX[i]+LX[i]+RY[j]+LY[j])%mod*q_pow(num,mod-2))%mod)%mod;
        ans^=tmp;
    }
    printf("%lld\n",ans);
    return 0;
}

原文地址:https://www.cnblogs.com/RRRR-wys/p/9216028.html

时间: 2024-10-10 00:06:35

Wannafly挑战赛18的相关文章

Wannafly挑战赛3

Wannafly挑战赛3 A    珂朵莉 B    遇见 水题 #include<bits/stdc++.h> using namespace std; #pragma comment(linker, "/STACK:102400000,102400000") #define rep(i,a,b) for (int i=a; i<=b; ++i) #define per(i,b,a) for (int i=b; i>=a; --i) #define mes(a

【Wannafly挑战赛4】F 线路规划 倍增+Kruskal+归并

[Wannafly挑战赛4]F 线路规划 题目描述 Q国的监察院是一个神秘的组织.这个组织掌握了整个帝国的地下力量,监察着Q国的每一个人.监察院一共有N个成员,每一个成员都有且仅有1个直接上司,而他只听从其上直接司的命令.其中1号成员是监察院的院长,这个庞然大物的主人.由于时代的进步,监察院议会决定升级组织的旧式通信器,安装最新的反侦测通信器.他们拿出了M组线路方案,其中第i组线路方案可以用一个四元组(x[i].y[i].k[i].w[i])描述,表示第x[i]号成员可以安装与y[i]号成员的直

牛客网 Wannafly挑战赛8 C-小C打比赛 (状压DP)

小C现在要参加一场wannafly挑战赛,一场挑战赛一共有n道题,一共有m分钟. 对于第i道题,小C解决它需要恰好j分钟的概率是pi,j. 小C每次会选择某一道没做完的题,然后把它解决(不能中途放弃),之后再决策下一道要做的题是哪道. 求小C在最优策略下,期望能做出几道题. 输入描述: 第一行两个正整数n,m接下来一共n行,每行有m个小数,第i行的第j个小数表示p i,j (这里假设不存在0分钟A题的dalao). 输出描述: 输出一个小数,表示期望能做出几道题,保留小数点后五位. 示例1 输入

Wannafly挑战赛11 D 白兔的字符串 Hash

Wannafly挑战赛11 D   白兔的字符串 白兔有一个字符串T.白云有若干个字符串S1,S2..Sn. 白兔想知道,对于白云的每一个字符串,它有多少个子串是和T循环同构的. 提示:对于一个字符串a,每次把a的第一个字符移动到最后一个,如果操作若干次后能够得到字符串b,则a和b循环同构. 所有字符都是小写英文字母 输入描述: 第一行一个字符串T(|T|<=10^6)第二行一个正整数n (n<=1000)接下来n行为S1~Sn (|S1|+|S2|+…+|Sn|<=10^7),max(

Wannafly挑战赛19

Wannafly挑战赛19 A. 队列Q 需要支持把一个元素移到队首,把一个元素移到队尾,移到队首就直接放到队首前面那个位置,原位置标为0,队尾同理. #include <bits/stdc++.h> #define rep(i,a,b) for(int i=a;i<=b;++i) typedef long long ll; const int N = 30000200; using namespace std; int n,m; int q[N],hd,x,ed,P[N]; char

Wannafly挑战赛22游记

Wannafly挑战赛22游记 幸福的人都是相似的,不幸的人各有各的不幸. --题记 A-计数器 题目大意: 有一个计数器,计数器的初始值为\(0\),每次操作你可以把计数器的值加上\(a_1,a_2,\ldots,a_n\)中的任意一个整数,操作次数不限(可以为\(0\)次),问计数器的值对\(m\)取模后有几种可能. 思路: 由裴蜀定理易得,答案即为\(\frac m{\gcd(m,a_1,a_2,\ldots,a_n)}\). 源代码: #include<cstdio> #include

Wannafly挑战赛25游记

Wannafly挑战赛25游记 A - 因子 题目大意: 令\(x=n!(n\le10^{12})\),给定一大于\(1\)的正整数\(p(p\le10000)\)求一个\(k\)使得\(p^k|x\)并且\(p^{k+1}\not|x\)的因子. 思路: 枚举\(p\)的每一个质因数\(q\),求出它在\(n!\)出现次数\(/p\)中出现次数,取\(\min\)即可.对于一个质因数\(q\),在\(n!\)中出现的次数等于\(\sum_{i=1}^{\inf}\frac n{q^i}\).

Wannafly挑战赛3 record

B 遇见 时间限制:C/C++ 1秒,其他语言2秒空间限制:C/C++ 65536K,其他语言131072K64bit IO Format: %lld题目描述A和B在同一条路上,他们之间的距离为 k 米.A现在想见到B,所以A开车以 x km/h的速度朝着B的方向行驶,同时B也以 y km/h的速度朝着A的方向走去.A的车有 n 个档位,每个档位有不同的速度.现在假设A开车去见B,求他最快和最慢在几秒后能见到B.输入描述:一开始一行三个整数 n, m, k ,代表A的车的档位数.B行走的速度和A

Wannafly挑战赛4. B

题意:求子区间异或和,要求区间长度在l到r之间,并且为偶数 题解:对于每一位算贡献,可以分奇偶来记录,计算的时候只加上奇偶性相同的就保证了为偶数,从大于l的点开始每次++,从大于r的点每次--,记录二进制上所有权值和 代码: 1 #include<bits/stdc++.h> 2 #define db double 3 #define ll long long 4 #define vec vector<ll> 5 #define Mt vector<vec> 6 #de