题目链接

题意

给定一张边带权的无向图,求生成树的权值和是 k 的倍数的生成树个数模 p 的值。

\(n\leq 100,k\leq 100,p\mod k=1\)

Sol

看见整除然后 \(p\mod k=1\) ,那么可以套个单位根反演。

我们要求的东西就是:

\(\sum_{E}[k|(\sum_{e\in E}val_e)]\)

单位根反演一套:

\(\frac{1}{k}\sum_{E} \sum_{i=0}^{k-1} w_k^{(\sum_{e\in E}val_e)i}\)

然后又是常规操作:

\(\frac{1}{k}\sum_{i=0}^{k-1}\sum_{E} w_k^{(\sum_{e\in E}val_e)i}\)

\(\frac{1}{k}\sum_{i=0}^{k-1}\sum_{E} \prod_{e\in E} (w_k^{i})^{val_e}\)

把一条边的边权看作 \((w_k^{i})^{val_e}\) 矩阵树定理求一下就做完了。

code:

#include<bits/stdc++.h>
#define Set(a,b) memset(a,b,sizeof(a))
using namespace std;
int mod;
template <typename T> inline void init(T&x){
x=0;char ch=getchar();bool t=0;
for(;ch>'9'||ch<'0';ch=getchar()) if(ch=='-') t=1;
for(;ch>='0'&&ch<='9';ch=getchar()) x=(x<<1)+(x<<3)+(ch-48);
if(t) x=-x;return;
}
typedef long long ll;
template <typename T>inline void Inc(T&x,int y){x+=y;if(x>=mod) x-=mod;return;}
template <typename T>inline void Dec(T&x,int y){x-=y;if(x < 0) x+=mod;return;}
template <typename T>inline int fpow(int x,T k){int ret=1;for(;k;k>>=1,x=(ll)x*x%mod) if(k&1) ret=(ll)ret*x%mod;return ret;}
int Sum(int x,int y){x+=y;if(x>=mod) return x-mod;return x;}
int Dif(int x,int y){x-=y;if(x < 0 ) return x+mod;return x;}
const int N=101;
int n,m,k,p,g;
struct edge{
int u,v,c;
}E[N*N];
namespace Matrix_Tree{
int a[N][N];
inline void Build(int w){
Set(a,0);
for(int i=1;i<=m;++i) {
int u=E[i].u,v=E[i].v,c=E[i].c;
int val=fpow(w,c);
Dec(a[u][v],val),Dec(a[v][u],val);
Inc(a[u][u],val),Inc(a[v][v],val);
}return;
}
inline int Gauss(int n){
int f=0;
for(int i=1;i<=n;++i) {
int p=i;
for(int j=i;j<=n;++j) {if(a[i][j]) {p=j;break;}}
if(p!=i) f^=1,swap(a[p],a[i]);
int inv=fpow(a[i][i],mod-2);
for(int j=i+1;j<=n;++j){
if(!a[j][i]) continue;
int t=Dif(0,(ll)a[j][i]*inv%mod);
for(int k=i;k<=n;++k) Inc(a[j][k],(ll)a[i][k]*t%mod);
}
}
int ret=1;
for(int i=1;i<=n;++i) ret=(ll)ret*a[i][i]%mod;
if(f) ret=Dif(0,ret);return ret;
}
}
inline void Getroot(int mod){
int x=mod-1;static int pri[50],cnt=0;
for(int i=2;i*i<=x;++i) if(x%i==0) {pri[++cnt]=i,x/=i;while(x%i==0) x/=i;}
for(g=2;;++g){bool fl=1;
for(int i=1;i<=cnt;++i) if(fpow(g,(mod-1)/pri[i])==1) {fl=0;break;}
if(fl)return;
}
}
int main()
{
init(n),init(m),init(k),init(p);
mod=p;Getroot(mod);int u,v,c;
for(int i=1;i<=m;++i){init(u),init(v),init(c);E[i]=(edge){u,v,c};}
int W=fpow(g,(mod-1)/k);
int w=1,ans=0;
for(int i=0;i<k;++i,w=(ll)w*W%mod) {
Matrix_Tree::Build(w);
Inc(ans,Matrix_Tree::Gauss(n-1));
}
ans=(ll)ans*fpow(k,mod-2)%mod;
cout<<ans<<endl;
return 0;
}

