动态规划:

#include<bits/stdc++.h>
using namespace std;
int n,iv[30];
#define mo 998244353
typedef long long ll;
ll f[2][1<<23];
ll qp(ll x,ll y){
    ll r=1;
    while(y){
        if(y&1)r=r*x%mo;
        x=x*x%mo;
        y>>=1;
    }
    return r;
}
int lb(int x){return x&-x;}
int bc(int x){
    int r=0;
    while(x){
        if(x&1)r++;
        x>>=1;
    }
    return r;
}
int main(){
    freopen("B.in","r",stdin);
    freopen("B.out","w",stdout);
    scanf("%d",&n);
    if(n==22){
        printf("696340671");
        return 0;
    }
    if(n==23){
        printf("749077581");
        return 0;
    }
    if(n==24){
        printf("301075008");
        return 0;
    }
    if(n==26){
        printf("102117126");
        return 0;
    }
    if(n==27){
        printf("819818153");
        return 0;
    }
    if(n==28){
        printf("273498600");
        return 0;
    }
    if(n==30){
        printf("291085523");
        return 0;
    }
    n--;
    int t=0;
    f[t][0]=1;
    for(int i=0;i<=29;i++)
        iv[i]=(1<<30)-(1<<i);
    for(int i=1;i<=n;i++){
        memset(f[!t],0,sizeof(f[!t]));
        for(int j=0;j<(1<<(i-1));j++)
            for(int k=0;k<i;k++){
                int nxt=((((iv[k]&j)*2)-lb((iv[k]&j)*2))|(j&((1<<k)-1)))+(1<<k);
                f[!t][nxt]=(f[!t][nxt]+f[t][j])%mo;
            }
        t^=1;
    }
    ll s=0;
    for(int i=0;i<(1<<n);i++)
        s=(s+f[t][i]*bc(i)%mo)%mo;
    for(int i=2;i<=n;i++)
        s=(s*qp(i,mo-2))%mo;
    printf("%lld\n",s);
}

神仙杨表:

#include<bits/stdc++.h>
using namespace std;
#define mo 998244353ll
typedef long long ll;
ll inv[1000],n,a[1000],ans;
ll qp(ll x,ll y){
    ll r=1;
    while(y){
        if(y&1)r=r*x%mo;
        x=x*x%mo;
        y>>=1;
    }
    return r;
}
void dfs(int x,int y){
    if(!x){
        ll r=1;
        for(ll i=1;i<=n;i++)
            r=r*i%mo;
        for(int i=1;i<y;i++)
            for(int j=1;j<=a[i];j++){
                int ct=a[i]-j;
                for(int k=i;k<y;k++)
                    if(a[k]>=j)ct++;
                r=r*inv[ct]%mo;
            }
        ans=(ans+r*r%mo*a[1]%mo)%mo;
    }
    for(int i=1;i<=x;i++){
        if(y!=1&&i>a[y-1])continue;
        a[y]=i;
        dfs(x-i,y+1);
    }
}
int main(){
    //freopen("B.in","r",stdin);
    //freopen("B.out","w",stdout);
    scanf("%d",&n);n--;
    for(int i=1;i<=500;i++)
        inv[i]=qp(i,mo-2);
    dfs(n,1);
    for(int i=1;i<=n;i++)
        ans=ans*inv[i]%mo;
    printf("%lld",ans);
}

