Song

原创 ssh免密登陆

ssh使用公钥密钥方式登陆，既可以免去输入密码，又可以防止手动输入账户密码时密码被截获，既方便又安全。方法：在本机 生成公钥密钥ssh-keygen -t [rsa|dsa]rsa或dsa可选其一如果选择rsa方式，则在用户目录/.ssh/目录下会生成 id_rsa，id_rsa.pub两个文件。id_rsa为密钥，id_rsa.pub为公钥。将id_rsa.pub拷贝到目标服务...

原创 poj 1159 Palindrome LCS一维滚动数组优化

原题地址 http://poj.org/problem?id=1159最长公共子序列问题LCS：给出两个子序列A,B, 如 A B C B D A B B D C A B A 最长公共子序列为4: BCBA 递推公式： dp(i,j)为A,B当前长度为i,j的LCSdp(i,j)=dp(i-1,j-1)+1 while A[i]==B[j]d

转载 [转]矩阵奇异值分解(SVD)及其应用

机器学习中的数学(5)-强大的矩阵奇异值分解(SVD)及其应用版权声明：    本文由LeftNotEasy发布于http://leftnoteasy.cnblogs.com, 本文可以被全部的转载或者部分使用，但请注明出处，如果有问题，请联系wheeleast@gmail.com。也可以加我的微博: @leftnoteasy前言：    上一次写了关于PCA与LD

原创 Python in Windows -- build--''error: Unable to find vcvarsall.bat''

1) Install Microsoft Visual C++ Compiler for Python 2.7link is https://www.microsoft.com/en-gb/download/details.aspx?id=44266No need to install any other Visual Studio package2) Visual C++ compiler wou

翻译 Bash on Windows 重装，卸载

https://www.howtogeek.com/261188/how-to-uninstall-or-reinstall-windows-10s-ubuntu-bash-shell/卸载 cmd 或 powershell中输入 lxrun /uninstall 保留Ubuntu设置文件 输入 lxrun/uninstall /full 完全卸载 再次安装 输入lxrun /insta

转载 WEKA使用教程(经典教程转载)

WEKA使用教程目录 1. 简介2. 数据格式3.数据准备4. 关联规则（购物篮分析）5. 分类与回归6. 聚类分析 1. 简介 WEKA的全名是怀卡托智能分析环境（Waikato Environment for Knowledge Analysis），它的源代码可通过http://www.cs.waikato.ac.nz/ml/weka得到。同时weka也

转载 opencv中camshift例子的应用

贴一篇对opencv里camshiftdemo注释比较详细的博客原地址http://www.cnblogs.com/tornadomeet/archive/2012/03/15/2398769.html目标跟踪学习笔记_1(opencv中meanshift和camshift例子的应用)     在这一节中，主要讲目标跟踪的一个重要的算法Camshift，

原创 在iPhone6s系统邮箱添加QQ邮箱

在网上找了好多方法，都有一定的问题。下面是最简单的方法吧，可以接受和发送QQ邮件 网页进入QQ邮箱，设置-账户-POP3/IMAP/SMTP/Exchange/CardDAV/CalDAV服务，下面开启服务都开启 下面有一个方框：温馨提示：登录第三方客户端时，密码框请输入“授权码”   点击这里的生成授权码 然后复制，这个授权码就是在iPhone邮箱账户里的密码

原创 ML 逻辑回归 note

翻译 ML: Octave&Matlab Tutorial

https://share.coursera.org/wiki/index.php/ML:Octave_TutorialBasic Operations %% Change Octave prompt PS1('>> ');%% Change working directory in windows example:cd 'c:/path/to/desired/direc

原创 UVA 12307 旋转卡壳

http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=24011给定n个点3~1e5，求出包含所有点的最小行#include #include #include #include #include #include #include #include #include #include #include #

原创 hihoCoder 1231 Border Length(北京赛区(2015)网络赛 E)

求圆与多边形公共部分的周长http://blog.csdn.net/s_h_r/article/details/48623111由于除法根号用的有点多，精度被卡到1e-3才能过#include #include #include #include #include #include #include #include #include #inclu

原创 DLUTOJ 1330 GCD 【莫比乌斯反演+组合】

http://acm.dlut.edu.cn/problem.php?id=13301330: GCDTime Limit: 1 Sec  Memory Limit: 128 MBSubmit: 22  Solved: 6[Submit][Status][Web Board]DescriptionHow many non decreasing s