【牛客Wannafly挑战赛23】F 计数的更多相关文章

  1. 牛客Wannafly挑战赛23 B.游戏

    游戏 题目描述 小N和小O在玩游戏.他们面前放了n堆石子,第i堆石子一开始有ci颗石头.他们轮流从某堆石子中取石子,不能不取.最后无法操作的人就输了这个游戏.但他们觉得这样玩太无聊了,更新了一下规则. ...

  2. 牛客wannafly 挑战赛14 B 前缀查询(trie树上dfs序+线段树)

    牛客wannafly 挑战赛14 B 前缀查询(trie树上dfs序+线段树) 链接:https://ac.nowcoder.com/acm/problem/15706 现在需要您来帮忙维护这个名册, ...

  3. 牛客网练习赛23 F 托米的游戏

    链接:https://www.nowcoder.com/acm/contest/156/F 来源:牛客网 题目描述 题目背景编不下去了 托米有一棵有根树 T, 树根为1,每轮他会在剩下的子树中等概率一 ...

  4. 牛客~~wannafly挑战赛19~A 队列

    链接:https://www.nowcoder.com/acm/contest/131/A来源:牛客网 题目描述 ZZT 创造了一个队列 Q.这个队列包含了 N 个元素,队列中的第 i 个元素用 Qi ...

  5. 牛客Wannafly挑战赛23F 计数(循环卷积+拉格朗日插值/单位根反演)

    传送门 直接的想法就是设 \(x^k\) 为边权,矩阵树定理一波后取出 \(x^{nk}\) 的系数即可 也就是求出模 \(x^k\) 意义下的循环卷积的常数项 考虑插值出最后多项式,类比 \(DFT ...

  6. 牛客 Wannafly 挑战赛26D 禁书目录 排列组合 概率期望

    原文链接https://www.cnblogs.com/zhouzhendong/p/9781060.html 题目传送门 - NowCoder Wannafly 26D 题意 放一放这一题原先的题面 ...

  7. 牛客Wannafly挑战赛26E 蚂蚁开会(树链剖分+线段树)

    传送门 题面描述 一颗n个节点的树,m次操作,有点权(该节点蚂蚁个数)和边权(相邻节点的距离). 三种操作: 操作1:1 i x将节点i的点权修改为x.(1 <= i <= n; 1 &l ...

  8. [牛客Wannafly挑战赛27D]绿魔法师

    description newcoder 给你一个空的可重集合\(S\). \(n\)次操作,每次操作给出\(x\),\(k\),\(p\),执行以下操作: \(opt\ 1\):在S中加入x. \( ...

  9. 牛客Wannafly挑战赛11E 白兔的刁难

    传送门 如果大力推单位根反演就可以获得一个 \(k^2logn\) 的好方法 \[ans_{t}=\frac{1}{k}\sum_{i=0}^{k-1}(w_k^{-t})^i(w_k^i+1)^n\ ...

随机推荐

  1. kaptcha Spring 整合

    转自:http://my.oschina.net/CandyDesire/blog/209364 生成验证码的方式有很多,个人认为较为灵活方便的是Kaptcha ,他是基于SimpleCaptcha的 ...

  2. 阶段3 1.Mybatis_04.自定义Mybatis框架基于注解开发_1 今日课程内容介绍

  3. MinGW GCC 9.1 2019年5月3日 出炉啦

    GNU 2019-05-03 发布了 GCC 9.1 https://gcc.gnu.org/onlinedocs/9.1.0/ 有详细的说明MinGW 上可用的 GCC 9.1 版本下载地址 [ m ...

  4. docker nginx ssl

    docker run -d -p 80:80 -p 443:443 -v $(pwd)/conf/conf.d:/etc/nginx/conf.d -v $(pwd)/cert:/etc/nginx/ ...

  5. 通过vue-router实现组件间的跳转

    三.通过VueRouter来实现组件之间的跳转提供了3种方式实现跳转:①直接修改地址栏中的路由地址 <!doctype html> <html> <head> &l ...

  6. Servlet 响应 响应相关与重定向 请求 获取表单数据2种方法

    一.HttpServletResponse  (响应) 包括下面三个: 1.响应消息行  HTTP/1.1  200 OK 200是HTTP状态码, 代表请求已成功. (查httpservletres ...

  7. 阿里云 Centos7 部署 Django 项目

    前期准备 阿里云服务器 mysql数据库 已经本地运行成功的项目 阿里云服务器的环境配置 Git #代码管理 Gitlab #代码托管,要求服务器内存不低于2G,我选择放弃 Mysql #连接数据库 ...

  8. fiddler过滤机制讲解

    1.User Fiters启用 2.Action Action:Run Filterset now是否运行,Load Filterset加载,Save Filterset保存: 3.Hosts过滤 Z ...

  9. 【HBase】四、HBase的安装及命令行接口

      通过前面的介绍,对HBase数据模型,运行机制等理论基本了解,接下来从实践的角度介绍HBase的安装以及其各种接口的使用方法. 1.HBase的安装   HBase安装很简单,和所有的Hadoop ...

  10. 初学node.js-nodejs中实现用户注册路由

    经过前面几次的学习,已经可以做下小功能,今天要实现的事用户注册路由. 一.users_model.js  功能:定义用户对象模型 var mongoose=require('mongoose'), S ...