C++ com函数代码示例

int main()
{
	char s[11];
	long n;
	int i,j;
	while(scanf("%s",s)!=EOF)
	{
		for(i=0;s[i+1];i++) 
		{
			if(s[i]>=s[i+1]) break;
		}
		if(s[i+1]) printf("0\n");
		else
		{
			i=strlen(s);
			for(--i,n=0;i;i--) 
				n+=com(26,i);
			for(i=0;s[i];i++)
			{
				if(i==0)j='a';
				else j=s[i-1]+1;
				for(;j<s[i];j++)
				{
					n+=com('z'-j,strlen(s)-i-1);
				}
			}
			printf("%ld\n",n+1);
		}
	}
	return 0;
}

void dgBody::ApplyImpulsesAtPoint (dgInt32 count, dgInt32 strideInBytes, const dgFloat32* const impulseArray, const dgFloat32* const pointArray)
{
	dgInt32 stride = strideInBytes / sizeof (dgFloat32);

	dgMatrix inertia (CalculateInertiaMatrix());

	dgVector impulse (m_veloc.Scale3 (m_mass.m_w));
	dgVector angularImpulse (inertia.RotateVector (m_omega));

	dgVector com (m_globalCentreOfMass);
	for (dgInt32 i = 0; i < count; i ++) {
		dgInt32 index = i * stride;
		dgVector r (pointArray[index], pointArray[index + 1], pointArray[index + 2], dgFloat32 (0.0f));
		dgVector L (impulseArray[index], impulseArray[index + 1], impulseArray[index + 2], dgFloat32 (0.0f));
		dgVector Q ((r - com) * L);

		impulse += L;
		angularImpulse += Q;
	}

	dgMatrix invInertia (CalculateInvInertiaMatrix());
	m_veloc = impulse.Scale3(m_invMass.m_w);
	m_omega = invInertia.RotateVector(angularImpulse);

	m_sleeping	= false;
	m_equilibrium = false;
	Unfreeze ();
}

void div(int a[],int b[],int ans[])
{
    int i,j,k;
    int vvv[28][35]={0};
    int chengji[35]={0};
    int aans[35]={0};
    for(i=1;i<28;i++) vvv[i][i]=1;

    for(i=1;i<28;i++)
    {
        for(;;)
        {
            add(aans,vvv[i],aans);
            for(j=0;j<30;j++) chengji[j]=0;
            mul(aans,b,chengji);
//        cout<<i<<"  aans: ";coutt(aans);cout<<endl;
//            cout<<"chengji: ";coutt(chengji);cout<<endl;

            if(!com(a,chengji))
            {
                for(k=0;k<30;k++) ans[k]=aans[k];

            }
            else
            {
                sub(aans,vvv[i],aans);
                break;
            }
        }
    }
}

int main()
{
    int j,t,sc[11],i,k,co,cho,temp;
    int ans;
    scanf("%d",&t);
    while(t--)
    {
        for(i=0;i<11;i++)
            scanf("%d",&sc[i]);
        for(i=1;i<11;i++)
        {
            temp=sc[i];
            j=i-1;
            while((temp<sc[j])&&(j>=0))
            {
                sc[j+1]=sc[j];
                j=j-1;
            }
            sc[j+1]=temp;
        }
        scanf("%d",&k);
        co=0;
        for(i=0;i<11;i++)
            if(sc[i]==sc[11-k])
                co++;
        i=11-k;
        cho=0;
        while(i<11&&sc[i]==sc[11-k])
            {cho++;i++;}
        printf("%ld\n",com(co,cho));
    }
    return 0;
}