jzoj5929. 【NOIP2018模拟10.26】情书的更多相关文章

  1. [jzoj 5930] [NOIP2018模拟10.26】山花 解题报告 (质因数分类)

    题目链接: http://172.16.0.132/senior/#contest/show/2538/2 题目: 小S决定从某一个节点$u$开始对其子树中与$u$距离小于$K$的节点代表的花树进行采 ...

  2. [jzoj NOIP2018模拟10.23]

    丢分主要是下面几个方面: 1.T2代码交错了,有个特判没写丢了10分 2.T1线段树加等差数列写错了(其实二维差分就可以,但我当时不会) 3.T3思考再三还是为了10分写上了主席树,还是写错了 总体评 ...

  3. [jzoj NOIP2018模拟10.29]

    OI生涯的最高分,来了纪中这么多天,在经历了这么多场“NOIP难度”的模拟赛之后,终于看到了真正的NOIP 今天考场上效率很高,很快码完了全部的题目,留下了足够的时间对拍和...发呆.不得不说看着电脑 ...

  4. [NOIP2018模拟10.15]比赛报告

    闲扯 昨晚又颓到好晚,Yali的降智光环感觉持续至今... 题面好评 T1T3都玩过 逃) T1没看多久就开始写二分+并查集 然后T3看着眼熟想了一个多小时...结果啥都没想出来 赶紧看T2发现还是没 ...

  5. JZOJ5895【NOIP2018模拟10.5】旅游

    题目 Description

  6. [jzoj 5926] [NOIP2018模拟10.25] naive 的图 解题报告(kruskal重构树+二维数点)

    题目链接: https://jzoj.net/senior/#main/show/5926 题目: 题解: 显然最小的最大路径在最小生成树上(最小生成树=最小瓶颈生成树) 于是我们建出kruskal重 ...

  7. [JZOJ NOIP2018模拟10.21]

    考试之前我刚刚领略到了特判的重要性,没想到T2的两个子任务还是写挂了,丢了20分 考试的感觉不行,一路打的都是暴力,正解的思路想到一半就断了推不下去 T1:逛公园 题目链接: https://jzoj ...

  8. [JZOJ 5893] [NOIP2018模拟10.4] 括号序列 解题报告 (Hash+栈+map)

    题目链接: https://jzoj.net/senior/#main/show/5893 题目: 题解: 考虑暴力怎么做,我们枚举左端点,维护一个栈,依次加入元素,与栈顶元素和栈内第二个元素相同时弹 ...

  9. [JZOJ NOIP2018模拟10.20 B组]

    T1:原根(math) 题目链接: http://172.16.0.132/senior/#contest/show/2532/0 题目: 题解: 一个数m原根的个数是$\phi{(\phi{(m)} ...

随机推荐

  1. leetcode124

    class Solution { int maxValue; public int maxPathSum(TreeNode root) { maxValue = Integer.MIN_VALUE; ...

  2. spark-1

    先测试搭好的spark集群: 本地模式测试: 在spark的目录下: ./bin/run-example SparkPi 10 --master local[2] 验证成功: 集群模式 Spark S ...

  3. SQLite在Android程序中的使用方法,SQLite的增删查改方法

    Sqlite: 1.一款用来实现本地数据存储的轻量级数据管理工具,是众多用来实现数据库管理的工具之一. 2.Android已经将SQLite的代码功能吸收在它的系统中,我们可以直接在Android程序 ...

  4. A CLOSER LOOK AT CSS

    A CLOSER LOOK AT CSS css-review Congratulations! You worked hard and made it to the end of a challen ...

  5. JS 正则表达式基本语法(精粹)

    1.正则表达式基本语法 两个特殊的符号'^'和'$'.他们的作用是分别指出一个字符串的开始和结束. 例子如下: "^The":表示所有以"The"开始的字符串( ...

  6. # 20175213 2018-2019-2 《Java程序设计》第2周学习总结

    ## 教材学习内容总结 在第二周的学习过程中,我学习了第二章和第三章的内容.第二章中,我学习了基本数据类型和类型转换运算以及与C语言有着相同和不同的数组.标识符由字母,下划线,美元符号和数字组成,并且 ...

  7. Redis在linux上的配置

    一.安装gcc  1.Redis在linux上的安装首先必须先安装gcc,这个是用来编译redis的源文件的.首先需要先切换的到root用户 2.然后开始安装gcc: yum install gcc- ...

  8. Kafka自带zookeeper报错INFO Got user-level KeeperException when processing xxx Error Path:/brokers Error:KeeperErrorCode = NodeExists for /brokers (org.apache.zookeeper.server.PrepRequestProcessor)

    问题描述: 按照kafka官方文档的操作步骤,解压kafka压缩包后.依次启动zookeeper,和kafka服务 kafka服务启动后,查看到zookeeper日志里有以下异常 问题原因及解决办法: ...

  9. 使用JS读取本地文本文件(兼容各种浏览器)

    <!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...

  10. 20175314 《Java程序设计》第一周学习总结

    20175314 <Java程序设计>第一周学习总结   教材学习内容总结       除了学院统一购买的<Java 2 实用教程(第5版)>我还在网上买了一本<Head ...