2018年01月_费马最后定理

原创 MATLAB-图像处理

s={};for i=1:47 name11=num2str(i); name12='.png'; % name11=strcat('D:\111\',name11); name21=strcat('D:\112\',name11); name22='.txt'; name1=strcat(name11,name12); na

2018-01-31 12:49:10 262

原创 Python-文件夹下所有文件批处理重命名，修改格式

#coding:utf8import os;def rename(): i=0 path=(r'C:\Users\25459\AppData\Local\Programs\Python\Python37\ceshi');#一定要加r转义单引号 filelist=os.listdir(path)#该文件夹下所有的文件（包括文件夹）

2018-01-31 09:35:43 1744

原创 MATLAB-字符串连接（以循环变量为名创建文本）

for i=1:3 d=strcat(num2str(i),'abc')t=fopen(d,'w');fprintf(t,d);end

2018-01-30 23:30:14 6303

原创 MATLAB-txt读取,输出txt（以矩阵形式）

fid=fopen('D:\新建文件夹\数据\0.txt','w');%写入文件路径a=int8(dd);[m,n]=size(a); for i=1:1:m for j=1:1:n if j==n fprintf(fid,'%g\r\n',a(i,j)); else fprintf(fid,'%g',a(i

2018-01-30 16:22:36 12967 1

原创 MATLAB-支持向量机-非线性-p213

非线性内核（首先使用线性内核发现判错率比较高）%原始数据cancerdata.txt可在网上下载，数据中的B替换成1，M替换成-1，X替换成2，删除了分割符*,替换后的数据命名成cancerdata2.txtclc,cleara=load('cancerdata2.txt');a(:,1)=[]; %删除第一列病例号gind=find(a(:,1)==1);

2018-01-29 20:07:20 794

原创 MATLAB-支持向量机模板

注意每个数据项为一行不是一列与下图表格一致就行 b0与dd0数据的分量即列根据问题选择调整group为每项已知数据的类别clc, cleara0=load('fenlei.txt'); %把表中x1...x8的所有数据保存在纯文本文件fenlei.txt中a=a0'; b0=a(:,[1:27]);dd0=a(:,[28:end]); %提取已分类和待分类的数据[b,ps

2018-01-29 19:32:21 1070

原创 C++-并查集模板

并查集的模板 find()函数的定义连接的操作需要注意的是初始化和p[y]=x的顺序不能颠倒!#include#includeusing namespace std;const int maxn=1000+10;int a[maxn]; //测试数组 int p[maxn]; //使用时要初始化 //find模板 int find(int x){ return x==p[x]?

2018-01-28 21:00:57 939

原创 C++-图论-边的结构体模板(带权可排序，直接赋值)

//图论常用的边的结构体模板这里用平面上的点和距离作例子 #include#include#include#includeusing namespace std;const int maxn=1000+10;int x[maxn],y[maxn]; //边的模板 struct Edge{ int u,v,d; //u，v为两个分量 d为权值结构体的大小由权值排序 Ed

2018-01-28 20:39:15 488

原创 11.2.2-uva1151-kruskal算法综合运用，贪心

#include#include#include#includeusing namespace std;const int maxn=1000+10;int p[maxn],x[maxn],y[maxn]; //x于y表示点的坐标 int n;//点的个数 vector subn[10];//用vector作容器存储套餐数组每个数组包含该套餐所包含的点 int cost[10

2018-01-28 20:26:21 236

原创 11.3.4-uva10048-floyd变形-两点之间路径上最大边的最小值

思路看主算法#include#include#includeusing namespace std;const int maxn=100+5;int d[maxn][maxn];const int INF=1000000000;int main(){ int n,m,Q,kase=0; while(cin>>n>>m>>Q&&n){ //初始化 for(int i=

2018-01-28 16:43:47 291

原创 C++-ID()函数模板-vector,map

以vector为容器每次放入新元素时在容器中进行循环寻找是否有相同元素如果有返回相应编号没有则放入容器创建新编号时间复杂度较大每次查询复杂度O（n^2）#include#include#include#includeusing namespace std;vector names;int n;int ID(string s){ for(int i=0;i<names.s

2018-01-28 15:50:41 469

原创 MATLAB-元胞自动机-森林火灾

%CA driver%%forest fireclfclear alln=100;Plightning = .001;Pgrowth = .0001; z=zeros(n,n);o=ones(n,n);veg=2*o;sum=z;imh = image(cat(3,z,veg==2,z));set(imh, 'erasemode', 'none')axis equal

2018-01-27 17:13:44 11225 2

原创 MATLAB-元胞自动机-激发介质（ BZ reaction or heart）

%CA driver%excitable mediaclf; clear alln=128;z=zeros(n,n);cells=z;cells = (rand(n,n))<.1 ;%cells(n/2,n*.25:n*.75) = 1;%cells(n*.25:n*.75,n/2) = 1;sum=z;imh = image(cat(3,cells,z,z));set(im

2018-01-27 15:33:42 911

原创 MATLAB-常用函数-r2016a

5.2.8-171-交互式绘图-[x,y]=ginput(n) gtext[arg]cat() 构造多维数组>> c = cat(3,a,b) c(:,:,1) = 8 1 6 3 5 7 4 9 2 c(:,:,2) = 1 1 1 1

2018-01-27 11:37:38 456

原创 MATLAB-元胞自动机-生命游戏-GUI

%Conway's life with GUIclfclear all%===========================%build the GUI%define the plot buttonplotbutton=uicontrol('style','pushbutton',...'string','Run', ...'fontsize',12, ...'position

2018-01-27 11:13:59 1029

原创 MATLAB-在线雇佣-仿真模拟法

a=[];s=10000;for t=1:99 k=0;for i=1:s x=randperm(100); best=0; best1=0; for j=1:t if x(j)>best best=x(j); end end for j=t+1:100

2018-01-26 21:10:32 1317 1

原创 MATLAB-sparse()与full()-稀疏矩阵

sparse函数功能：创建稀疏矩阵用法1：S=sparse(X)—将矩阵X转化为稀疏矩阵的形式，即矩阵X中任何零元素去除，非零元素及其下标（索引）组成矩阵S。如果X本身是稀疏的，sparse(X)返回S。举例如下：>> a=[1,0,2;0,0,1;0,0,6];>> aa = 1 0 2 0 0 1

2018-01-24 17:13:03 447

原创 8.6-247-uva12174-Shuffle-滑动窗口

#include#include#includeusing namespace std;const int maxn=100000+5;int s,n,x[maxn*3],cnt[maxn],ok[maxn*2];int main(){ int T; cin>>T; while(T--){ int n,s; cin>>s>>n; fill(x,x+n+2*s,-1);

2018-01-24 16:49:55 169

原创 8.5-246-uva11093-just finish it up-模拟法

#include#include#include#includeconst int maxn=100000+5;int p[maxn],q[maxn],n;using namespace std;int go(int s){ int fuel=p[s]-q[s]; for(int i=(s+1)%n;i!=s;i=(i+1)%n){ if(fuel<0)

2018-01-23 17:01:52 397

原创 8.5-239-uva11573- Unique snowflakes-滑动窗口

a[l~r]是可行解则a[l+1~r]也是可行解 r不能增大时，只需将r右滑一位，不必减少r，即r在整个过程只增不减每次要么增加r要么增加l 故增加次数为O(n)判断是否出现用set()#include#include#include#includeconst int maxn=1000000+5;int a[maxn];using namespace std;int m

2018-01-23 16:04:39 143

原创修改圈近似算法

function maina=zeros(6);a(1,2)=56;a(1,3)=35;a(1,4)=21;a(1,5)=51;a(1,6)=60;a(2,3)=21;a(2,4)=57;a(2,5)=78;a(2,6)=70; a(3,4)=36;a(3,5)=68;a(3,6)=68; a(4,5)=51;a(4,6)=61;a(5,6)=13; a=a+a'; L=size(a,1)

2018-01-23 14:46:41 2457 9

原创模拟退火

clc, clearsj0=load('sj.txt'); %加载100个目标的数据，数据按照表格中的位置保存在纯文本文件sj.txt中x=sj0(:,[1:2:8]);x=x(:);y=sj0(:,[2:2:8]);y=y(:);sj=[x y]; d1=[70,40]; sj=[d1;sj;d1]; sj=sj*pi/180; %角度化成弧度%数据预处理这里sj为n*2矩阵

2018-01-23 11:38:26 315

原创 8.5-237-uva11134-Fabled Rooks-问题分解-贪心

行和列是无关的，当且仅当行列都满足条件，则存在解不失一般性现在只考虑x坐标（x为列的编号）题目等价于找这样的一个序偶满足：用一系列闭区间[xli,xri]覆盖整数点1,2,3.....n 且每个整数对应于唯一的一个区间（区间与整数点是一对一的）故可用贪心法从第一个整数点1开始对闭区间进行遍历找到覆盖1的区间中右端点最小的一个（找最小的b与找最小值的方法类似，即先假设最小，再循环

2018-01-22 19:02:10 164

原创 8.5-237-uva1152-中途相遇法-二分法

合为0的四个值给定4个n元素集合a,b,c,d，分别取一个元素出来使其和为0，有几种取法？思路：直接枚举需要四个循环时间复杂度为O(n^4);现在考虑中途相遇法用两个循环遍历a+b 放入另一个数组sums1[]中并排序循环时间复杂度为O（n^2）对sums[]中的n*n个元素排序为n^2*log(n^2) 总的复杂度为O（n^2*logn）同理循环遍历c+d并用二分法查找-c-d 时间复...

2018-01-22 15:00:34 202

原创 8.6-244-uva714-Copying books-二分法-最大值尽量小

题目：把一个包含m个正整数的序列划分（划分表面不相交）成k个非空连续子序列设第i个序列的和为S（i）找一个划分使最大的S（i）（记为maxp）的尽量小多种方案优先取S（1）最小，，S（1）相同S（2）最小，，，依此类推。思路：记P（x）为将序列划分为m个，maxp=x，满足m最终找到的x即为maxpk表示对应需要‘/’的最小个数对于5个元素的数组 1 1 1 1

2018-01-21 19:10:42 423

原创 0-1背包

#include#includeusing namespace std;int v[]={0,1,2,1,1};int w[]={0,1,2,3,4};int c=3;int d[5][4]={0};int main(){ for(int j=0;j<=c;j++) d[0][j]=0; //边界初始化 //下面两个for循环可调换位置 for(int i=1;i<=4

2018-01-19 21:42:27 130

qq_41093189的博客