void KVINDRAOnlineDataAnalyser::HandleCommands(TString& ordre, TSocket* theSocket)
{
   // Gestion des ordres
   // Returns kTRUE if command was understood & executed
   // Returns kFALSE if command unknown
   //
   // STOP   will stop the analysis (calls EndRun() and EndAnalysis())
   // CLEAR  RAZ de tous les spectres
   //
   // DUMP=yes  start dumping fired data parameters
   // DUMP=no   stop dumping fired data parameters
   //
   // SAVE=toto.root          will save all histograms in '${ANALYSIS_RESULTS}/toto.root' (does not stop analysis)
   // SAVE=/home/toto.root    will save all histograms in '/home/toto.root' (does not stop analysis)
   // SAVE                    will save all histograms in '${ANALYSIS_RESULTS}/[class name]_Run[run number]_[date].root' (does not stop analysis)

   KVString com(ordre.Data());
   com.Begin("=");
   ordre = com.Next();
   ordre.ToUpper();
   if (ordre == "STOP") {
      // arrete l'analyse, en appelant EndRun() et EndAnalysis()
      fGoEventLoop = kFALSE;
      theSocket->Send("ok");
      return;
   } else if (ordre == "SAVE") {
      TString tmp_fil = com.Next();
      TString fil;
      bool send_name = kFALSE;
      if (tmp_fil.Contains("/")) fil = tmp_fil;
      else if (tmp_fil == "") {
         TDatime save_time;
         save_time.Set();
         TString date = save_time.AsString();
         date.ReplaceAll(" ", "_");
         fil = Form("${ANALYSIS_RESULTS}/%s_Run%d_%s.root", GetName(), GetRunNumber(), date.Data());
         send_name = kTRUE;
      } else fil = Form("${ANALYSIS_RESULTS}/%s", tmp_fil.Data());
      SaveSpectra(fil.Data());
      if (send_name) theSocket->Send(fil.Data());
      else theSocket->Send("ok");
      return;
   } else if (ordre == "CLEAR") {
      ClearAllHistos();
      theSocket->Send("ok");
      return;
   } else if (ordre == "DUMP") {
      TString tmp_fil = com.Next();
      tmp_fil.ToLower();
      if (tmp_fil == "yes")fDumpEvents = kTRUE;
      else if (tmp_fil == "no")fDumpEvents = kFALSE;
      else {
         theSocket->Send("unknown command");
         return;
      }
      theSocket->Send("ok");
      return;
   }
   theSocket->Send("unknown command");
}

int com(int n, int r)
{
	if(r==0 || r==1)
	return n;

	return ((n-r+1)*com(n,r-1)/r);
}

/*
 * prepares the communication maps for remote elements
 * should be called immediately after load balancing,
 * before any work transfer.
 *
 */
void
NonlocalMaterialWTP :: init(Domain *domain)
{
    int ie, gie, nelem = domain->giveNumberOfElements();
    EngngModel *emodel = domain->giveEngngModel();
    Element *elem;
    int nproc = emodel->giveNumberOfProcesses();
    int myrank = emodel->giveRank();
    CommunicatorBuff cb(nproc, CBT_dynamic);
    Communicator com(emodel, &cb, myrank, nproc, CommMode_Dynamic);
    this->nonlocElementDependencyMap.clear();

    // build nonlocal element dependency array for each element
    for ( ie = 1; ie <= nelem; ie++ ) {
        elem = domain->giveElement(ie);
        if ( ( elem->giveParallelMode() == Element_local ) ) {
            gie = elem->giveGlobalNumber();
            this->giveElementNonlocalDepArry(nonlocElementDependencyMap [ gie ], domain, ie);
        }
    }

    /* send and receive nonlocElementDependencyArry of migrating elements to remote partition */
    com.packAllData(this, domain, & NonlocalMaterialWTP :: packMigratingElementDependencies);
    com.initExchange(MIGRATE_NONLOCALDEP_TAG);
    com.unpackAllData(this, domain, & NonlocalMaterialWTP :: unpackMigratingElementDependencies);
    com.finishExchange();
}

void RenderCenterOfMass (NewtonWorld* const world)
{
	glDisable (GL_LIGHTING);
	glDisable(GL_TEXTURE_2D);

	glColor3f(0.0f, 0.0f, 1.0f);

	glBegin(GL_LINES);
	for (NewtonBody* body = NewtonWorldGetFirstBody(world); body; body = NewtonWorldGetNextBody(world, body)) {
		dMatrix matrix;
		dVector com(0.0f); 
		
		NewtonBodyGetCentreOfMass (body, &com[0]);
		NewtonBodyGetMatrix (body, &matrix[0][0]);

		dVector o (matrix.TransformVector (com));

		dVector x (o + matrix.RotateVector (dVector (1.0f, 0.0f, 0.0f, 0.0f)));
		glColor3f (1.0f, 0.0f, 0.0f);
		glVertex3f (o.m_x, o.m_y, o.m_z);
		glVertex3f (x.m_x, x.m_y, x.m_z);

		dVector y (o + matrix.RotateVector (dVector (0.0f, 1.0f, 0.0f, 0.0f)));
		glColor3f (0.0f, 1.0f, 0.0f);
		glVertex3f (o.m_x, o.m_y, o.m_z);
		glVertex3f (y.m_x, y.m_y, y.m_z);

		dVector z (o + matrix.RotateVector (dVector (0.0f, 0.0f, 1.0f, 0.0f)));
		glColor3f (0.0f, 0.0f, 1.0f);
		glVertex3f (o.m_x, o.m_y, o.m_z);
		glVertex3f (z.m_x, z.m_y, z.m_z);

	}
	glEnd();
}

