Ngshily

摘要： 先写了一道傻逼题poj3041 1 #include 2 #include 3 #define maxn 505 4 int n,k,ma[maxn][maxn],used[maxn],bel[maxn]; 5 bool find(int x){ 6 for(int y=1;y 2 usin... 阅读全文

摘要： 我们计查分数列为$A$,$A[i]=a[i]-a[i-1]$对于每一个差分序列,由于$(k-1)*m$n-\sum{A[i]}$总共有$m^{k-1}$个差分数列,所以-->$ans=\sum...\sum\left(n-\sum{A[i]}\right)$ =$m^... 阅读全文

摘要： 精度BOOM，然而可以用对数的性质把乘法罗成加法 1 #include 2 using namespace std; 3 4 #define esp 1e-6 5 #define inf 0x3f3f3f3f 6 #define maxn 500005 7 int cnt,v[maxn<<1],n... 阅读全文

摘要： 对于一个盘子(dep[u]中，当且仅当 ①lca(u,v)==u时 设w为u的是v的祖先的儿子，则 1 2 using namespace std; 3 #define maxn 40005 4 #define maxt 32000005 5 6 int c... 阅读全文

摘要： 首先对于一张拓扑图，生成树的个数ans-->除根节点以外的节点的入度之积但是加了一条边以后就变得不拓扑了，所以要减去所产生的不合法的情况，即统计由Y到X的路径条数对于包含这样的路径的不合法的"树"，会有-->所有不包含于这条路径的点的入度的乘积种即-->ans/du[所有属于路径的节点]我们用f[i... 阅读全文

摘要： 很容易想到思路，然而考试的时候勺的呢，没想到用priority_qq倒序输出反图的字典序最大拓扑序 1 #include 2 using namespace std; 3 4 #define maxn 100005 5 #define maxm 1000... 阅读全文

摘要： 神一般的DP，f[i][j]表示前i张牌，还剩j轮的概率对于第i张，如果它在后j轮中都没有被出出去，则有f[i+1][j]r会比较快，233 1 #include 2 using namespace std; 3 double f[250][150],p[250]; 4 ... 阅读全文

摘要： 大半夜罗了奇怪的做法，发现Ｍ啦ＱＡＱ，内心死亡，然而第二天早上写完另一种发现交的代码是nn=4，内心死亡again 首先从题意分析就是，线段树上罗，然后求欧拉函数 两种做法的区别在于求的公式不同 第一种是，phi(x)=sigma(pn-1)pn^(an-1)，其中pn为对x分解质因数 ... 阅读全文

摘要： 首先想到dp的状态f[i][j]，当前处理到第i位，并有末尾j位是与原串匹配上的转移就应该是f[i][j]=sigma(f[i-1][k]，k->匹配到第k位，能加上一位使得匹配长度变成j位)；也可以写成f[i][j]=sigma(f[i-1][k]*a[k][j])，其中a[k][j]表示能使匹配... 阅读全文

摘要： 安利论文　《弦图与区间图》---cdq任意一个长度＞３的环都有弦（即非连接相邻点的边），满足弦图的性质，对弦图进行最小染色然而并不能想得明白，维萨最小染色是从从后往前，最大独立集是从前往后．．． 1 #include 2 using namespace std; 3 #define maxn 100... 阅读全文