4.9 省选模拟赛 生成树求和 变元矩阵树定理 生成函数 iDFT 插值法
有同学在loj上找到了加强版 所以这道题是可以交的。LINK:生成树求和 加强版
对于30分 爆搜 可实际上我爆搜只过了25分 有同学使用按秩合并并茶几的及时剪枝通过了30分。
const int MAXN=45;
int n,m;
struct wy
{
int x,y,z;
}t[MAXN];
int w[MAXN],f[MAXN];ll ans;
inline int getfather(int x){return x==f[x]?x:getfather(f[x]);}
inline int add(int x,int y)
{
int cnt=0;
int p=1;
while(x||y)
{
int w1=x%3;
int w2=y%3;
int cc=(w1+w2)%3;
cnt+=cc*p;
p=p*3;x/=3;y/=3;
}
return cnt;
}
inline void dfs(int x,int v)
{
if(m-x+1<n-v)return;
if(v==n)
{
int cnt=0;
rep(1,n,i)f[i]=i;
rep(1,n-1,i)
{
int xx=getfather(t[w[i]].x);
int yy=getfather(t[w[i]].y);
if(xx==yy)return;
f[xx]=yy;
cnt=add(cnt,t[w[i]].z);
}
ans=(ans+cnt)%mod;
return;
}
w[v]=x;
dfs(x+1,v+1);
dfs(x+1,v);
}
int main()
{
freopen("sum.in","r",stdin);
freopen("sum.out","w",stdout);
get(n);get(m);
rep(1,m,i)
{
int get(x);int get(y);int get(z);
t[i]=(wy){x,y,z};
}
if(m<=35){dfs(1,1);putl(ans);}
return 0;
}
之所以这样做是因为虽然这道题很容易想到矩阵树定理,按位处理 但是依然无法做出来。
矩阵树除了能求方案还可以求边上的权值积的和,即变元矩阵树定理 具体理解可以把边权想象成有边权条重边。
但是对于这道题来说还是不行 求的是生成树边权和的和 把积变和只有是在指数上了。
而且还有三进制的加法 非常的不可做 但是可以按位处理。
把边权变到指数上只能是 生成函数了 列出生成函数 每一条边的权值用生成函数表示。
利用高斯消元解行列式 最后求出的多项式就好了。
而直接使用多项式无法方便的求出行列式 中间由于矩阵树定理还无法取模 所以有多项式乘法除法啥的 所以变的更不可做了。
但是考虑最后的次数 其实并不大最多2n次 可以取2n+1个数带入 然后利用插值求解原多项式。
复杂度\(n^4log\)
为了通过加强版 我们考虑一种比较简单的方式 如果能找到\(x^3=1(mod P)\)
那么我么最后的多项式就只有三项且是正确的了 考虑到复平面单位根。
由于最后有三项 所以 把复平面平分三份的复平面单位根即可。
考虑最后的插值 利用复平面单位根的相加相乘操作可以轻松的插值出来。不需要再次高斯消元或者拉格朗日插值。
非常的巧妙。
const ll MAXN=100010,inv2=(mod+1)>>1,inv3=(mod+1)/3,g3=82062379,maxn=110;
inline ll ksm(ll b,ll p,ll cc){ll cnt=1;while(p){if(p&1)cnt=cnt*b%cc;p=p>>1;b=b*b%cc;}return cnt;}
struct wy
{
ll a,b;
inline wy(ll aa=0,ll bb=0){a=aa;b=bb;}
inline wy operator +(wy w){return wy((a+w.a)%mod,(b+w.b)%mod);}
inline wy operator -(wy w){return wy((a-w.a+mod)%mod,(b-w.b+mod)%mod);}
inline wy operator *(ll y){return wy(a*y%mod,b*y%mod);}
inline wy operator *(wy y){return wy((a*y.a%mod-b*y.b%mod+mod)%mod,(a*y.b+b*y.a)%mod);}
inline wy operator /(ll y)
{
ll inv=ksm(y,mod-2,mod);
return wy(a*inv%mod,b*inv%mod);
}
inline wy operator /(wy y)
{
return wy(a,b)*wy(y.a,mod-y.b)/((y.a*y.a+y.b*y.b)%mod);
}
}a[maxn][maxn],cs[3],inv,w[3],invw[3];
int pow3[maxn],eu[maxn*maxn],ev[maxn*maxn],ec[maxn*maxn];
ll n,m;
inline wy det()
{
wy ret=wy(1,0);
rep(1,n-1,i)
{
ll p=i;
rep(i,n-1,j)if(a[j][i].a||a[j][i].b){p=j;break;}
if(i!=p)
{
rep(1,n-1,k)swap(a[i][k],a[p][k]);
ret=ret*(mod-1);
}
rep(i+1,n-1,j)
{
inv=a[j][i]/a[i][i];
rep(i,n-1,k)a[j][k]=a[j][k]-inv*a[i][k];
}
ret=ret*a[i][i];
}
return ret;
}
inline void idft(wy *p)
{
wy c[3];
c[0]=p[0];c[1]=p[1];c[2]=p[2];
p[0]=c[0]+c[1]+c[2];
p[1]=c[0]+invw[1]*c[1]+c[2]*invw[2];
p[2]=c[0]+invw[2]*c[1]+c[2]*invw[1];
rep(0,2,i)p[i]=p[i]*inv3;
}
inline ll calc(ll p)
{
wy w0,w1,w2,c;
rep(0,2,i)
{
memset(a,0,sizeof(a));
w0=wy(1,0);
w1=w[i];
w2=w[(i+i)%3];
rep(1,m,j)
{
ll k=ec[j]/pow3[p]%3;
c=(k==1?w1:k==2?w2:w0);
a[eu[j]][eu[j]]=a[eu[j]][eu[j]]+c;
a[ev[j]][ev[j]]=a[ev[j]][ev[j]]+c;
a[eu[j]][ev[j]]=a[eu[j]][ev[j]]-c;
a[ev[j]][eu[j]]=a[ev[j]][eu[j]]-c;
}
cs[i]=det();
}
idft(cs);
return (cs[1].a+cs[2].a*2)*pow3[p]%mod;
}
signed main()
{
freopen("1.in","r",stdin);
//freopen("sum.out","w",stdout);
get(n);get(m);ll ans=0;
w[0]=wy(1,0);w[1]=wy(mod-inv2,g3*inv2%mod);
w[2]=wy(mod-inv2,mod-g3*inv2%mod);
invw[0]=w[0]/w[0];
invw[1]=w[0]/w[1];
invw[2]=w[0]/w[2];
rep(1,m,i)get(eu[i]),get(ev[i]),get(ec[i]);
pow3[0]=1;
rep(1,10,i)pow3[i]=pow3[i-1]*3;
rep(0,10,i)ans=(ans+calc(i))%mod;
putl((ans+mod)%mod);return 0;
}
4.9 省选模拟赛 生成树求和 变元矩阵树定理 生成函数 iDFT 插值法的更多相关文章
- P6624-[省选联考2020A卷]作业题【矩阵树定理,欧拉反演】
正题 题目链接:https://www.luogu.com.cn/problem/P6624 题目大意 \(n\)个点的一张图,每条边有权值,一棵生成树的权值是所有边权和乘上边权的\(gcd\),即 ...
- 省选模拟赛 4.26 T1 dp 线段树优化dp
LINK:T1 算是一道中档题 考试的时候脑残了 不仅没写优化 连暴力都打挂了. 容易发现一个性质 那就是同一格子不会被两种以上的颜色染.(颜色就三种. 通过这个性质就可以进行dp了.先按照左端点排序 ...
- 6.6 省选模拟赛 线段 二维数点问题 树套树 CDQ分治
LINK:线段 还是太菜了 没看出这道题真正的模型 我真是一个典型的没脑子选手. 考虑如何查询答案. 每次在一个线段x的状态被更改后 可以发现有影响的是 和x相连那段极长连续1子段. 设这个子段左端点 ...
- 5.21 省选模拟赛 luogu P4207 [NOI2005]月下柠檬树 解析几何 自适应辛普森积分法
LINK:月下柠檬树 之前感觉这道题很鬼畜 实际上 也就想到辛普森积分后就很好做了. 辛普森积分法的式子不再赘述 网上多的是.值得一提的是 这道题利用辛普森积分法的话就是一个解析几何的问题 而并非计算 ...
- 5.20 省选模拟赛 T1 图 启发式合并 线段树合并 染色计数问题
LINK:图 在说这道题之前吐槽一下今天的日子 520 = 1+1+4+514. /cy 这道题今天做的非常失败 一点分都没拿到手 关键是今天的T3 把我整个人给搞崩了. 先考虑 如果得到了这么一张图 ...
- 【洛谷比赛】[LnOI2019]长脖子鹿省选模拟赛 T1 题解
今天是[LnOI2019]长脖子鹿省选模拟赛的时间,小编表示考的不怎么样,改了半天也只会改第一题,那也先呈上题解吧. T1:P5248 [LnOI2019SP]快速多项式变换(FPT) 一看这题就很手 ...
- 省选模拟赛第四轮 B——O(n^4)->O(n^3)->O(n^2)
一 稍微转化一下,就是找所有和原树差距不超过k的不同构树的个数 一个挺trick的想法是: 由于矩阵树定理的行列式的值是把邻接矩阵数值看做边权的图的所有生成树的边权乘积之和 那么如果把不存在于原树中的 ...
- CJOJ 1943 【重庆八中模拟赛】寻找代表元(二分图最大匹配)
CJOJ 1943 [重庆八中模拟赛]寻找代表元(二分图最大匹配) Description 八中一共有n个社团,分别用1到n编号. 八中一共有m个人,分别用1到m编号.每个人可以参加一个或多个社团,也 ...
- loj6271 「长乐集训 2017 Day10」生成树求和 加强版(矩阵树定理,循环卷积)
loj6271 「长乐集训 2017 Day10」生成树求和 加强版(矩阵树定理,循环卷积) loj 题解时间 首先想到先分开三进制下每一位,然后每一位分别求结果为0,1,2的树的个数. 然后考虑矩阵 ...
随机推荐
- MVC + EFCore 项目实战 - 数仓管理系统3 - 完成整体样式风格配置
上次课程我们新建了管理员的模板页. 本次我们就完善这个模板页,顺便加入样式和一些基本的组件,配置好整个项目的UI风格. 一.引入 共用的css和js文件 后端库用nuget, 前端库用libman ...
- keyring源码加密解密函数分析
Encrypt the page data contents. Page type can't be FIL_PAGE_ENCRYPTED, FIL_PAGE_COMPRESSED_AND_ENCRY ...
- [ArcEngine二次开发]为Feature的属性赋值
在创建FeatureClass之后,需要为FeatureClass添加Features,在为Feature的字段赋值时,代码大致如下: 在这里赋值的时候,出现了一个错误: The operation ...
- 隐藏Layer中的Iframe内部元素
隐藏Layer中的Iframe内部元素 1.Layer:一个web弹出层/窗解决方案 2.隐藏Layer弹出层中Iframe内部的元素,例子如下: 这个就是在Layer弹出层内的内容,加载完毕之后,获 ...
- bzoj4631踩气球
bzoj4631踩气球 题意: 有一个序列和一个区间集合,每次将序列中的一个数-1,求此时集合里有多少个区间和为0.序列大小≤100000,区间数≤100000,操作数≤100000. 题解: 此题解 ...
- bzoj3211花神游历各国&&bzoj3038上帝造题的七分钟2*
bzoj3211花神游历各国 题意: n个数的序列,m个操作,操作两种:区间开根(向下取整)和区间求和.n≤100000,m≤200000,序列中的数非负且≤109. 题解: 一个≤109的数开6次根 ...
- bzoj2697特技飞行*
bzoj2697特技飞行 题意: N个单位时间,每个单位时间可以进行一项特技动作,可选的动作有K种,每种动作有一个刺激程度Ci.每次动作的价值为(距上次该动作的时间)*Ci,若为第一次进行该动作,价值 ...
- commons-fileload图片文件上传工具 , servlet文件图片上传案列
本案列是java maven工程小项目,提供个大家学习! 1.在pom.xml文件中导入依赖: <!--文件上传依赖--><dependency> <groupId&g ...
- [Qt2D绘图]-05绘图设备-QPixmap&&QBitmap&&QImage&&QPicture
这篇笔记记录的是QPainterDevice(绘图设备,可以理解为一个画板) 大纲: 绘图设备相关的类:QPixmap QBitmap QImage QPicture QPixmap ...
- Web应用程序安全与风险
一.Web应用程序安全与风险 更多渗透测试相关内容请关注此地址:https://blog.csdn.net/weixin_45380284 1.web发展历程 静态内容阶段(HTML) CGI程序阶段 ...