void  momentofinertiamatrix(int type, double *MofImatrix)
{
  int i,j,count;
  double p1[3],massi;
  std::vector<double> com (3);
  count=0;
  updatePartCfg(WITHOUT_BONDS);
  for(i=0;i<9;i++) MofImatrix[i]=0.;
  com = centerofmass(type);
  for (j=0; j<n_part; j++) {
    if (type == partCfg[j].p.type) {
      count ++;
      for (i=0; i<3; i++) {
      	p1[i] = partCfg[j].r.p[i] - com[i];
      }
      massi= (partCfg[j]).p.mass;
      MofImatrix[0] += massi * (p1[1] * p1[1] + p1[2] * p1[2]) ; 
      MofImatrix[4] += massi * (p1[0] * p1[0] + p1[2] * p1[2]);
      MofImatrix[8] += massi * (p1[0] * p1[0] + p1[1] * p1[1]);
      MofImatrix[1] -= massi * (p1[0] * p1[1]);
      MofImatrix[2] -= massi * (p1[0] * p1[2]); 
      MofImatrix[5] -= massi * (p1[1] * p1[2]);
    }
  }
  /* use symmetry */
  MofImatrix[3] = MofImatrix[1]; 
  MofImatrix[6] = MofImatrix[2]; 
  MofImatrix[7] = MofImatrix[5];
  return;
}

void sub(int a[],int b[],int ans[])
{
    int jinwei=0;
    int i,j,k;
    if(com(a,b))
    {
        ans[35]=1;
        sub(b,a,ans);
        return ;
    }
    for(i=29;i>=0;i--)
    {
        ans[i]=jinwei+a[i];
        if(ans[i]<b[i])
        {
            ans[i]+=10;
            ans[i]-=b[i];
            jinwei=-1;
        }
        else
        {
            ans[i]-=b[i];
            jinwei=0;
        }
    }
}

int main(){
 char buf[1024];
 while(scanf("%s",buf) != EOF){
     printf("%d\n",com(buf));
 }

 return 0;
}

const arma::vec& CylindricalBirefringentMaterial::centerOfMass()
{
  com.set_size(3);
  com.fill(0.0);
  double mass = 0.0;
  for ( unsigned int z=0;z<Nz();z++ )
  for ( unsigned int y=0;y<Ny();y++ )
  for ( unsigned int x=0;x<Nx();x++ )
  {
    com(0) += x*get(x,y,z);
    com(1) += y*get(x,y,z);
    com(2) += z*get(x,y,z);
    mass += get(x,y,z);
  }
  com /= mass;
  return com;
}

 Vector3d Skeleton::getWorldCOM() {
     assert(mMass != 0);
     Vector3d com(0, 0, 0);
     const int nNodes = getNumNodes();
     for(int i = 0; i < nNodes; i++) {
         BodyNode* node = getNode(i);
         com += (node->getMass() * node->getWorldCOM());
     }
     return com / mMass;
 }

 vector<vector<int> > combine(int n, int k) {
     // Start typing your C/C++ solution below
     // DO NOT write int main() function
     vector<vector<int> > ret;
     if(n < k || n <= 0 || k <= 0)
       return ret;
     vector<int> com(k); 
     combine(n, k, 0, com, ret);
     return ret;
 }

void loop()
	{
	float result;
	int key=resist();
	if(key==15)
		{
		//LCM CLEAR
		wlcmd_4('c',0x01);
		_delay(1600);
		//student num
		student_num();
		//parameter reset
		aa=0;bb=0;_status0=0;
		lcm_num[0]=0;
		lcm_num[1]=0;
		lcm_num[2]=0;
		lcm_num[3]=0;
		lcm_num[4]=0;
		lcm_num[5]=0;
		lcm_num[6]=0;
		return;
		}
	if((key>=0)&&(key<=9)&&(_status0!=4))
		{
			wlcmd_4('d',key);
			if(_status0==0)
			{
				aa=aa*10+key;
			}
			if(_status0==2)
			{
				bb=bb*10+key;
			}
		return;	
		}
	if((key>=10)&&(key<=13)&&(_status0!=4))
		{
			if(_status0==0)
			{
				wlcmd_4('d',key);
				operand=key;
				_status0=2;
				return;
			}
		}
	if((key==14)&&(_status0!=4))	
		{
			wlcmd_4('d',key);
			result=com(aa,bb,operand);
			divide(result);
			_status0=4;
			return;	
		}
	return;
	}