原创 zoj 3729 2013长沙regional 求斐波那契数列模n循环节长度

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=5075ArnoldTime Limit: 8 Seconds      Memory Limit: 65536 KBDo you know Vladimir Arnold? He's a mathematician who demonstrat

转载 根据日期判断星期几（使用基姆拉尔森计算公式）

//基姆拉尔森计算公式根据日期判断星期几void CalculateWeekDay(int y, int m,int d){ if(m==1||m==2) m+=12,y--; int iWeek = (d+2*m+3*(m+1)/5+y+y/4-y/100+y/400)%7; switch(iWeek){ case 0: printf("星期一\n"); break; c

原创 二分图带权最大独立集

问题：在一个二分图中，每个点都有一个权值，在保证选中的点不相交的情况下，求选出的点的最大权值和解法：【最小割】前提是能把图分成二分图，比如矩阵的格子点，同一侧的各个点不会有边相连。源点向X集合连流量为权值的边，Y集合各点向汇点连流量为权值的边。X与相连的Y连流量为INF的边。对于任意一个割集，与割相连的点都删掉就会得到一个独立点集。（割边只能是源点到X，和Y到汇点的边）因此

原创 网络流习题

一、POJ 2391    【拆点+二分+最大流】【题目大意】给定一个无向图，点 i 处有 Ai 头牛，点 i 处的牛棚能容纳 Bi 头牛，求一个最短时间 T 使得在 T 时间内所有的牛都能进到某一牛棚里去。 （1 1500, 0 【建模方法】将每个点 i 拆成两个点 i’, i’’，连边(s, i’, Ai), (i’’, t, Bi)。二分最短时间 T，若 d[i][j

原创 网络流

一、求最小割集http://paste.ubuntu.com/12211833/二、无源汇上下界可行流/*无源汇上下界可行流：边side(u,v)有流量上界up,下界low,构造的图使该边流量为up-low,u节点出流量low,v节点入流量low,用du[i]表示 i 节点入流之和与出流之和的差加上附加源点 S 汇点 T，如 du[i]>0，S->i 连容量为 du

原创 二分图最大匹配

讲解：http://www.renfei.org/blog/bipartite-matching.html（该博客最大独立集）补充定义和定理：最大匹配数：最大匹配的匹配边的数目最小点覆盖数：选取最少的点，使任意一条边至少有一个端点被选择最大独立数：选取最多的点，使任意所选两点均不相连最小路径覆盖数：对于一个 DAG（有向无环图），选取最少条路径，

原创 各种gcd

1.  求∑gcd(i, N) 1k=gcd(i,n)   , gcd(i/k,n/k)=1       ans=∑k*phi(n/k)[poj2480]2.求LCM（a,b,c,d,e,.......）依次分解质因数，保留最大次项3.求LCM（∑C(n,k)）  (1hdu 5407ans的因子：小于等于n+1的质数p，以及p的多次幂   最后

原创 hdu5365 判断四个点是否构成正方形

问题描述小花是一个热爱健身的姑娘，这天她下载了一个跑步软件，这个软件可以记录下小花跑步的轨迹。小花决定去公园跑步。公园里有许许多多的座椅，小花希望在一些座椅休息一下，并且她在两条座椅之间只跑直线。小花是一个完美主义者，她希望自己最后的轨迹是一个正三边形或者正四边形或者正五边形或者正六边形。小花会从某条座椅开始打开跑步软件，并在回到这个座椅后关闭。请问小花有多少种跑法。注：若两种跑法经过的座椅

原创 二分回顾

基础的东西很重要OvO1.STL里面有upper_bound()以及lower_bound() lower_bound(a,a+n,value);返回值为数组a中≥value的位置。如果都比value小，则指向a+n，注意该位置越界。 upper_bound(a,a+n,value);返回值为数组a中＞value的位置。如果都比value小，则指向a+n，注意该位置越界。区别在于如

原创 CF 552-C. Vanya and Scales

C. Vanya and Scalestime limit per test1 secondmemory limit per test256 megabytesinputstandard inputoutputstandard outputVanya has a scales for weighing loads and weights

原创 CF 551DGukiZ and Binary Operations 矩阵/位运算

D. GukiZ and Binary Operationstime limit per test1 secondmemory limit per test256 megabytesinputstandard inputoutputstandard outputWe all know

原创 codeforces/#305 Div1/547/C Mike and Foam 【容斥】

C. Mike and Foamtime limit per test2 secondsmemory limit per test256 megabytesinputstandard inputoutputstandard outputMike is a bartender at Rico's bar. At Rico

原创 hdu 5072 容斥

CoprimeTime Limit: 2000/1000 MS (Java/Others)    Memory Limit: 262144/262144 K (Java/Others)Total Submission(s): 1123    Accepted Submission(s): 453Problem DescriptionThere are n people st

原创 高效位运算 __builtin_系列函数

•int __builtin_ffs (unsigned int x) 返回x的最后一位1的是从后向前第几位，比如7368（1110011001000）返回4。 •int __builtin_clz (unsigned int x) 返回前导的0的个数。 •int __builtin_ctz (unsigned int x) 返回后面的0个个数，和__builtin_clz相对。 •in

原创 codeforces/#305 Div1/547/B Mike and Feet

B. Mike and Feettime limit per test1 secondmemory limit per test256 megabytesinputstandard inputoutputstandard outputMike is the president of country What-The-F

原创 LCM lightoj1215

Problem:  http://lightoj.com/volume_showproblem.php?problem=1215 LCM(a,b,c)=L give you a , b and L, please find the minmum c.solve: m=lcm(a,b), then lcm(m,c) is L m分解为p1^e1 * p2^e2 * p3^e3….. c分解

原创 大数a*b%m a^b%m

ll motil(ll a, ll b , ll m) //a*b%m{ ll ans=1; while(b){ if(b&1){ ans=(ans+a)%m; } b>>=1; a=(a<<1)%m; } return ans;}ll motil_pow(ll a, ll b , ll m)//a^b %m{ ll ans=1;

原创 一元线性同余方程组【数论

http://poj.org/problem?id=2891 数论P56#include "cstdio"#include <iostream>#include <cstring>#include <cmath>using namespace std;#define ll long longll exgcd(ll a,ll b,ll &x, ll & y){ if(b==0){

原创 【扩展欧几里得】练习题

1 .poj 1061青蛙的约会 最基础的一道 http://poj.org/problem?id=1061 题意：有两只青蛙，一只在坐标x，另一直在坐标y，青蛙x一次跳跃可以前进m单位距离，青蛙y一次跳跃可以前进n单位的距离，两青蛙都在同一纬度，该纬度长度为L。两只青蛙同方向同时跳啊跳，问你最少跳多少次，它们才可以相遇，如果不能相遇，输出impossble#include<cstdio>#

原创 summary of SET/MAP

set 一些归纳 set能够在O(logN)添加，删除，检查一个元素是否存在。 O(1)统计不一样的元素的个数,也就是set里元素的个数。 注意set每种元素只存放一次，也就是说s.insert(1)后，再s.insert(1)将不改变set的内容。写法： 头文件 set s;//可以是char，string，pair等等，要求可以比较 插入： s.insert(1); //插入元

原创 ZOJ 3865 Superbot

SuperbotTime Limit: 2 Seconds      Memory Limit: 65536 KBSuperbot is an interesting game which you need to control the robot on an N*M grid map.As you see, it's just a simple game: there i

转载 poj1837 Balance 天平和砝码

dp 背包http://blog.csdn.net/lyy289065406/article/details/6648094/题目大意：有一个天平，天平左右两边各有若干个钩子，总共有C个钩子，有G个钩码，求将钩码全部挂到钩子上使天平平衡的方法的总数。其中可以把天枰看做一个以x轴0点作为平衡点的横轴输入：2 4 //C 钩子数 与 G钩码数

原创 背包问题总结

一：01背包 题目描述：N件物品放入容量为V的包裹中，第i件物品体积Ci，价值Wi。怎么装可以使价值和最大。 方法：时间复杂度O(VN)，空间复杂度O(V) 定义F[v]为前i件物品放入容量v的背包中得到的最大价值。通过优化F[i,v]中仅保留F[v]即可。 代码memset(0,F,sizeof(F));for(int i=1;i<=N;i++){

原创 数论

欧拉函数Phi(x)等于不超过

原创 析取范式变主析取范式

#include #include #include #include using namespace std;const int n = 4; //n为变元个数int p[n];bool func1(){ return p[0]&&p[1]||p[2]&&p[3];}void change(){ for(int i=n-1;i>=0;i--){

原创 AC自动机模板 LA4670

出现次数最多的子串 n个小写字母组成的字符串和一个

原创 Trie字典树【模板

刘汝佳版本字典树#include #include #include using namespace std;const int maxnode = 4000 * 100 + 10;const int sigma_size = 26;// 字母表为全体小写字母的Triestruct Trie { int ch[maxnode][sigma_size]; int val